{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "480765be",
   "metadata": {},
   "outputs": [],
   "source": [
    "from astropy.io import fits\n",
    "import HAWCtoHerschel as HtoH"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "d5477674",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Source information\n",
    "SourceName = 'Fil10'\n",
    "FullName = 'Filament 10'\n",
    "# Loading the fits files\n",
    "SourceT = fits.open('../FIELDMAPS_Plots/'+SourceName+'/Data/'+SourceName+'_Tdust_10_10.fits')\n",
    "SourceN = fits.open('../FIELDMAPS_Plots/'+SourceName+'/Data/'+SourceName+'_N_10_10.fits')\n",
    "SourcePol = fits.open('../FIELDMAPS_Plots/'+SourceName+'/Data/'+SourceName+'_Gal.fits')\n",
    "SourcePlanck = fits.open('../FIELDMAPS_Plots/'+SourceName+'/Data/'+SourceName+'_Planck.fits')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "8fbf97a2",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Plot parameters\n",
    "ScaleLength = 0.0197572 # 0.0179049 # 1 pc in degrees for the target\n",
    "\n",
    "# Region location, width, and height\n",
    "l_center = 332.4266828 # Position in galactic longitude (degrees)\n",
    "b_center = -0.0840824 # Position in galactic latitude (degrees)\n",
    "l_width = 0.6019563 # Width in galactic longitude (degrees)\n",
    "b_height = 0.2932983 # Height in galactic latitude (degrees)\n",
    "Region = [l_center,b_center,l_width,b_height]\n",
    "\n",
    "# Intensity contours for the plot\n",
    "Icontours = [20, 40, 60, 80, 100] # in mJy/arcsec^2\n",
    "# Ranges for plots\n",
    "Nscale = [0,2.0e+22] # Range for column density in cm^-2\n",
    "Tscale = [15.0,30.0] # Range for temperature map in K\n",
    "Iscale = [0.0,50.0] # Range for intensity map in mJy/arcsec^2\n",
    "# Vector scaling fraction in polarization plots\n",
    "Pvscale = 0.6\n",
    "\n",
    "# Figure size\n",
    "MapFigSize = [10.67,6] # in inches\n",
    "\n",
    "# Effective HAWC+ beam size to use\n",
    "HAWCBeam = 0.00519 # 18.7'' in Band E"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "a2a163b6",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Beam Size of polarization data\n",
      "18.2\n",
      "Beam size of target pojection\n",
      "36.395999999999994\n",
      "Pixel size of the polarization data\n",
      "4.550000000000001\n",
      "Pixel size of the target projection\n",
      "14.0000004\n",
      "Target standard deviation for smoothing in arcseconds\n",
      "13.384761981071023\n",
      "Target standard deviation for smoothing in pixels\n",
      "2.941705929905719\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "WARNING: nan_treatment='interpolate', however, NaN values detected post convolution. A contiguous region of NaN values, larger than the kernel size, are present in the input array. Increase the kernel size to avoid this. [astropy.convolution.convolve]\n",
      "c:\\Users\\simon\\OneDrive\\Python\\FIELDMAPS_Testing_Herschel\\HAWCtoHerschel.py:96: RuntimeWarning: invalid value encountered in divide\n",
      "  dPI = ((StokesQ*dQ)**2.0 + (StokesU*dU)**2.0)**0.5/PI_biased\n",
      "c:\\Users\\simon\\OneDrive\\Python\\FIELDMAPS_Testing_Herschel\\HAWCtoHerschel.py:106: RuntimeWarning: invalid value encountered in divide\n",
      "  P_biased = 100.0*PI_biased/StokesI\n",
      "c:\\Users\\simon\\OneDrive\\Python\\FIELDMAPS_Testing_Herschel\\HAWCtoHerschel.py:109: RuntimeWarning: invalid value encountered in divide\n",
      "  dP = np.absolute(P_biased*((dPI/PI_biased)**2.0 + (dI/StokesI)**2.0)**0.5)\n",
      "c:\\Users\\simon\\OneDrive\\Python\\FIELDMAPS_Testing_Herschel\\HAWCtoHerschel.py:120: RuntimeWarning: invalid value encountered in divide\n",
      "  dO = (0.5*180.0/math.pi)*((StokesQ*dU)**2.0 +\n"
     ]
    }
   ],
   "source": [
    "# Reprojecting the HAWC+ Stokes data to the Herschel pixel scale\n",
    "NewPol = HtoH.ReprojectCube(SourcePol,SourceN, Smooth=True)\n",
    "#NewPol.writeto('test2.fits', overwrite=True)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "31eb9087",
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "c:\\Users\\simon\\OneDrive\\Python\\FIELDMAPS_Testing_Herschel\\HAWCtoHerschel.py:228: RuntimeWarning: invalid value encountered in divide\n",
      "  imask_01 = np.where(Pdata[0].data/Pdata[1].data < IdI) # Total intensity SNR threshold\n"
     ]
    }
   ],
   "source": [
    "# Creating downsampled vector catalogs (independent vectors)\n",
    "# Creating the mask for the region\n",
    "DownsampleMask = HtoH.DownsampleVectors(NewPol,Step=2)\n",
    "# Creating the new catalogs for the region\n",
    "CatI, CatdI, CatQ, CatdQ, CatU, CatdU, CatP, CatdP, CatO, CatB,CatdO, CatPI, CatdPI, CatN, CatT = HtoH.MakePCats(SourceN, SourceT, SourcePol, NewPol, CatMask=DownsampleMask)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "2c4fb0cb",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "===\n",
      "Curve fit test\n",
      "Log scale fit: (-0.6869272610108983,43.82548180212998)\n",
      "Direct power law fit: (-0.9000362632107898,82.68448839872504)\n",
      "===\n",
      "\n",
      "Power law index: -0.6869272610108983 ± 0.03044164162002954\n",
      "Coefficient: 43.82548180212998 ± 4.695540881055554\n",
      "Chi-Squared: 1202.298583289453\n",
      "Number of elements: 187\n",
      "Reduced Chi-Squared: 6.4989112610240705\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAfUAAAF2CAYAAAB3bAFnAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/P9b71AAAACXBIWXMAAA9hAAAPYQGoP6dpAAB54ElEQVR4nO3deXhMZ/vA8e9MViEJEdlE7Pu+V2n5oVQX2qp6dUEStBp7aXmrFG0tb2tXSiO0pZa2dKO0XktftRNLEUUQQiLIIiSRzPn9MZ1pJplJZiaTzCTuz3WdK5Mz5zznnsnJ3HOe8ywqRVEUhBBCCFHqqe0dgBBCCCFsQ5K6EEIIUUZIUhdCCCHKCEnqQgghRBkhSV0IIYQoIySpCyGEEGWEJHUhhBCijJCkLoQQQpQRktSFEEKIMuKhSeqZmZm8//77ZGZm2jsU4WDk3BCmyLkhjHHk88IhkvrSpUtp1qwZXl5eeHl50aFDB7Zu3ap/PiMjg4iICCpXrkyFChXo27cvCQkJFh0jMzOTadOmOeQfQdiXnBvCFDk3hDGOfF44RFIPDg5m1qxZHDlyhMOHD9O1a1f69OnDn3/+CcDYsWP58ccf2bhxI7t37yY+Pp4XXnjBzlHb3pIlSxy+fGvKsGQfc7YtbJuCni/u97i4yLlR9HOjLJ4XIOeGLbYrU+eG4qAqVaqkfP7550pycrLi4uKibNy4Uf/cmTNnFEDZt2+f2eWlpKQogJKSklIc4dpEw4YNHb58a8qwZB9zti1sm4KeN/acnBsPx7lhzXNybjj+uWGL7crSZ4azXb9RGJGTk8PGjRtJT0+nQ4cOHDlyhAcPHtC9e3f9Ng0aNCAkJIR9+/bxyCOPGC0nMzPToGokOTkZgJSUlGKNvyhycnJITU116PKtKcOSfczZtrBtCnre2HO634vzvS8qOTeKfm5Y85ycG45/bthiO0f7zFAUhVu3buHj44Na/U+FupubG25uboXu7BBOnDihlC9fXnFyclK8vb2Vn3/+WVEURVmzZo3i6uqab/u2bdsqb7/9tsnypk6dqgCyyCKLLLLIUiaWqVOnFppLHeZKvX79+kRHR5OSksI333zDoEGD2L17t9XlTZo0iXHjxul/T0lJISQkhLi4OLy8vGwRsmW+/BJGjIBWrWDnzpI/vg1lZ2dz7do1qlatCqB/7Oyc/3TKva2zs3O+34ty7Lz7F7VsIYRwBKmpqVSrVo0rV67g7e2tX1/oVTrgMJ98rq6u1KlTB4DWrVtz6NAhFixYQP/+/cnKyiI5OZmKFSvqt09ISCAgIMBkeaaqKXQt7EvcCy9ok/rRo3D/Pvj7l3wMNpKdnU1qaqr+fdQ91iVSjUajX5f7sS6p593e2mMbS+pFKVsIIRyJt7e3xfnKIVq/G6PRaMjMzKR169a4uLiwY8cO/XMxMTFcuXKFDh062DFCCwUGaq/SAXJ11yuLNBoNycnJaDQae4cihBAPFYe4nJk0aRK9evUiJCSEtLQ01q5dy65du9i2bRve3t6Eh4czbtw4fHx88PLyYuTIkXTo0MFkIzmH9fTT2iv1LVtg8GB7RyOEEKKMcYiknpiYyMCBA7l+/Tre3t40a9aMbdu28cQTTwAwb9481Go1ffv2JTMzk549e/Lpp5/aOWorPP00zJgB27bBgwfg4mLviIQQQpQhDpHUIyMjC3ze3d2dJUuWlL5BAPJq0wZ8fSEpCfbuhS5d7B2REEKIMsRh76mXSU5O8NRT2sc//GDfWIQQQpQ5ktRLWp8+2p/ffw+KYt9YhBBClCmS1Etajx7g5gYXL8LfY9sLIYQQtiBJvaRVqAC6IW+//96+sQghhChTJKnbQ+4qeCGEEMJGJKnbw7PPan8eOgTx8faNRQghRJkhSd0eAgKgfXvtY2kFL4QQwkYkqduLVMELIYSwMUnq9qJL6v/9L6Sl2TcWIYQQZYIkdXtp2BDq1IGsLPjlF3tHI4QQogyQpG4vKpVUwQshhLApSer2pEvqP/+sneBFCCGEKAJJ6vb06KPaCV6Sk+H33+0djRBCiFJOkro9OTnBM89oH0sVvBBCiCKSpG5vMsGLEEIIG5Gkbm9PPAHu7nD5Mpw4Ye9ohBBClGKS1O2tfHltYgepghdCCFEkktQdga4KXoaMFUKUkFWrVqFSqVCpVPp1Xbp0QaVSMXjwYPsFJopEkrojeOYZbb/1I0fg6lV7RyOEKOV0ydnYsnnzZgCqVKlC+/btaa+bh6IMGDx4MCqVii5dupi1/ciRI2nevDnOzs6oVCoCAgKMbnfkyBGefPJJvLy88PDwoFOnTvz22282jNx2nO0dgAD8/aFDB/jjD20VfESEvSMSQpQBrq6utGzZ0mCdj48PAE8//TRPP/20PcJyGF9++SWurq74+Phw8+ZNo9ucOHGCxx9/nHv37uHr64uXlxd79+7lySefZMuWLfTo0aOEoy6YXKk7iuef1/789lv7xmGG+Ph4FixYQLxMGyuEQwsMDGT//v0Gy+OPPw4Yr343ZtCgQdStWxdPT09cXV2pXr06o0aNIjU1Vb+N7gq5Ro0arF69mmrVquHl5cWYMWO4f/8+Y8aMwdvbm5CQEJYtW2ZQfnx8PGFhYQQFBeHq6kqtWrWYMWMG2dnZ+m10NQ8DBw5k6tSpBAYGUqlSJV599VXS/p47Q3dsgN27d+tf265du0y+tpMnT5KYmMhTTz1lcpvJkydz7949atSowcWLF7l06RLt27cnJyeH8ePHF/je2YMkdUfRt6/25+7dkJho31gKERkZydmzZ4mMjLR3KEKIYvb9999z584dateuTbVq1bhy5QqLFi0iPDw837bXr1/nzTffxM3NjbS0NBYsWEDr1q2JiorC29ubuLg4IiIiOHPmDAC3bt3ikUceISoqirt379KwYUPi4uKYMmUKw4YNy1f+unXrmDdvHuXKlSM5OZk1a9Ywa9YsAFq2bImvry8Anp6e+lsLXl5eJl9btWrVCnzt2dnZ+mr2Hj164OnpibOzM7179wa0Xwoc7eJGkrqjqFkTWrcGjQY2bbJ3NAUKDw+nQYMGRv+phRCO4/Lly/nuqVtq9+7dJCUlER0dzYULF3j33XcB2Lx5MxkZGQbbZmVlsX37ds6dO6dPmOfPn+f48eOcPn0ad3d3NBoNu3fvBmDx4sXExcXh7+/PhQsXOH78ON988w2grUk4f/68Qfnu7u6cOXOG8+fP07p1awB27NgBwKZNm/S3E1q1aqWvmWjVqpXFr1knKSmJ+/fvA+Dn56df7+/vr3985coVq8svDnJP3ZH066dtLPfNN/D66/aOxqSgoCBGjx5NUFCQvUMRQhTA2D11S/3222+88sorXLhwwSCJZ2dnc/PmTYOr3UqVKtGxY0cAQkJCiIuLo0mTJtSoUQPQJsYrV66QkJAAwMGDBwFISEgwSJoAiqJw4MAB6tSpo1/XtWtXqlatCkCDBg04cuSIvqySpDjwQGGS1B3Jiy/CxImwcyckJWnHhS9jdPfjx44dS0hIiL3DEaJM091Tt9aaNWv0940DAwOpVq0aSUlJXLx4EYCcnByD7XNXdTs7O+dbp6spyJsUPT09adSoUb7je3h4GPxesWLFfOUXZ4L19fWlXLly3L9/n8Rct0VzP3a0zzGpfncktWtDy5aQkwN/dzspa+R+vBClh+4LgaenJ7GxsRw4cMCmrb3btm0LaBP0unXr9FXmv/76K2+++SbP6xoQm0n3JSA9Pd0m8Tk7O9OtWzcAtm/fTlpaGtnZ2fzw95giTZs2dbgaS0nqjqZfP+3PjRvtG4eZLG0JL/fjhSg9mjVrBkBaWhq1atWiVq1abNiwwWblR0REULVqVe7cuUP9+vVp0aIFtWvXpnLlygwaNMji8ho0aADA4cOHadq0KY888oj+nrgxXbp0oU6dOnz33XeA9h56nTp1qFOnDgcOHADggw8+oFy5cly6dIlatWpRo0YNDhw4gJOTE3PmzLHiVRcvSeqO5sUXtT937IBbt+wbixksvfKW+/FClB7h4eGMGzcOX19f0tLS6NKlC9OnT7dZ+VWqVGH//v2EhoZSuXJl/vzzT+7fv89jjz3GvHnzLC4vLCyMvn374u3tzalTpzhw4EC+WwS5Xbp0iQsXLui7xeXk5HDhwgUuXLig/zLQvHlzdu/ezRNPPEFGRga3bt3i0UcfZcuWLTz55JPWvfBipFIc+Y6/DaWmpuLt7U1KSkqBXRwcQosWcPw4REZCWJi9o8knOzubq1evEhwcTHx8PPPmzTO4R577eUD/2NnZ2eA53T0xa4+dd/+ili2EEI6gKPlKrtQdke5qvRRUwee+8tZoNCQnJ6PRaOwdlhBCPJQkqTsi3X31336DO3fsG4sFJKkLIYR9SVJ3RPXrQ9OmkJ0t07HmIUPUCiGEaZLUHVUpqoIviK2TsHSJE0II0ySpOypdFfyvv0Jysl1DKQpbJ2HpEieEEKZJUndUDRtCo0bw4AH8PdBBaWTrJCxd4oQQwjRJ6o5Md7X+9wQHpZEkYSGEKDmS1B2Z7r76tm2QkmLfWIQQQjg8SeqOrHFjbRV8VpbDT8cqhBDC/hwiqc+cOZO2bdvi6emJn58fzz33HDExMQbbdOnSJd+8wG+88YadIi4hKhW8/LL28dq19o1FCPFQWLJkCTVq1MDd3Z327dvrp0ctyLVr13j11VepXLky5cqVo2nTphw+fFj/fFpaGmPGjKF69eqUK1eORx99lEOHDjlM/JbsM2vWLFQqFWPGjLFh1LbjEEl99+7dRERE6GfnefDgAT169Mg3087QoUO5fv26fnHEwfRtbsAA7c8dO+DGDfvGYkMyUI0Qjmf9+vWMGzeOqVOncvToUZo3b07Pnj0NphrN686dO3Ts2BEXFxe2bt3K6dOn+eSTT6hUqZJ+myFDhvDrr7/y5ZdfcvLkSXr06EH37t25du1aoTF16dKFVatWFVv8luxz6NAhPvvsM/1ENw5JcUCJiYkKoOzevVu/rnPnzsro0aOtLjMlJUUBlJSUFBtEWMIeeURRQFHmz7d3JIqiKMqDBw+U2NhY5cGDB0Yf379/3+jzuffNvY21x7bkOSHs6cCBA0rnzp0Vd3d3pX79+sqhQ4eUzz77THn22WftHZqBdu3aKREREfrfc3JylKCgIGXmzJkm93nnnXeUTp06mXz+3r17ipOTk/LTTz8ZrG/VqpXy7rvvFhpT586dlaioqMKDV6yL39x90tLSlLp16yq//vprkfNRYYqSrxziSj2vlL8bhfn4+BisX7NmDb6+vjRp0oRJkyZx7949k2VkZmaSmppqsJRaUgUvRKm1f/9+OnfuzNNPP82JEydo2LAh06dPZ/bs2UybNs1g248++ogKFSoUuFy5cqVY4szKyuLIkSN0795dv06tVtO9e3f27dtncr8ffviBNm3a0K9fP/z8/GjZsiUrVqzQP5+dnU1OTg7u7u4G+5UrV47//e9/do3fkn0iIiJ4+umnDbYtbnlzWGZmZqH7OFxS12g0jBkzho4dO9KkSRP9+pdffpmvvvqKnTt3MmnSJL788kteffVVk+XMnDkTb29v/VKtWrWSCL94vPQSqNVw8CCcP2/vaAolQ7kK8Y9x48bRr18/JkyYQN26dRkwYAA///wzzZs3p2XLlgbbvvHGG0RHRxe4FFf30KSkJHJycvD39zdY7+/vz40Cbv1dvHiRpUuXUrduXbZt28bw4cMZNWoUq1evBsDT05MOHTowY8YM4uPjycnJ4auvvmLfvn1cv349X3l5v9j8/vvvvPHGG4V+sbEmfnP3WbduHUePHmXmzJkm34fiUK1aNYM8Zs7xHW5+yoiICE6dOpXvG9ywYcP0j5s2bUpgYCDdunXjwoUL1K5dO185kyZNYty4cfrfU1NTS29i9/eH7t1h+3b4+mt47z17R1SgqKgo/Shy7zl4rEIUp6tXr7Jv3z4+/vhj/TpnZ2cURcl3lQ7a2sm8NZSWmjhxIrNnzy5wmzNnztCgQYMiHUdHo9HQpk0bPvroIwBatmzJqVOnWLZsGYMGDQLgyy+/JCwsjKpVq+Lk5ESrVq0YMGAAR44cyVfeG2+8wUsvvaT//ZVXXqFv37688MIL+nUlOe5FXFwco0eP5tdff81X21ASx8499aqbm1uh+zhUUh8xYgQ//fQTe/bs0c/FbUr79u0BOH/+vNGk7ubmZtYbUGq8/LI2qa9ZA5Mna1vGO6jQ0FDu3r0rQ7mKh96ZM2cAaNWqlX5dTEwM7dq1o2nTpvm2/+ijj/TJ0ZTTp08TEhJi8vm33nqLwYMHF1hGrVq18q3z9fXFycmJhIQEg/UJCQkEBASYLCswMJBGjRoZrGvYsCHffvut/vfatWuze/du0tPTSU1NJTAwkP79+xuNI+8Xm3LlyuHn50edOnUKfE3WxG/OPkeOHCExMdHgb5iTk8OePXtYvHgxmZmZODk5FRibtby8vCyeT90hkrqiKIwcOZJNmzaxa9cuatasWeg+0dHRgPaEeig8/zy8/jrExEB0NOSptnMkMoqcEFopKSk4OTmh+vtL+O3bt/n4449p3ry50e3zXqUaU9j/VZUqVahSpYrFsbq6utK6dWt27NjBc889B2ivwnfs2MGIESNM7texY8d8XZDPnTtH9erV821bvnx5ypcvz507d9i2bZtNezBZE785+3Tr1o2TJ08a7BcaGkqDBg145513ii2hW83mzfasMHz4cMXb21vZtWuXcv36df1y7949RVEU5fz588r06dOVw4cPK7Gxscr333+v1KpVS3n88cfNPkapbv2u8+KL2lbw48fbNQxHbf1++fJlZcyYMcrly5dt8jqFKKq//vpLAZTp06crZ86cUXr06KG0atVK8ff3Vy5dumTv8PJZt26d4ubmpqxatUo5ffq0MmzYMKVixYrKjRs39NssWrRI6dq1q/73gwcPKs7OzsqHH36o/PXXX8qaNWsUDw8P5auvvtJv88svvyhbt25VLl68qGzfvl1p3ry50r59eyUrKytfDGlpaQZ5wNiSnZ1tVfx5Yzf3NeflyK3fHSKpA0YXXTeGK1euKI8//rji4+OjuLm5KXXq1FEmTJhg0QsuE0n9u++0Sb1qVUXJybFbGI6a1KdMmaI8+eSTypQpU2zyOoWwhenTpyuVK1dW3N3dlcGDBytJSUlKq1atlAYNGtg7NKMWLVqkhISEKK6urkq7du2U/fv3Gzw/depUpXr16gbrfvzxR6VJkyaKm5ub0qBBA2X58uUGz69fv16pVauW4urqqgQEBCgRERFKcnKy0eNPnTrVZE7QLbGxsVbFbyx2c15zXo6c1FWKoiglXDlgF6mpqXh7e5OSkmLxPQqHkZEBAQHaceB37YLOne0SRnZ2NlevXtW3e8j7OCAggBs3buR73tnZWb9v7m2cnc2/C5T72Hn3u3LlCvPmzWPs2LEF3nMUQghHVpR85XBd2kQB3N2hb1/tY+mzno/cyxdCPOyKlNQfPHhAXFwcMTEx3L5921YxiYLoBqLZuFE70YsQQgjxN4uTelpaGkuXLqVz5854eXlRo0YNGjZsSJUqVahevTpDhw4ttoH6BdClCwQGwp078PPP9o5GCCGEA7Eoqc+dO5caNWoQFRVF9+7d2bx5M9HR0Zw7d459+/YxdepUsrOz6dGjB08++SR//fVXccX98HJyglde0T7+e8QmIYQQAizsp37o0CH27NlD48aNjT7frl07wsLCWLZsGVFRUfz+++/UrVvXJoGKXAYNgo8/1l6p37wJVvRJFUIIUfZYdKX+9ddfm0zoubm5ufHGG28QFhZmdWCiAE2aQOvWkJ2tHTZWCCFyOXbsGC4uLnTp0qXYj2VqfvGlS5fSrFkz/ahoHTp0YOvWrSbLycnJ4b333qNmzZqUK1eO2rVrM2PGDIqzg5Y5c6PXqFEDlUqVb4mIiCi2uIrCJq3fk5KS+Pnnn/nhhx+MDtAvioFuGEgz5xkWQjw8Ro0axfjx4zl+/HixHqeg+cWDg4OZNWsWR44c4fDhw3Tt2pU+ffrw559/Gi1r9uzZLF26lMWLF3PmzBlmz57NnDlzWLRokVmxWDLvemGx593u+vXr+uXXX38FoF+/fmYfqyQVOal/++231KlTh2nTpjF16lRq165NVFSULWITBRkwAFxc4NgxyDOEoRDi4bV27VoqVapEREQEycnJXLp0qViOc/fuXV555RVWrFhBpUqV8j3/7LPP8tRTT1G3bl3q1avHhx9+SIUKFdi/f7/R8v744w/69OnD008/TY0aNXjxxRfp0aMHBw8eLPHYc6tSpQoBAQH65aeffqJ27dp0ttM4IYWxOKnfvXvX4Pdp06Zx8OBBDh48yLFjx9i4cSPvvvuuzQIUJlSuDM8+q30sDeaEEEB6ejr//ve/mT17NsHBwXh7e+vnycirqHO3WzK/eE5ODuvWrSM9PZ0OHToY3ebRRx9lx44dnDt3DoDjx4/zv//9j169ehX+wi1k7dzoWVlZfPXVV4SFhenH83c0Fk/o0rp1a+bMmUOfPn20BTg7k5iYSL169QDt7Daurq62jVIYN2gQfPcdfPUVzJoFFozMVlZoNBpSU1NL7yiBQtjQRx99xJNPPknDhg0BaNSoEdHR0frJSnIryuQxuvnFC+u+fPLkSTp06EBGRgYVKlRg06ZN+WZ005k4cSKpqak0aNAAJycncnJy+PDDD3lF19snj7wz2t2/f5/9+/cbTN5ibEY7c2M3ZvPmzSQnJxc6C549WZwFtm3bRkREBKtWrWLJkiUsWLCA/v37k5OTQ3Z2Nmq12qL7GqIIevXStnxPSIBt2+Dpp+0dUYnTaDQkJydToUIFe4cihF1dvHiRzz77jFOnTunXNWnSxOSVurVzt1syv3j9+vWJjo4mJSWFb775hkGDBrF7926jiX3Dhg2sWbOGtWvX0rhxY6KjoxkzZgxBQUH6edlzs2be9aLOjR4ZGUmvXr0ce9RKawecX7t2rVKnTh1l4cKFyr1795QTJ04ox44dU+7fv29tkcWqTEzoYsyYMdpJXl58scQO6UgTupiKRYiHTe/evRVAcXJy0i8qlcroBCaKoigffvihUr58+QIXYzMebtq0Kd9xAEWlUilOTk4mZ1BTFEXp1q2bMmzYMKPPBQcHK4sXLzZYN2PGDKV+/fpmvf7OnTvrJwEzpSixX7p0SVGr1crmzZvNiqcoipKvrK6vHTBgAL169WL8+PF06dKF5cuX06JFC1t8zxCWGDQI5s+HH36A27fBim/eQojSbfv27ezdu5djx44ZTHR06NAhwsLCSE5OpmLFigb7WFv9XpT5xTUaDZmZmUafu3fvHmq1YTMvJycnNBpNgTFaoiixR0VF4efnx9MOXiNqVVLfsmULZ86coXnz5nz++efs3r2bV155hV69ejF9+nTKlStn6ziFKS1aQPPmcPw4rF8Pw4fbOyIhRAl68OABY8aMYcKECfkurHRtTaKjo/P1Wbe2+t3T05MmTZoYrCtfvjyVK1c2WD9p0iR69epFSEgIaWlprF27ll27drFt2zYAFi9ezKZNm9ixYwegbS3/4YcfEhISQuPGjTl27Bhz5841Od7J3bt3DRpur1u3DoAbN27o11WpUsUgUZsbe97YNBoNUVFRDBo0yKJZJe3B4tbvb731FqGhoRw6dIjXX3+dGTNm0LlzZ44ePYq7uzstW7YscIABUQx095ukK6EQD53Fixdz69YtgwZiOtWqVcPDw8PkffXilJiYyMCBA6lfvz7dunXj0KFDbNu2jSeeeALQjm9y4cIF/faLFi3ixRdf5M0336Rhw4aMHz9en2OM+fjjjwkMDCxwiYuLsyr2vLH99ttvXLlypVQMqGbxfOqVK1dm+/bttG7dmtu3b/PII4/ouyCAtrXh66+/zu+//27zYIuiTMynbkpiIlStqh1h7sQJaNq0WA/nSPOpm4rF0b9NCyGEKSU6n3r58uWJjY0FtC0J87YgbNSokcMl9DLPzw/+7mLI55/bNxYhhBB2Y3FSnzlzJgMHDiQoKIjOnTubrBoRJWzIEO3PL7+EjAz7xiKEEMIuLE7qr7zyCnFxcXz//fdcunRJPwiNsLMnnoCQEO086999Z+9ohBBC2IFFSV03ZGDlypVp27Ztvi4SeV27ds3qwISFnJxA14hjxQr7xiKEEMIuLErqbdu25fXXXy9weL2UlBRWrFhBkyZN+Pbbb4scoLBAWBioVLBrF/z1l72jEUIIUcIsaiJ8+vRpPvzwQ5544gnc3d1p3bo1QUFBuLu7c+fOHU6fPs2ff/5Jq1atmDNnDk899VRxxS2MqVYNnnwStm6FyEjtePBCCCEeGhZdqVeuXJm5c+dy/fp1Fi9eTN26dUlKSuKvv68KX3nlFY4cOcK+ffskodvL0KHan6tWwYMHdg1FCCFEybKqM2+5cuV48cUXefHFF20djyiqZ54Bf3/tJC8//QTPP2/viIQQgri4OF577TUSExNxdnbmvffeo1+/fvYOq8yxuPW7cHAuLqCbFrAMNphTq9VUrFgx3xjRQgjH5uzszPz58zl9+jTbt29nzJgxpKen2zusMkc+Gcui8HDtz19+gb97LJQVktSFKJ0CAwP1Y9MHBATg6+vL7du37RtUGSSfjGVR3brwf/8HilIsV+uSWIUQRXHkyBFycnKoVq2azcu+du0ar776KpUrV6ZcuXI0bdqUw4cPm9z+/fffR6VSGSwNGjTIt92SJUuoUaMG7u7utG/fnoMHD9o8dluQT+WySjdb24oVkJVl06IlqQshrHX79m0GDhzI8uXLbV72nTt36NixIy4uLmzdupXTp0/zySefUKlSpQL3a9y4MdevX9cv//vf/wyeX79+PePGjWPq1KkcPXqU5s2b07NnTxITE23+GorKqk/ly5cv8+OPPzrsNxUBPPccBARoG8xt3mzvaIQQZVy/fv2oUqWKQbI+cOAArq6ubN++HYDMzEyee+45Jk6cyKOPPmrzGGbPnk21atWIioqiXbt21KxZkx49elC7du0C93N2diYgIEC/+Pr6Gjw/d+5chg4dSmhoKI0aNWLZsmV4eHiwcuVKm7+GorI4qX/99dfUq1ePPn360KFDB9q0acPNmzeLIzZRFC4u/3Rv+/RT+8YihCjzFi5cSN++fZk+fTqgne/81VdfZfjw4fTo0QNFURg8eDBdu3bltddeM1nORx99RIUKFQpcrphoK/TDDz/Qpk0b+vXrh5+fHy1btmSFGbcg//rrL4KCgqhVqxavvPKKQflZWVkcOXKE7t2769ep1Wq6d+/Ovn37zH17SozFSX3atGm8/PLLnD17Vv/ta+LEiTYPTNjAsGHa4WN374Y//7R3NEKIMiwwMJAxY8Zw7do1bt26xahRo3Bzc2P27NkA7N27l/Xr17N582ZatGhBixYtOHnyZL5y3njjDaKjowtcgoKCjMZw8eJFli5dSt26ddm2bRvDhw9n1KhRrF692mTc7du3Z9WqVfzyyy8sXbqU2NhYHnvsMdLS0gDt3Oo5OTn4+/sb7Ofv78+NGzesfbuKjcX91C9evMgvv/xCjRo1qFevHl999RWtW7cmMjKyOOITRREcDL17w6ZNsGwZLFpk74iEEGVYvXr18PDwYMqUKaxZs4aDBw/qp+fu1KkTGo2m0DJ8fHzw8fGx6vgajYY2bdrw0UcfAdCyZUtOnTrFsmXLGDRokNF9evXqpX/crFkz2rdvT/Xq1dmwYQPhup5EpYjFV+rZ2dl4eHjof2/QoAEajcYhv7EI/mkwt3o13L1r31iEEGWaWq2madOmfPrpp3zwwQc0b97c4jKKUv0eGBhIo0aNDNY1bNjQ5PbGVKxYkXr16nH+/HkAfH19cXJyIiEhwWC7hIQEAgICLHx1xc+qhnKrV6/mjz/+4O7fScLZ2Zl79+7ZNDBhI926abu4paXB2rX2jkYIUYYpigJAq1ateOutt0xul5OTY/K5olS/d+zYkZiYGIN1586do3r16ma/hrt373LhwgUCAwMBcHV1pXXr1uzYsUO/jUajYceOHXTo0MHsckuKxUn9scce44MPPqBTp05UrFiRunXrkpGRQWRkJDt37tTfhxAOQq3+52p9yRJt33UhhCgG8+fP58CBA2g0mnxdXnv37s2bb75J27ZtiYqKMlmGj48PderUKXBxdjZ+53js2LHs37+fjz76iPPnz7N27VqWL19OREQEAIsXL6Zbt24G+4wfP57du3dz6dIl/vjjD55//nmcnJwYMGCAfptx48axYsUKVq9ezZkzZxg+fDjp6emEhoZa+1YVH8VK586dU77++mtlwoQJSrdu3ZRKlSopKpVKcXJyUho0aGBtscUmJSVFAZSUlBR7h1Lybt1SFHd3RQFF+d//bFr0gwcPlNjYWOXBgwf6x/fv38+37sGDBwbb597G1scvaplCCMudOHFCcXNzU958803F1dU13/9hjRo1lPnz5xd7HD/++KPSpEkTxc3NTWnQoIGyfPly/XNTp05VqlevbrB9//79lcDAQMXV1VWpWrWq0r9/f+X8+fP5yl20aJESEhKiuLq6Ku3atVP2799fbK+hKPlKpSi2u3SLjY3l8OHDHDt2TN9QwVGkpqbi7e1NSkoKXl5e9g6n5IWHw8qV0K8fbNhgs2Kzs7O5evUqwcHBAFy9epWAgABu3LhBcHAwarWa1NRUvLy8UKvV+u1zb2PqW7e1xy9qmUIIy2RkZNC2bVtatWrF/Pnz8fHx4eTJkzRp0gSAtLQ0GjdubNG97YdZUfKVTYcEq1mzJv369XO4hC6A0aO1P7/9Fi5fLrHD5h19TkajE6LsmThxIunp6SxevJhKlSpRvXp15s+fT3x8PAB//vlnsQw2I/JziE/WmTNn0rZtWzw9PfHz8+O5557L19ghIyODiIgIKleuTIUKFejbt2++1oiiAM2aQdeuoNFo763biaVJXaPRkJycbFZXGCFEydu+fTtLlizhq6++wtPTE4DJkyezefNm/b3skydP0rRpU3uG+dBwiKS+e/duIiIi2L9/P7/++isPHjygR48eBtPyjR07lh9//JGNGzeye/du4uPjeeGFF+wYdSk0Zoz254oVUEqmPJSkLoRj69GjBw8ePDC4Eh8yZAhJSUls2rQJkKRekhzixuMvv/xi8PuqVavw8/PjyJEjPP7446SkpBAZGcnatWvp2rUrAFFRUTRs2JD9+/fzyCOP2CPs0ufpp6F2bbhwAb744p9W8XYQHx/PggULGDt2LCEhIWbto9FoDO7PCyFKh4ULF9o7hIeGQ34ypqSkAOhHFTpy5AgPHjwwGHu3QYMGhISEmBx7NzMzk9TUVIPloadWw6hR2scLFmir4u0kMjKSs2fPWjQSobGrdt2XA929OyGEKCvy5rDMzMxC9ylSUt+xYwf//ve/GTJkCGFhYQaLtTQaDWPGjKFjx476lpM3btzA1dWVihUrGmxb0Ni7M2fOxNvbW78Ux7y9pVJoKHh5QUwMbNtm06J1CfbGjRuF3jcPDw+nQYMGZg3DWFAVvDVfDoQQojSoVq2aQR6bOXNmoftYndSnTZtGjx492LFjB0lJSdy5c8dgsVZERASnTp1i3bp1VpcBMGnSJFJSUvRLXFxckcorMzw9td3bAObPt2nRugQbFRVVaFIPCgpi9OjRJkeGyk1RFFJSUowmdUu+HAghRGkSFxdnkMcmTZpU6D5W31NftmwZq1atKnAKPUuNGDGCn376iT179uj7HQMEBASQlZVFcnKywdV6QWPvurm54ebmZrPYypSRI7XV79u3w+nTkGesZGuFh4eTmppaognWki8HQghRmnh5eZVcP/WsrCyb9TtUFIURI0awadMm/vvf/1KzZk2D51u3bo2Li4vB2LsxMTFcuXLFIcfedXg1a2pnbwOwYQMWSbBCCGFfVif1IUOGsNZGE4RERETw1VdfsXbtWjw9Pblx4wY3btzg/v37AHh7exMeHs64cePYuXMnR44cITQ0lA4dOkjLd2vpurd98QXcumXXUIQQQtiG1dXvGRkZLF++nN9++41mzZrh4uJi8PzcuXPNLmvp0qUAdOnSxWB9VFQUgwcPBmDevHmo1Wr69u1LZmYmPXv25NNPP7U2fPH449CiBURHa/utT5xo74iEEEIUkdVJ/cSJE7Ro0QKAU6dOGTynUqksKsuc4efd3d1ZsmQJS+w4GlqZolJpr9YHD4ZFi2DsWJA2CEIIUapZndR37txpyziEPfzrXzBpEsTHa+dad8RpBIUQQpjNIQefESXEze2fe+tz5th1MBohhBBFV6RhYpOTk4mMjOTMmTMANGrUiPDwcLy9vW0SnCgBr78OH34IZ8/Cjz9Cnz72jkgIIYSVrL5SP3z4MLVr12bevHncvn2b27dvM2/ePGrXrs3Ro0dtGaMoTt7e/4wBP3s2mNG+ITeZSlUIIRyH1Z/EY8eOpXfv3ly6dInvvvuO7777jtjYWJ555hnG6Kp0RekwejS4usK+ffC//1m0qyMldUeKRQgh7KFIV+rvvPMOzs7/1OA7Ozvz9ttvc/jwYZsEJ0pIYCAMGqR9PHu2fWMpAknqQoiHndWffl5eXly5ciXf+ri4ODw9PYsUlLCDCRO03dx+/hnydFEUQghROlid1Pv37094eDjr168nLi6OuLg41q1bx5AhQxgwYIAtYxQloW5d6NtX+3jOHPvGIoQQwipWt37/+OOPUalUDBw4kOzsbABcXFwYPnw4s2bNslmAogS98w588w18/TV88AGEhNg7IiGEEBaw+krd1dWVBQsWcOfOHaKjo4mOjta3gJfZ0UqpNm2ga1fIzgYLhvkVQgjhGIrcosjDw4OmTZvStGlTPDw8bBGTsKd33tH+XLECbt60byxCCCEsYlH1+7hx45gxYwbly5dn3LhxBW5ryYQuwoE88YT2iv3wYe3V+syZ9o5ICCGEmSxK6seOHePBgwf6x6ZYOqGLcCAqFUyZop1vffFiGD8eKle2d1RCCCHMYFFSzz2Jy+rVqwkODs7XJ1hRFOLi4mwTnbCPZ575Z1rW+fNhxgy7haLRaLh7967c2hFCCDNYfU+9Zs2aJCUl5Vt/+/ZtatasWaSghJ3prtYBFi6EO3fsFoouqZszPa8QQjzsrE7qpj5k7969i7u7u9UBCQfRpw80bQqpqdrELoQQwuFZ3E9d10BOpVIxZcoUg2rRnJwcDhw4QIsWLWwWoLATtRreew9eeklbBT9mjHbyFyGEEA7L4qSuayCnKAonT57E1dVV/5yrqyvNmzdn/PjxtotQ2E/fvtCoEZw+DYsWweTJ9o5ICCFEASxO6rrGcqGhoSxYsAAvLy+bByUchFqtTeQvv6zt3jZ6NMi4/kII4bCsvqceFRUlCf1h8NJLUK+etrHc4sX2jkYIIUQBrE7qM2fOZOXKlfnWr1y5ktmlePpOkYeT0z/V7v/5D6Sk2DceO9FoNCQnJ6PRaOwdihBCmGR1Uv/ss89o0KBBvvWNGzdm2bJlRQpKOJiXX4aGDbVX6w/pSIGS1IUQpYHVSf3GjRsEBgbmW1+lShWuX79epKCEg3FygunTtY/nzgUj4xM4Ckm+QoiHmdVJvVq1auzduzff+r179xIUFFSkoIQDeuEFaNkS7t4FB769IkldCPEwszqpDx06lDFjxhAVFcXly5e5fPkyK1euZOzYsQwdOtSWMQpHoFZr51gHbYO5+Hj7xiOEECIfi7u06UyYMIFbt27x5ptvkpWVhaIolCtXjnfeeYeJEyfaMkbhKHr1gg4dYN8++OgjaQ0vhBAOxuordZVKxezZs7l58yb79+/nxIkT3L59mylTpsgsbWWVSgUffqh9vHw5XLpk13AA1Go1FStWzDexkBBCPIysvlLXuXLlCrdu3SIrK4uLFy/q1/fu3buoRQtH9H//B926wY4d2sZzRro1liRdUhdCCFGEpH7x4kWef/55Tp48iUql0k/wortKz8nJsU2EwvF8+KE2qa9eDW+9BY0bF+vhkpKSWL9+PWPGjJErciGEKIDVn5CjR4+mZs2aJCYm4uHhwalTp9izZw9t2rRh165dNgxROJz27eH550GjgXfesWnRxlqv//DDD8TExBAVFWXTYwkhRFljdVLft28f06dPx9fXF7VajZOTE506dWLmzJmMGjXKljEKRzRrFjg7w88/w9/zAdiCsaTeu3dv6tevT2hoqM2OI4QQZZHVST0nJwfPvyf38PX1Jf7vLk7Vq1cnJibGNtEJx1WvHrz+uvbxhAnaq/Zi4uvry4gRI2T8AyGEKITVSb1JkyYcP34cgPbt2zNnzhz27t3L9OnTqVWrls0CFA5syhTtrG1HjsC6dfaORgghHnpWJ/XJkyfrq0inT59ObGwsjz32GFu2bGHhwoU2C1A4MD+/f+6p//vfkJlZrIeLj49nwYIF+lohIYQQhqxO6j179uSFF14AoE6dOpw9e5akpCQSExPp2rWrzQIUDm7sWAgKgsuXUX36qc2Lj4+PZ+XKlSQkJBAVFcXZs2eJjIy0+XF0ZJhZIURpZlVSf/DgAd26deOvv/4yWO/j4yMDzzxsPDxgxgwA1B99hDo52abFR0VFceHCBTZu3EhoaCgNGjQgPDy8yOWaSt6S1IUQpZlVSd3FxYUTJ07YOhZRWg0aBE2bokpOxnvJEpsWHRoaSu3atenXrx9BQUGMHj3aJg3mJHkLIcoiq6vfX3311WKtBhWliJMTzJkDgNcXX0BsrM2KDgoKIiwsDH9/f5uVKYQQZZXVST07O5ulS5fSpk0bXn/9dcaNG2ewWGLPnj08++yzBAUFoVKp2Lx5s8HzgwcPRqVSGSxPPvmktaGL4tCzJ5pu3VBlZaG28YA0QgghzGP1MLGnTp2iVatWAJw7d87gOUvvq6enp9O8eXPCwsL0je/yevLJJw1GFHNzc7MwYlGsVCo0H3+Mqk0b1Js2aYeR7dbN3lEJIcRDxeqkvtOGo4j16tWLXr16FbiNm5sbAQEBNjumsC2NRkNycDAur7yC9xdfwOjREB2tHXVOCCFEibC4+v3ixYv6yVtK0q5du/Dz86N+/foMHz6cW7duFbh9ZmYmqampBosoPhqNhpSUFG6PHo1SuTL8+ScsXWrvsIQQotTKm8MyzRgLxOKkXrduXW7evKn/vX///iQkJFhajEWefPJJvvjiC3bs2MHs2bPZvXs3vXr1KnAmuJkzZ+Lt7a1fqlWrVqwxCi2NtzeaadO0v0yZAklJ9g1ICCFKqWrVqhnksZkzZxa6j8VJPe9V+pYtW0hPT7e0GIv861//onfv3jRt2pTnnnuOn376iUOHDhU4G9ykSZNISUnRL3FxccUao/iHMmQING8Oycnw3nv2DscmZDQ7IURJi4uLM8hjkyZNKnSfUjk5da1atfD19eX8+fMmt3Fzc8PLy8tgESXEyQkWLdI+/uwz7b31PNRqNRUrVrR6fvSS7mceGRlZ7KPZCSFEbnlzmDkNxC3+RNV1Kcu7riRdvXqVW7duERgYWKLHFRZ47DH4179AUWDUKO3PXEpbUg8PD7fZaHZCCFFcLG6arCgKgwcP1n9jyMjI4I033qB8+fIG23333Xdml3n37l2Dq+7Y2Fiio6Px8fHBx8eHadOm0bdvXwICArhw4QJvv/02derUoWfPnpaGL0rSnDnw/ffw+++wZg28+qq9I7KaLUezE0KI4mLxZdKgQYPw8/PT37h/9dVXCQoKMriZ7+3tbVGZhw8fpmXLlrRs2RKAcePG0bJlS6ZMmYKTkxMnTpygd+/e1KtXj/DwcFq3bs3vv/8ufdUdXbVq/9xTf+stuHPHvvGYSYaQFUKUVhZfqeceAMZWunTpUmA3uW3bttn8mKKEvPUWfPklnDmjnZ61CN3cNBoNd+/eRaPRWF1tbw5FUUhJSSnS7QEhhLAH+cQSxcvVFZYt0z7+7DPYv9+i3XNPvZo7qQshhMhPkroofo8/rp3JTVHgjTcgO9voZsaqvXNPvSqEEKJgktRFyfjPf6BSJTh+HBYvNrqJblS63Ldick+9KoQQomCS1EXJqFJFPz0r770HV6+atZtMvSqEEOaTpC5KTlgYPPoo3L0LI0fm67suhBCiaIo0hdaOHTvYsWMHiYmJ+RovrVy5skiBiTJIrdY2mmvVCjZvho0b4aWX7B2VEEKUGVZfqU+bNo0ePXqwY8cOkpKSuHPnjsEihFFNm2q7tgGMGCETvgghhA1ZfaW+bNkyVq1axWuvvWbLeMTD4N134dtvtdOzjhkDX31l74iEEKJMsPpKPSsri0cffdSWsYhSTqVS4e3tXfiALa6usHKltjp+zRr4+eeSCVAIIco4q5P6kCFDWLt2rS1jEaWcRZO0tGsHY8dqH7/+OqSkFG9wQgjxELC6+j0jI4Ply5fz22+/0axZM1xcXAyenzt3bpGDE2Xc9OnaCV/On4e334YlS4pcpG7e87FjxxISEmKDIIUQovSw+kr9xIkTtGjRArVazalTpzh27Jh+iTYyf7YQ+Xh4wOefax8vX45q164iFynzngshHmZWX6nv3LnTlnGIh1XnztqhY5ctQz1sGKpNm1AqVLC6uPDwcFJTU2XecyHEQ0kGnxH2N3s2hISgio3lwXvvsXDhQuLj460qSuY9F0I8zIo0+ExycjKRkZGcOXMGgEaNGhEeHm7xfOriIeflBatXo3Ttyroff+R8y5ZEVqzItGnTbHoYjUZDamoqXl5eNi1XCCEchdVX6ocPH6Z27drMmzeP27dvc/v2bebNm0ft2rU5evSoLWMUD4MuXVBGj2Yo0ObsWcKfe07/lEqlwsvLy6q5zXPP/GZsFjghhChLrE7qY8eOpXfv3ly6dInvvvuO7777jtjYWJ555hnGjBljwxDFw0IzYwa+desy7f59gqdP148Nb6yrnK6Ve2HV9JYm8vj4eBYuXEhCQoL1L0QIIeykSFfq77zzDs7O/9TgOzs78/bbb3P48GGbBCceMu7uJHz8MYqLC+rNm+GLL0xuWtRW7qa+FERGRhITE8P69eutKlcIIezJ6qTu5eXFlStX8q2Pi4vD09OzSEGJh1dW48bcHjVK+8vIkXD5stHtwsPDadCggdWt3E19KQgPD6d+/fr079/fqnKFEMKerE7q/fv3Jzw8nPXr1xMXF0dcXBzr1q1jyJAhDBgwwJYxiodM8rBhKB06QFoarmFhkJ2db5uitnI39aUgKCiIUaNGyfztQohSyerW7x9//DEqlYqBAweS/feHrouLC8OHD2fWrFk2C1A8hJycyFm1Cuc2bVD/8QfBUVHQpo1ND2HplwKLhsAVQgg7sfoTytXVlQULFnDnzh2io6OJjo7Wt4B3c3OzZYziYVSrFnz2GQBVo6JQ795t13AkqQshSoMif0J5eHjQtGlTmjZtioeHhy1iEkJrwACyBw1CpSi4hIbK3OtCCFEIi6rfx40bx4wZMyhfvjzjxo0rcFuZ0EXYQvYnn/Bg1y7KXb6My7Bh8N579g5JCCEclkVJ/dixYzx48ED/2BSVSlW0qITQKV+ev2bMoOnQoTht3UpAvXrQsqXFxehGkzOnv3pCQgLr16/nrbfekpnehBClikVJPfckLqtXryY4ODjfPUZFUYiLi7NNdEIA9+rWJXv2bFzGjCFkyRIevPSSdj52E3QDzvj4+BisS0lJ0T9X0FCx69evJyYmhsjISJsPVSuEEMXJ6nvqNWvWJMnIPc7bt29Ts2bNIgUlSh9zR3izVs6wYeT07o06OxuX116D5GST2yqKok/gljyn079/f+rXry8zvQkhSh2rk7ry9xCeed29exd3d3erAxKlU7HPY65ScXnqVN6rUIEbFy+iDg2FYhrD3d/fn1GjRslMb0KIUsfifuq6BnIqlYopU6YYtHjPycnhwIEDtGjRwmYBitKhJOYxX7V5MwcbNcLpyBHe/+knvBs0gJkzi+14QghR2lic1HUN5BRF4eTJk7i6uuqfc3V1pXnz5owfP952EYpSoSTmMR8yZAhpaWmEv/ACTJxIxblz0XTvDj17FtsxhRCiNLE4qesay4WGhrJgwQKZm1qUGN0Xh8DgYDRnz6JetQr1q6/C0aNQrZq9wxNCCLuzepjYqKgoAE6fPs2VK1fIysoyeL53795Fi0yIAmgWLuTB4cO4nToFL74Ie/aAjGQohHjIWZ3UY2Njee655zh58iQqlUrfcE7XRz0nJ8c2EQphTLly3Pz0U6r26YPq4EEYNQqWLbNrSLq+8F5eXjKcrBDCLqz+5Bk1ahQ1a9YkMTERDw8P/vzzT/bs2UObNm3YtWuXDUMUDzO1Wk2FChWMJsnsatXQfPEFqFSwfDl8+qnNj6/r127OoDWWbCuEEMXB6qS+b98+pk+fjq+vL2q1GrVaTadOnZg5cyajdPNhC2GB+Ph4Fi5cSEJCgn5dQUkdQHnySdDNCjh6NKr//temMUmiFkKUJlYn9ZycHDw9PQHw9fXVDzpSvXp1YmJibBOdeKhERkYSExPD+vXrLdtxwgR47TXIyUH9r3/hfOlSscQnhBCOzuqk3qRJE44fPw5A+/btmTNnDnv37mX69OnUqlXLZgGKh0d4eDj169enf//+lu2oq35v3x7VnTsEvv466rS0fJslJCSwePFig5oAIYQoS6xO6pMnT9ZXSU6fPp3Y2Fgee+wxtmzZwsKFC20WoHh4BAUFMWrUKPz9/S3f2d0dNm1CqVoV1wsX8BszBv5urKkb933dunXW1QQIIUQpYXVS79mzJy+88AIAderU4ezZsyQlJZGYmEjXrl0tKmvPnj08++yzBAUFoVKp2Lx5s8HziqIwZcoUAgMDKVeuHN27d+evv/6yNnRRVgUGkvPtt2jc3Sm/ezfqiROBf8Z779evn3U1AUIIUUrYtN+Nj4+PVdOupqen07x5c5YsWWL0+Tlz5rBw4UKWLVvGgQMHKF++PD179iQjI6OoIQsHp1arqVixYr6GcqbW07o1ibNna7eZP5/rM2fqG9/5+/szYsQI62oChBCiFLA6qXft2tXotJR37tyx+Eq9V69efPDBBzz//PP5nlMUhfnz5zN58mT69OlDs2bN+OKLL4iPj893RS/KHouTOpD+9NPc+nuOgsipUzn/++9ERUXZ5X66tJ4XQpQkq5P6rl27WLx4Mc899xzp6en69VlZWezevdsmwYF2kJsbN27QvXt3/Tpvb2/at2/Pvn37TO6XmZlJamqqwSIcn0aj4e7du0VOgsnDh6MJC2OYotD25Emck5KIiYkhKiqKhQsXFtsUsXlJUhdCWCtvDsvMzCx0nyJVv//222/cuHGDRx55hEvF1I3oxo0bAPmqTP39/fXPGTNz5ky8vb31SzUZG7xUUBSlyEk9ISGBhYsWcfXddwno2ZP3s7OZ/PvvNAwIANAnd0m4QghHVq1aNYM8NtOMWSmLlNQDAwPZvXs3TZs2pW3btg41ktykSZNISUnRL3FxcfYOSZSQ9evXExMTQ+Tq1Wi+/prMRo0ISU5mzp49hD//PPXr1yc0NFTfgK6oSV2+HAghikNcXJxBHps0aVKh+1id1HUN4tzc3Fi7di2jR4/mySef5FMbD9UZ8PfVVd57oQkJCfrnjHFzc8PLy8tgEQ+H/v37U79+fe3c7p6eXP/8cx4EBeF66RKtpk5l9NChNp0itqhJXb4UCCGMyZvD3MyYtMrqpK6bwEVn8uTJrFmzhk8++cTaIo2qWbMmAQEB7NixQ78uNTWVAwcO0KFDB5seS5QN/v7+jBo1Sp+4c/z8uLZiBTleXpQ7dgz/UaPgwQM7R/kPSepCCFuxOqnHxsbi6+trsK5v377s37+flStXWlTW3bt3iY6OJjo6Wl92dHQ0V65cQaVSMWbMGD744AN++OEHTp48ycCBAwkKCuK5556zNnzxkMmqU4drS5dq+7Dv3InLkCH6wWnsRZK5EMLWrJ56tXr16kbXN2nShCZNmlhU1uHDh/m///s//e/j/u6ONGjQIFatWsXbb79Neno6w4YNIzk5mU6dOvHLL7/g7u5ubfjiIZTRpg3XFy4k6M03cdqwAV9nZ5KMdMssKbqkLuexEMJWLErq48aNY8aMGZQvX16feE2ZO3eu2eV26dIlX3V+biqViunTpzN9+nSzyxSlk0qlKnBWtqJK79yZxI8/xm/sWLzXrkXj7Q25hjVWq9V4e3tLF0ghRKlkUVI/duwYD/6+F3ns2DGT21kzqpwQUPhUq7Zw95ln8HFywmXkSCotXUpOjRrw9tv641esWJG7d+8W2/GFEKK4WJTUd+7cCcCDBw9Qq9UsW7aMunXrFktgQljD3CttZdgw0hMSKP/BBzi98w6ULw8RERbd59ZoNKSmpuLh4WGr8IUQokisuhxycXHhxIkTto5FiCIraPjYvNu5TZ1K8vDh2hUjRsCyZWRnZxMfH29WbYE0dBNCOBqr6zhfffVVIiMjbRmLECUuecIENLr2IcOHo1qxgnv37uHl5VWstwCEEKI4WN36PTs7m5UrV/Lbb7/RunVrypcvb/C8JQ3lhLAblQrNrFmoFQXmzcNlxAiqTJoELVqUeCi6K/+CvlDoqvzlS4cQwhirk/qpU6do1aoVAOfOnTN4ThrKiVJFpYJPPgGNBhYsoNasWWTXrAlDh5ZoGLqkXlDVvznbCCEeXlYndV2jOSHKBJUK5s0jOzMT52XLcH79de26wYPNLsJWM8wJIYS15Ku+KDXMbQSXV0JCgnlzqatUZM+dy42+fVEpCgwdimrBArOPc+3aNRYsWMDVq1ctik8IIWylyEn99OnT/PLLL/zwww8GixC2Zm1S183a9sMPPxjcGoqPj2fBggWGc6urVFx66y2yx44FIGH8eJa9/DLx164VepxVq1Zx/vx5oqKiLIqvJOmq77Ozs6XlvhBlkNXV7xcvXuT555/n5MmTqFQq/Yhwug/NHDuPqy2ETv/+/Vm/fj3Dhw83GLkwMjKSs2fPEhkZybTcw8WqVGR/+CHOlSqxYsoU/tq3j5WvvML7v/+urZI3YfDgwSQmJhIaGlqcL6dIcg9NK/fmhSh7rP5vHj16NDVr1iQxMREPDw/+/PNP9uzZQ5s2bRxqXnUh8s7aphMeHk6DBg20U7TmpVLBe+8R9u67NAOG7d2LMnIkybdvm7y6DQwMJDw83OJpXY3WGAghhBWsTur79u1j+vTp+Pr6olarUavVdOrUiZkzZzJq1ChbxigeUrYe3EWtVlOpUiUqVaqEWq0mKCiI0aNHF5iEA99/n4gPPyRQpUK1ZAnOYWFoMjIAw2RsTay6/RcuXMjZs2cNqu1lYBshhDWsTuo5OTl4enoC4Ovrq7/KqF69OjExMbaJTjzUrE1spu69q9VqfHx88PHxsajK+e7LL6NZtQrF2ZkK33+PundvSEvTV99HRUWhKAppaWkmJyYy9lp0+wM0aNDAoNpekroQwhpW31Nv0qQJx48fp2bNmrRv3545c+bg6urK8uXLqVWrli1jFKWAtY3YijMWwGZJUXn5ZTS+vqj69UO9Ywd06UL4Z5+RmprKU089xcKFC3nhhRdwdjb+L2Wsf3l4eDipqamMHDkStVpNQEAAN27csEm8QoiHk9WfwJMnT9Z/YE6fPp3Y2Fgee+wxtmzZwsJcU1mKh4MjJfXiovTowY2vv0apUgWOHqXaSy/xVu/ebNmyhZiYGDZu3GiwfWFX28aq/21xhS5X+UI8vKz+BO7ZsycvvPACAHXq1OHs2bMkJSWRmJhI165dbRagEEWhm7XNVqMcZjVrRs6ePVC7NqrYWAJffJGwli2pX78+/fr1M9hWo9Fw7tw5ixrBKYpCSkpKkRJydnY2V69eJTs72+oyhBClk00vq3x8fGSIWFEkCQkJLFy40GYtwS2tQdBoNNy/f7/gpFqnDuzdi9KqFU63b1PjlVeYVLs2fn5++TbV9ZEvyb7rMrKdEA8vi+6pj9PNZmUGmdBFWEOXBCMjI3nvvfdK/PjXrl0jMjKS4OBg6tata3pDf39yduwgq29fPP77XwJGj8YpIoLb//qXwWa6PvKO3HddCFF2WHSlfuzYMbOW6OjoYgpXlHX9+/enfv36xvuOl4BVq1YRGxvLqlWrCt+4QgUSly8ne/RoAKosWUK9adPg/n39Jqb6yFurqH3adfsfO3ZM+sYLUQZZdKUuk7iI4mbrJGipXr168eOPP9KrVy/zdnByInvWLG77+VFl6lSq/PYbmp494YcfwNfX5vGZHAUvD11jOT8/P4NbD7r9J06ciIeHR6HlCCFKF6u7tAEkJycTGRnJmTNnAGjcuDFhYWF4e3vbJDjxcNE1aktNTbVbDFu3bsXT05OtW7fSqVMns/dL69+f+0FBBI0ahcuhQ9CuHWzcWOTEnneOdV03uMJqMhRFITU1VT84lI5u/5deeokNGzbYrUZECFE8rG4od/jwYWrXrs28efO4ffs2t2/fZu7cudSuXZujR4/aMkbxkHCEbnGDBw+mZs2aDLZgylWde+3bc2L5cjR160JcHE6dO+P57bdm7x8fH8/ChQsNZpPL2z3NnFHwQJvUjTWW0+3fsmVLs8oRQpQuVn96jh07lt69e3Pp0iW+++47vvvuO2JjY3nmmWcYM2aMDUMUouQEBgby8ssvExgYaNX+GdWqkbVnDzzzDKrMTPzeeQffKVMgK6vQfaOiooiJiWH9+vUG663pd24qqRdE+rcLUfoV6Ur9nXfeMRhBy9nZmbfffpvDhw/bJDghHFnumoWEhASWLFlCUlISVKwI339PztSpKCoV3mvX4vrEEzgVMlpcaGgo9evXp3///gbrzem7XlhCNidhS1IXovSzOql7eXlx5cqVfOvj4uL0Y8ILUZbkbXmeO6mvX7+eEydOMG3aNO0EL8CdESOI/+wzcry8UB88SNAzz8CePSbLDwoKYtSoUVSpUsWs5Jo7CdsiqVuynRDCMVmd1Pv37094eDjr168nLi6OuLg41q1bx5AhQxgwYIAtYxSiWFh6Dz93y/O8+vfvT0ZGBm5ubkRFRaHRaEhJSSG9c2eubtqEpmlTnG7dwumJJ2D6dMjJMXkcXSM3U5PD6BRHArbFiHZCCPuxOql//PHHvPDCCwwcOJAaNWpQo0YNBg8ezIsvvsjs2bNtGaMQxcLSpB4eHq4fDjZv0vP39+c///kPjRo1yjfQTHb16lz6+mvebdiQ6xoNTJ0K3btDAX3EizKyXnx8PFFRUdpbAcjVtxAPE6uTuqurKwsWLODOnTtER0cTHR3N7du3mTdvHm5ubraMUQiHEBQUxMiRIylXrpzRBOnn50dYWJjRFuVR69ZxtHp1PuvTB8qXh127cGrdGt9Dh4x+qdi4caPVw8tGRUVx4cIFfvrpJwB9rUFhV/7FSb5YCFEyLE7qGo2G2bNn07FjR9q2bcv06dOpU6cOTZs2xcPDozhiFKJILLkiv379OmvXruX69es2TUShoaE0aNCA8IUL4ehRaNECVVISFV56CZd//5uKHh4G8fXr14/69etbNbxsaGgotWvX5plnngGKNha8rd4DSepClAyLk/qHH37Iv//9bypUqEDVqlVZsGABERERxRGbeMgVdUhUHUuSeu5hYm15hWvQv7xePdi3D0aMAED1ySdU7NUL9ZkzJCQksHjxYgCrR9YLCgoiNDQU378HvtEl9ZycHJKTk8nOziY1NVU/oE1BJBkLUbpYnNS/+OILPv30U7Zt28bmzZv58ccfWbNmjfzTC5srqGFacSnK4DMWcXeHRYtg0ybw8YGjR3Fq357v33mHmLNn883NbqmkpCRWrVpl8IVI1wBPo9GQlpZmNKnbeoY3+VIgRMmyOKlfuXKFp556Sv979+7dUalUMjGEsLnw8HBtlXUJDmVa1MFnLPbcc3DqFDz1FKrMTEb/73+0PX2aARYMUWvMjz/+yIULFwzuyWs0GtLT06lQoYLJKZKtGbSmIJLUhShZFif17Oxs3N3dDda5uLjw4MEDmwUlBJg/JKqj0lWl5x721ajAQPjpJ3KWLiXAw4MPrl2jXVgYnt98A1ZU/cfHx3Pv3j0CAwMN7skrisK9e/dMVrvHx8ezePFifat5IUTpY/GELoqiMHjwYIMW7hkZGbzxxhuUL19ev+67776zTYRC2JFKpdInwcKuNlUqFRUqVECtVhMfH8+4ceNwc3NDpVLRqlWrwg6EMmQIV+vVw2/CBModPYrfxInk/Pe/OL/7LtlVqwKG7QNMxRMVFcXly5fJyMiw6LVGRkYSExNDeno6PXv2tGjf3DQajf6evRCiZFl8pT5o0CD8/Pzw9vbWL6+++ipBQUEG64QojdRqNeXKldNfyVrSyE6tVuuTemRkJO7u7mRkZOQb9rWgKuns6tWJ++orbo4fj8bVlYTt24nq3p37ixdDTo5Z8YSGhhoMhGOu8PBwqlatyr1794p0O83Y68vd6NGSKnmpvhfCMhZfqVvTb1aI0iJvUjd3H29vb+7cuaNfFx4eTnJyMi+++CL+/v4G2+sSle4LQF4JSUkszszktdWr+fnttzkZF4dq/nym7NsHkZHQuHGB8QQFBTFlyhR+/PHHfF3iCkqSQUFBVKhQgYsXLxIVFcWMGTPMfg8KExUVpW/0+N577xX4+nPHevv2bVJTUwvdVgihJf8lQhSRsavngIAAwsLCqFKlisXl6Qae+WrfPgbu3Ent//s/wjw8UB04AC1bakeky8gwaKmet/ufr68vgwcPztceQVEU0tLSTHbT69evH7Vr17aqf7wxurieeuopixs9mtulUK7mhfiHJHUhioEu4VrTx1038Ez//v0JqlqV11esIGf7djTPPgsPHmjHjm/cGH7+WZ/UC+r+Fx8fz8qVK7lx4wYZf38ZMMXf399gVLyijhWgi2vLli3F1uhRkroQ/ygVSf39999HpVIZLA0aNLB3WELo6argTXUVs4S/vz8jRozQV9snJCQwb8MGri5aBBs2QFAQXLxI0gsvsLl/f64fPFhg97+oqChOnz7NpEmTuHbtmkVfNIo6VoBuvPy+fftK0hWiBJSKpA7QuHFjrl+/rl/+97//2Tsk8RAyNdGKpZPDWGL9+vXExMQQuXIl9OsHZ8/ChAksV6v568oVVnfrRtWVKxnzxhtGr4SfeuopLl26hFqtZvv27QUeKyEhgZUrV+pfX1HHCtCNl+/h4WHXsed1jF3Vy5W+KEtKTVJ3dnYmICBAv+iGwBTClnK3YDf23Pfff69NsGZeuapUKry9vYuU7Pv370/9+vX/SayenjBnDoO2bqVeYCDDsrNxmjaNoB49UH33Xb6+7Vu2bKFGjRpkZGRw//59EhMTTR5r48aNBoPWWDNWgK2G9wXzZquz5HiS1EVZV2qS+l9//UVQUBC1atXilVde4cqVKwVun5mZSWpqqsEiHF9xXvGae/zcfc1zJwu1Wk1ERAQNGzY0+8rV2OsxlYRMVeH7+/vnGwdeo9Hg1qQJ3b/4gsqrV6MEBeFy5QpO/fvj+sQTVLl0SV9OaGgojRo1om7duhw5coRly5aZTJS2aChXWJW9uUlY19f/+PHjBX6JssdwwkKUhLw5LDMzs9B9SkVSb9++PatWreKXX35h6dKlxMbG8thjj5GWlmZyn5kzZxr0m69WrVoJRiysZe+knpuxZGHtKHe5X5epJGTJa9e1DL93/z4PXniBnD//JHnkSJRy5VDv3UuDgQNpMGMGXL5MUFAQYWFh+lkUz507Z7K2IW9DOWvGgjdWZZ+7xsLcJJy7r39BX6KKeovAljULQthStWrVDPLYzJkzC93H/p+cZujVqxf9+vWjWbNm9OzZky1btpCcnMyGDRtM7jNp0iRSUlL0S1xcXAlGLEqr3InVlmPPW1quxdX2FSqQPG4cOadPo3ntNRSViiq//opb8+Y4T56MU1oar732Gp06deLf//63YXU+/4wCl3dceGuSurEvPsZef2hoaIHV3uHh4TRv3py5c+cW+CXKki9axhK4XOkLRxUXF2eQxyZNmlToPqUiqedVsWJF6tWrx/nz501u4+bmhpeXl8EiRGFyJx9Lr8pz73vjxg2++uoro/evzSlXrVbj5eWln1XNbMHBaFau5Ormzdx/5BFUmZk4f/IJbfr1o8mPP/LkY4+xe/dubXe5PNX5ycnJNm/MplaruX//PosWLdInUt3rDwgIKDCpBwUFMWrUKPz9/c26723ONsYSuD0mDhLCHHlzWO7h2U0plUn97t27XLhwoeRm0hLCDLmTelRUFLGxsQZTqFraIKso87lnNW5M/JdfkrNpE5pGjXC+e5fARYs4P24cl377jQ1r1+bbRzc1qy0Tu1qt1g+mU5Qr4YLeO937DhT6/hpL4KV94iAhcisVSX38+PHs3r2bS5cu8ccff/D888/j5OTEgAED7B2aEEbp+mfnHvfdnKRuq/7uKpUK74oVUfXuTdbBg5ybMoXMatUYnplJu3PnGLdhA6oVK7SD2eTZz9ZDslp6JWzpl5/cSb2wWwWSwEVZVyqS+tWrVxkwYAD169fnpZdeonLlyuzfv9+qITiFKAm5q44toVar8fHxwcfHp0iJ1aDRnZMTN3v0YM9nnzGjbVsG+/oScvMmTm++CfXqwdKl8PeMbnm79BXUxc9clibS3Eld17YACk/Ypu7/F/YlwZEaZwpRVKXiLF63bh3x8fFkZmZy9epV1q1bR+3ate0dlihmD+OHbe5pSwt73YmJiXz99ddcv37drLK/++EHop2cmPPCCyS99x6Kvz9cugRvvolTvXp4f/45qvR0g33UajUeHh4W3dvP3R4gOzvb6j7gub/g6Bq4Xb161eT28fHxREZG5mvFLkldPEzkLBYO62H8sM2dgK5evcrHH3/MjRs3jG67atUqfv/9dxYtWmRW2S+88AI1atSgb//+pAwaRM65c7BwIQQHo7p+Hd9Zs6jVtSuVFi+G5GT9foqikJKSYjQpGkuYeZO6te0C8rZROH/+fIGzRC5cuJDffvuNhQsXmlVmUciANcJRPTyflkI4CHMTS2RkJGfOnCmw66Yl5fv5+fGvf/0LPz8/7QoPDxg5Ei5cIOezz8iqXh2n5GR85s/HqVYtGD+e+EOHWLhwIQkJCUaPZU1ysyaxhoaGUqdOHf2gOJYcN3c3NmuTet7jSVIXjkqSuhAlzNzEkjuR5W04p6tibt26NSNHjjQoM2/5uvviJhvfubqihIURt20b1+fOJbthQ1RpafDJJ6zs2JFLGzbw7YIFZr02cxsDFvT6jfUlDwoKIjQ0FA8PDzQajdHjjBo1iu7duzNq1CiD8mzRDz3v8WTAGuGoJKkLUUzyDiBjaSIICgpiyJAh1K9fHzDsrhUZGcm1a9fw8PAgMDAwX6LMnYTMaewWHx/Pwk8/JbZ9e4iOhp9+gm7dGKYotLpxg7c2bMDp8cfhm28gO9tkOba4gjWVhHVd7kzdpw8KCiI8PDxfg7yi9kM39ppkwBrhqCSpC1FM8ibaoiSCvPe1w8PDqVevHj169NBvkzv55E3q3t7eODs7m2xVHxkZSUxMDOvXrwe1Gk2vXiR/8w1VDh1ibN++BLq4oNq3TztLXN26MHs23LxpVuy5W7AnJyfrk/LVq1eNfskpLAlb2n+/qN3YNBoN586dM4hVBqwRjkqSuhAlQKPR0K9fP6PDs1pzZRsUFMSIESMIDg5GrVaj0Wi4ffs2t2/fzleW7suFi4sL/v7+uLi45BuCNjw8nODgYO7evUt8fLw+rpwmTbg5ezaX9+xB8+674OurbTE/cSIEB6N+7TXcDh7MNzOcseODYVI39SUndxLWxXHt2jVWrlxp8t5+cTI2sUxwcDBTp04lODi4xOMRoiCS1IWwgLn3w41Vh5crV46RI0caHZ7V2i5f5cqV0yf1gq5edVfrarWaSpUq5btiDwoKonz58ly7ds1oTUJOlSpo3n8frlyBqCho1w6yslCvW0dg//44tWwJixdDSgpgeKvB1BeO0NDQQq92de9PVFQUly9f1o/QZ2xKVmdnZ4KCgnB2ds732s35m5n6WxibWOZh7JkhSgc5I4UwwdgHt7VJ3drjV6hQAWdnZzw9PUlPT7f6XrUuHmdnZ5Nx5Zu33Zhy5WDwYDT79pG2axc5YWFoypVD9eefMHIkTiEheI8eTdTUqZw9c4aFCxeyYMECYmJi8o0tb6pa3NggMmFhYdSoUYP09HSOHTtmdEpWZ2dnQkJCbJ7UzZ1YRghHIEldCBPscTWWt/uVLql7eXlx7949q0ZUM5exedvB+FWxRqPhVvXqXH73Xd4aPJi4adOgUSOu37vHwh9+4PmvvqLtvn2o9u/n3LFjBmPgWyN3TcLEiRPNmpLVHLkTuamGjNaODiiEPUhSF8KB5L7PbGlreUVR9Em9qF2ucu+/fv16kxOyREVFcTY2ls9zcuDUKZYOGsShwEB+cHNjemoqEw4c4JHff2f077/jtWlTvhHrjDHVWl9XkzBr1iybXTnnTuoFzXNfqVIlKlWqJNXtwuHJGSpECTD3Kjp3q2pLWsvHx8cza9YsFi9erB8uNe++ltQ86PaPiooqsFre4L64SkXolCkEPv00rxw8yM358/Hv1o0pKhW1o6MJmDiRGu3bk9i3L5+OHEl8bKxZr0v35aJKlSqEhYXRvHnzYrlyNtWi3Vbj8QtREuQMFaIEaDQa0tPT8fT0LDAx5L7PnDvJGBvXPPeMapGRkezZs4d9+/axcuVKoy3tcyf1whro6Y4dGhpqslo+b7y638PDwwmsVYv0Pn148NNPXN69m6SxY8msXh11RgZRW7Zw4aef+KJFC7xHj0bZskU/W1zeuHJ/uchdE1EcbD2Dm4w6J+xBkroQJcTYlXJB1eS5k4yxsc9zV1OHh4fz+OOP06FDBwYNGmS0pX1uhSUc3bGDg4MN+pgXdk8/JiaGzz//3GCSmXgnJ95OTmZg+/a89sQTPNa7N40rVGBYVhaVfvgBlz59ICAAXn8dZds2zv35p/49yf3lojgU58hwktSFPUhSF6IEmJpBzNwq9txDxupmcvP29tbf5w0KCmLixImMGDGCwMBAi+MzNY977j7muQe/MZYMdfelz58/z6pVq/Td7NatW8eePXvYtm0bPx48yJc+Pgw6ehTvbdu4+vzzKH5+cPs2LF+O01NPsa1HDy5t2MCykSMJKF9e/+WiYsWKeHp64uzsbJM558Hw/bd1w0gZSlbYgyR1IUqAqVnGrBmZTJfUK1asaHCfNyEhgZUrV5o9FWtulia03NXiufXr14/atWszePBg/Sh4/fr14/HHHyc4OJgKFSroDsiDdu04MWQI6efOwY4d2iv1gABez8qi1Y0bjPjhB1yqVsUvPBynqCgy4uKIiooiMTHRZve4c7//tk7qMpSssAdJ6kKUgLyzjOmYex/XnKlHN27cyIULF1i1apUtQi6QqWpxf39/wsLCqFq1qv5q2t/fn4kTJ7JkyRJ69OjB4MGDDQtzcoKuXWHZMnIuX0azYQMjw8PxrV4dVWYmHv/9L05vvMHKNm2IjYpi+YABcOAAFQtpn2AOW99Hz02GkhX2IEldiBJgarIRc5n6UqCjVqsJCwujTp06+ZOmFfJetebtq24qGeoa7+kGuXFycsLLywuVSoWfnx8DBw6kfPny+jHpdSPi6VyNj+eDHTs40Lcvl7dvJ/PoUZLHj0dp3ZphQMuUFN784w/Ujz4K/v7wyivw1Vdmj0NfkmQoWWEPktSFKAGm7qmbq7AvBWq1mgYNGjBkyBCz7qkbu99b0PStBfVVz13m4sWLSU9PN5j2NXc5N2/eZMmSJdy4ccNoUs89h7wCKA0bkhwRQc7+/fhcuMDLEyfi2bMnirc33LoFa9fCa69pE3y7djB5MprffiP5+nWDBmr2aLQmQ8kKe5CzTYgSYE71ua1cv36dRYsWcf/+fZMJxZJ+7Gq1mqFDh5rsw63bx2CmN/KP/56ens727du5ePEiGzduJDo6moiICI4dO6YvT1cj0bVrVxYvXsyxY8f++fIRFETiM89wfeFCcq5fhz17YNIklBYttBPKHDoEH36I+okn8K5RA7p1gw8/hH370GRmGk3qknhFWSNnshAloLDqc52CrpZ1TLWq1nVxW716NTExMWzcuNFksrLkfq+uFmDMmDH5agpyxxgeHk79+vXp378/kL+P+f379xkwYAB16tQhLCyMf//731y5coWJEyfqywsKCqJ3794sWrSIEydOMHHiRIMvH0lJSdoBdm7ehMceg48+IufQIeIOHCDn88/htde45u/P9KwsbuzaBZMnw6OP4uTnh19YGKp58+DAAcjKyhe/EGWBnMlClABz76mbk9QLGs60QoUKhIWFFZqwi2OgFQ8PD0aNGkW9evX0SV7XmC4xMZE1a9Zw69YtALy8vJg9ezYhISHMmjXLoKwffviBcuXKkZGRwaxZs/Sv5caNG3zwwQf5JnIByPHzQxk0CL74guXDhnGgUyc+69ULXnwRKlfm+t27/GfnThLefhseeQQqVoT/+z9t0v/lF/3scqZemyVV99I/XdiTc+GbCCHsLe9AM6mpqSaTti0Ttq77XIUKFQq8otX1SQf02+WO47PPPiM2Npb//Oc/eHh4EBUVxbvvvsuSJUto0aKFvpz4+Hju3btHrVq1GDp0KK1bt6Zhw4Z4eXkxbdo0ypcvn28iF13NxdixYwkJCSF8yBBS09IIHzsWQkJAo2Hp669zeOdONGo1027d0vaL37VLuwCoVNC0KXTqBB07au/P164NKpU+SRsbj97Ue2HJ9nnfay8vrwLf58K2EQ83OSuEKAGmJimxZCpX3f6mknZxVCXrEhRgULYlx1Kr1URERNCwYUOmT59e4G2IqKgorl+/Tvny5fH3989XvW9sIpe8NRf5Wp2r1YROmkTw//0fg375BW7eRHPqFPfmz0cZOBCldm3tPfkTJ+DTT7Ut6uvWBV9f6NmT62PHsuDNN7l69KhZ75mp2yOFXcFfvXqVadOmcfXqVZNlSy2AKIxcqQtRAgpL6ubsn3vc9oK2yc7OtkXIZsVjrqCgIMaMGYOvry8VK1Y0WYsQHh5OUlISL730ktEyRo0aZXSf3DUXhcamVqOpX5/E8uUJjogA4PqxYwTFxuK0bx/88QdER2uv5rdvZ8H27WwDnLZu5ePgYO1VfNu20KYNNG8OVaoYFJ/7S8a0adP06wu7gje1X255ayWMkav5h5v8xYUoBezVoMvcoU7j4+NZuHAhCQkJ+nXWxBwcHMz48eMJCAgwe5+AgADCw8ML3Kew3gc5VaqgvPAC6BrSpaXB4cPw6adoGjVC4+KCAnD1Knz3HUyaBE88AX5+ULUqPP00/PvfsGEDoZ07U6tmzXy1EYW9l+Y0XjRnlDpzrvhF2SVJXYgS4GitrM2Nx9yhTnXd2b7//nujVfS6x87OzoXea1apVFSsWDHf/OVqtRpPT0/S09Pz9UEvbPa28PBwGjVqpE+YuROs0WTr6gqtW8Pw4YzZupWnRoxg9OnTsHs3/Oc/8NJL2ip6lQri42HLFpg5E/r3p0a3bixctYqQ/v1h+HBYuhR27yZy0aIC30tz2kLYKvGLskuq34UoAZZWV9viWAUlTnPjKaxRXt7tIiIijB43960BXVLXtQ9wdnbOt62Pj4/R9RUrViQxMbHQuPMKDg5mwoQJ+vvseRNfQdXeISEhzJ07V/tLw4bw+OP/PHn3Lpw8qa2u/3u5duIEn9+/z9CDBwk6eFC/6esA5cox7M4d7Qh4jRpB48ban35+Zv3ddO0FvLy8TG5j7G8mVfIPD5WiKIq9gygJulmtUlJSCvyHEMLRleQHdHZ2NlevXiU4ODhfkrVmu8JiL+z5K1euMG/ePIN7ypcuXeI///kPEyZMoEaNGmbFl7scIF+ZRXmPp0yezIH//pdHAgKYVq8enDoFp09DbKzpnXx8tMm9fn1tDYBuqVMHypWz6PjGYjf2vgnHVZR8JUldCGGSrZN6UU2dOpWDBw/Srl07/VW1sXV55U10ueMF8sVelCRoal9NWhrpR49S/vJl1GfOaBP96dNw4YK29b0pwcGGib5uXahXD2rVAjc3s2Iy5z0SjqMo+Uqq34UQRWbrNgOmrpSNVS2bc4sg7+2GvPHmjd2cluimmKoiV3t64tm5c/4d7t+Hc+fgzz/hr7+0y7lz2p/JydrGeVevws6d+fcNCoIaNaBmTe3P3I+rVdO2DcD82yii9JMrdSGESfa6F2vqyt9YPMVRS3DlyhXmzp3LkCFDaNSokb4rYYm+F4qinbRGl+jzLmlpBe+vUkHVqig1a5IVGIimWjXc6tRBHRKivfqvWlVb7a9SFf9rERaR6nczSFIXovSwpPq7OJJtdnY2ly9fBqB69eo4OzuX2C0Gs+gSfmwsXLqkXXI/vnRJWwMAxAMrgKFAvnb17u7aBJ93qVr1n59+fto570WJkep3IUSZYkn1d3H0LFCr1dy/f5/PP/+ccePGERISYtbALyVGpdKOeOfrqx0IJy9FgcREuHSJz+bM4eipUygVKvB+UBBcu6atzr95EzIy4Px57WKKWq0dYMffHwICtEvux7l/lyt/u5OkLoRwOPa+B6xWq9m4caN+Dvlp06YV6T57iVOptInW35+hCxawZMkShkREaK++dTIytH3sr179J9HnXq5dg+vXQaOBhATtcuJEwcd1cdFe2QcEaL8I6L545F4qVzZ87OJSvO/FQ0aq34UQDscRqrrz3gIoqW5hV69eZcmSJURERPwzfr295ORAUhLcuKFdEhL+eZz399u3Cy3O6K0Ab+/8id/HBypV0s6mZ+xnpUraWwdltFZA7qmbQZK6EKWHIwyWkveLRUnFVJzdz4r1C0NWlrbKX5fkb93SfiHItbx/9Cgnbt+mmUbD+/fvF9yVrzCurgUnfi8v7eLpmf+n7rGHh0N+MZB76kKIMqUkR+ArLIa8w94Wt+K89bBixQqio6NZsWKF7W8huLr+09DOhCF/f6kYEhEBgYHaLnt5Ej83b8KdO9olOTn/z+RkbQ1CVtY/twWspVYbJnljid/TE8qX134BKF/e9OPc6+x4S0Gu1IUQwoEUZ42AQ1XtWxuPomiH5zWV9HVfCNLSIDXV+M+0NG1bgeLi7Fx48i9XzuSSCni/8UbZr35fsmQJ//nPf7hx4wbNmzdn0aJFtGvXzqx9JakLIYRjef/99zl69CitWrXi/fffL7kDKwrcu2c64eddl56u3T493fBx7nXp6doaBBtIBbyhbFe/r1+/nnHjxrFs2TLat2/P/Pnz6dmzJzExMfj5+dk7PCGEEBYaNmwYK1asYOjQoSV7YJXqnytoC6b5LVRWlukvAHm/DNy/b3pJS4MdO6x7aaXlSr19+/a0bduWxYsXA9oqqmrVqjFy5EgmTpxY6P5ypS6EEKI0KPMN5bKysjhy5AiTJk3Sr1Or1XTv3p19+/YZ3SczM5PMzEz97ykpKYD2zRJCCCEclS5P6fKWjpubG26FTOJTKpJ6UlISOTk5+Pv7G6z39/fn7NmzRveZOXOm0dad1apVK5YYhRBCCFvKOx7C1KlTC217UCqSujUmTZrEuHHj9L8nJydTvXp1rly5gre3tx0jM61t27YcOnTIocu3pgxL9jFn28K2Keh5Y8+lpqZSrVo14uLiHPbWjJwbRT83rHlOzg3HPzdssZ2jfWYoisKtW7fw8fEx6AFR2FU6lJKk7uvri5OTEwl5+iMmJCQQYKKRg6lqCm9vb4f953RycirW2GxRvjVlWLKPOdsWtk1Bzxf0nJeXl5wbJVxGSZ4b1j4Hcm448rlhi+0c8TPD2otP+wzVZCFXV1dat27NjlytATUaDTt27KBDhw52jMy2IiIiHL58a8qwZB9zti1sm4KeL+73uLjIuVH0c6Msnhcg54YttitL50apaf2+fv16Bg0axGeffUa7du2YP38+GzZs4OzZs/nutRsjrd+FKXJuCFPk3BDGOPJ5USqq3wH69+/PzZs3mTJlCjdu3KBFixb88ssvZiV00FbHT5061ax7EuLhIueGMEXODWGMI58XpeZKXQghhBAFKxX31IUQQghROEnqQgghRBkhSV0IIYQoIySpCyGEEGWEJHXgp59+on79+tStW5fPP//c3uEIB/L8889TqVIlXnzxRXuHIhxIXFwcXbp0oVGjRjRr1oyNGzfaOyThIJKTk2nTpg0tWrSgSZMmrFixokSP/9C3fs/OzqZRo0bs3LkTb29vWrduzR9//EHlypXtHZpwALt27SItLY3Vq1fzzTff2Dsc4SCuX79OQkICLVq04MaNG7Ru3Zpz585Rvnx5e4cm7CwnJ4fMzEw8PDxIT0+nSZMmHD58uMRyykN/pX7w4EEaN25M1apVqVChAr169WL79u32Dks4iC5duuDp6WnvMISDCQwMpEWLFgAEBATg6+vL7du37RuUcAhOTk54eHgA2tlCFUWhJK+dS31S37NnD88++yxBQUGoVCo2b96cb5slS5ZQo0YN3N3dad++PQcPHtQ/Fx8fT9WqVfW/V61alWvXrpVE6KKYFfXcEGWXLc+NI0eOkJOTIzNAlhG2ODeSk5Np3rw5wcHBTJgwAV9f3xKKvgwk9fT0dJo3b86SJUuMPr9+/XrGjRvH1KlTOXr0KM2bN6dnz54kJiaWcKSipMm5IUyx1blx+/ZtBg4cyPLly0sibFECbHFuVKxYkePHjxMbG8vatWvzTUZWrJQyBFA2bdpksK5du3ZKRESE/vecnBwlKChImTlzpqIoirJ3717lueee0z8/evRoZc2aNSUSryg51pwbOjt37lT69u1bEmEKO7D23MjIyFAee+wx5YsvviipUEUJK8rnhs7w4cOVjRs3FmeYBkr9lXpBsrKyOHLkCN27d9evU6vVdO/enX379gHQrl07Tp06xbVr17h79y5bt26lZ8+e9gpZlBBzzg3xcDLn3FAUhcGDB9O1a1dee+01e4UqSpg550ZCQgJpaWkApKSksGfPHurXr19iMZaaCV2skZSURE5OTr5JX/z9/Tl79iwAzs7OfPLJJ/zf//0fGo2Gt99+W1q+PwTMOTcAunfvzvHjx0lPTyc4OJiNGzeWqel+RX7mnBt79+5l/fr1NGvWTH/P9csvv6Rp06YlHa4oQeacG5cvX2bYsGH6BnIjR44s0fOiTCd1c/Xu3ZvevXvbOwzhgH777Td7hyAcUKdOndBoNPYOQzigdu3aER0dbbfjl+nqd19fX5ycnPI1UkhISCAgIMBOUQlHIOeGMEXODWFKaTg3ynRSd3V1pXXr1uzYsUO/TqPRsGPHDqlCfcjJuSFMkXNDmFIazo1SX/1+9+5dzp8/r/89NjaW6OhofHx8CAkJYdy4cQwaNIg2bdrQrl075s+fT3p6OqGhoXaMWpQEOTeEKXJuCFNK/blRYu3si8nOnTsVIN8yaNAg/TaLFi1SQkJCFFdXV6Vdu3bK/v377RewKDFybghT5NwQppT2c+OhH/tdCCGEKCvK9D11IYQQ4mEiSV0IIYQoIySpCyGEEGWEJHUhhBCijJCkLoQQQpQRktSFEEKIMkKSuhBCCFFGSFIXQgghyghJ6kIIIUQZIUldCCGEKCMkqQshhBBlhCR1IYQQJsXFxdGlSxcaNWpEs2bN2Lhxo71DEgWQCV2EEEKYdP36dRISEmjRogU3btygdevWnDt3jvLly9s7NGGEXKkLh9KlSxfGjBlj7zBKRFl7rfZ8Pbdu3cLPz49Lly7Z5fhlWWBgIC1atAAgICAAX19fbt++rX/+X//6F5988omdohN5SVIXVrl58ybDhw8nJCQENzc3AgIC6NmzJ3v37gVKR8JSqVQFLu+//75Z5RTnax08eDDPPfecRfvY673/7rvvmDFjhs1j6Ny5M2FhYQVu8+GHH9KnTx9q1KhRpGNZ834/TI4cOUJOTg7VqlXTr5s8eTIffvghKSkpdoxM6DjbOwBROvXt25esrCxWr15NrVq1SEhIYMeOHdy6dcveoZnt+vXr+sfr169nypQpxMTE6NdVqFDBHmGVWj4+PjYvU1EUjh07Rr9+/Uxuc+/ePSIjI9m2bZvNj18SsrKycHV1tWsMLVq0IDs7O9/67du3ExQUBMDt27cZOHAgK1asMNimSZMm1K5dm6+++oqIiIgSiVcUQBHCQnfu3FEAZdeuXUafHzRokAIYLLGxsYqiKEpGRoYycuRIpUqVKoqbm5vSsWNH5eDBg/p9O3furIwePVr/+08//aR4eXkpX331laIoipKTk6N89NFHSo0aNRR3d3elWbNmysaNGw2Ov3HjRqVJkyaKu7u74uPjo3Tr1k25e/duga8pKipK8fb2NvpcQTEX9Fq3bt2qdOzYUfH29lZ8fHyUp59+Wjl//rzJ12rqvezTp4/BPiNHjlQmTJigVKpUSfH391emTp1qsL2xeAp73wor15z3Vfd6TMWwevVqxcfHR8nIyDAot0+fPsqrr75q9PXHxMQogPLHH3+YfI82btyoVKlSxWBd586dlREjRiijR49WKlasqPj5+SnLly9X7t69qwwePFipUKGCUrt2bWXLli0m329r4lUU8/7uERERyujRo5XKlSsrXbp0URRFe27Pnj1bqV27tuLq6qpUq1ZN+eCDD/SvsaD3vrC/b0FlmyMjI0N57LHHlC+++MLo89OmTVM6depkdnmi+EhSFxZ78OCBUqFCBWXMmDH5PvAURVGSk5OVDh06KEOHDlWuX7+uXL9+XcnOzlYURVFGjRqlBAUFKVu2bFH+/PNPZdCgQUqlSpWUW7duKYpimOjWrFmjeHp6Kj/++KO+7A8++EBp0KCB8ssvvygXLlxQoqKiFDc3N/0XjPj4eMXZ2VmZO3euEhsbq5w4cUJZsmSJkpaWVuBrKiipFxRzQa/1m2++Ub799lvlr7/+Uo4dO6Y8++yzStOmTZWcnJx8r9UUY0ndy8tLef/995Vz584pq1evVlQqlbJ9+/YC3/vC3rfCyjXnfdW9HlMx3Lt3T/H29lY2bNig3ychIUFxdnZW/vvf/xp9/WvXrlWcnJyU9PR0k+/RqFGjlCeffNJgXefOnRVPT09lxowZyrlz55QZM2YoTk5OSq9evZTly5cr586dU4YPH65UrlzZoOzc77c18SqKeX/3ChUqKBMmTFDOnj2rnD17VlEURXn77beVSpUqKatWrVLOnz+v/P7778qKFSvMeu8L+/uaKtscGo1G+de//pXvS15uW7duVVxdXY1+HoiSJUldWOWbb75RKlWqpLi7uyuPPvqoMmnSJOX48eP6540lrLt37youLi7KmjVr9OuysrKUoKAgZc6cOQb7LV68WPH29jaoDcjIyFA8PDzyXbWFh4crAwYMUBRFUY4cOaIAyqVLlyx6PaaSuiUxF+bmzZsKoJw8edLs/Ywl9bxXRG3btlXeeecdg21yl2vO+1ZYuea8r7mPa+q1DR8+XOnVq5f+908++USpVauWotFojJY5fvx4pVGjRiaPqSjaK+ewsLB8seR+PdnZ2Ur58uWV1157Tb/u+vXrCqDs27dPvy7v+21pvMYY+7u3bNnSYJvU1FTFzc3NaKIt7L0v7O9bUNnm+P333xWVSqU0b95cv5w4ccJgm+PHj1v1fydsT+6pC6v07duXp59+mt9//539+/ezdetW5syZw+eff87gwYON7nPhwgUePHhAx44d9etcXFxo164dZ86c0a/75ptvSExMZO/evbRt21a//vz589y7d48nnnjCoNysrCxatmwJQPPmzenWrRtNmzalZ8+e9OjRgxdffJFKlSpZ9TrNjdmYv/76iylTpnDgwAGSkpLQaDQAXLlyhSZNmlgVD0CzZs0Mfg8MDCQxMdHk9ua8b4WVa6v3dejQobRt25Zr165RtWpVVq1axeDBg1GpVEa3P3r0KK1bty6wzPv37+Pu7p5vfe7X4+TkROXKlWnatKl+nb+/P0CB752l8YJ5f/e8r+nMmTNkZmbSrVu3fOUV9t4X9vctqGxzdOrUSf8aTClXrhygbd8g7EuSurCau7s7TzzxBE888QTvvfceQ4YMYerUqSaTurlatmzJ0aNHWblyJW3atNF/gN69exeAn3/+mapVqxrs4+bmBmg/vH/99Vf++OMPtm/fzqJFi3j33Xc5cOAANWvWLFJclnr22WepXr06K1asICgoCI1GQ5MmTcjKyipSuS4uLga/q1SqAj90zXnfCivXVu9ry5Ytad68OV988QU9evTgzz//5Oeffza5/dGjR3n22WcLLNPX15c7d+7kW2/s9eRepzuvCnrvLI0XzPu75+3jrUuKxhT23hf2901OTi4wXlvQdXGrUqVKsR9LFEy6tAmbadSoEenp6QC4urqSk5Nj8Hzt2rVxdXXVd3sDePDgAYcOHaJRo0YG2+3cuZPvv/+ekSNHGpTv5ubGlStXqFOnjsGSu4uNSqWiY8eOTJs2jWPHjuHq6sqmTZusek3mxGzstd66dYuYmBgmT55Mt27daNiwodHEUxzyxmPu+1YYS95XY++JzpAhQ1i1ahVRUVF0797dZAwXL14kOTm50Cv1li1bcvr0abNfh6XMjRes/7vXrVuXcuXKsWPHDqPPF/TeF/b3LaxsWzh16hTBwcH4+voW2zGEeeRKXVjs1q1b9OvXj7CwMJo1a4anpyeHDx9mzpw59OnTB4AaNWpw4MABLl26RIUKFfDx8aF8+fIMHz6cCRMm4OPjQ0hICHPmzOHevXuEh4cbHKNevXrs3LmTLl264OzszPz58/H09GT8+PGMHTsWjUZDp06dSElJYe/evXh5eTFo0CAOHDjAjh076NGjB35+fhw4cICbN2/SsGFDq16rOTEbe62VKlWicuXKLF++nMDAQK5cucLEiROL9sabyVg8hb1vhbH0fTUWg1qtvYZ4+eWXGT9+PCtWrOCLL74wecwjR46gUqn0A5+Y0rNnTyZNmsSdO3esvs1SEHPjBaz+u7u7u/POO+/w9ttv4+rqSseOHbl58yZ//vknTZo0KfC9N+f/wlTZef/vrPX777/To0cPm5QlisjeN/VF6ZORkaFMnDhRadWqleLt7a14eHgo9evXVyZPnqzcu3dPURRtV6RHHnlEKVeunEE3r/v37ysjR45UfH19zerSdvr0acXPz08ZN26coijalrjz589X6tevr7i4uChVqlRRevbsqezevVu/fc+ePfXdz+rVq6csWrSo0NdUUOv3wmI29Vp//fVXpWHDhoqbm5vSrFkzZdeuXQqgbNq0yehrNcZYQ7m8+/Tp00cZNGhQgfEU9r4VVq4572vuMky9Jzqvvfaa0e5iuU2cOFGpV69ege+PTrt27ZRly5YZjUWnevXqyrx58wzW5f576OLq27dvvvLNiVfH2r97Tk6O8sEHHyjVq1dXXFxclJCQEOWjjz4y670v7O9rqmxbuH//vuLt7W3Q4FDYj4z9LoQocd26daNx48YsXLjQJuX9/PPPTJgwgVOnTulrBKzx5JNPUqdOHRYvXmyw3tbxliVLly5l06ZNbN++3d6hCKT6XQhRgu7cucOuXbvYtWsXn376qc3Kffrpp/nrr7+4du2aRe0Ecse1d+9edu3axRtvvFHs8ZYlLi4uLFq0yN5hiL/JlboQosTUqFGDO3fu8N577zF+/Hh7h6P3/PPPc+jQIQYNGsQHH3ygbxnvqPEKYYokdSGEEKKMkC5tQgghRBkhSV0IIYQoIySpCyGEEGWEJHUhhBCijJCkLoQQQpQRktSFEEKIMkKSuhBCCFFGSFIXQgghyghJ6kIIIUQZIUldCCGEKCMkqQshhBBlhCR1IYQQooz4f7iLlWYmWTcQAAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 512x384 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "PvI_plot_Step2 = HtoH.PvI(CatP, CatdP, CatI, Iscale=[1.0,1500.0], Pscale=[0.0,30.0], showfit='true', weighted='false', errorbars='true', Source=FullName)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "59006181",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Power law index: -1.0027705327275422 ± 0.08046783662639521\n",
      "Coefficient: 5.2333296111371495e+22 ± 2.1357245803848217e+23\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAgAAAAF2CAYAAAAY6yC7AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/P9b71AAAACXBIWXMAAA9hAAAPYQGoP6dpAABtfUlEQVR4nO3deVhUZfvA8e8MqyCgiAIj4L6+Ki6plVbmbmWWZrQqiC1mbqSWLaJoUVaaGulbIVhZkr/MNrPFXMpccsHcVxSUxRQBcUFgzu+PeWdiZ2aYYQa4P9d1rpgzZ865Z2zmuc+zqhRFURBCCCFEnaK2dQBCCCGEqH6SAAghhBB1kCQAQgghRB0kCYAQQghRB0kCIIQQQtRBkgAIIYQQdZAkAEIIIUQdJAmAEEIIUQdJAiCEEELUQZIACFED5OXlMWfOHPLy8mwdihDCyqrr+24XCcCyZcvo0qULnp6eeHp6ctttt/Hjjz8anr9x4wYTJ06kUaNG1K9fn1GjRpGRkWHDiIWoXnl5ecydO1cSACHqgOr6vttFAhAQEMCbb77Jnj172L17N/3792fEiBEcOnQIgGnTpvHdd9+xZs0atmzZQmpqKiNHjjT5OjExMZYO3apsHa+1r2/p81vifFU5hzmvtfW/cW1W0z5bW8cr33frv9bW/8alKHaqYcOGyscff6xkZWUpTk5Oypo1awzPHTlyRAGU7du3m3TODh06WDpMq7J1vNa+vqXPb4nzVeUc5rzW2NdkZ2crgJKdnW3yNeoqW39/TGXreOX7bv3X2tv33dHG+UcphYWFrFmzhqtXr3LbbbexZ88e8vPzGThwoOGY9u3bExQUxPbt27n11lvLPE9eXl6x6hOtVkteXh7Z2dmoVCqrvw9LKCwsJCcnp9Ze39Lnt8T5qnIOc15r7Gv0x9jy/4eaxtbfH1PZOl75vlv/tca+Jjs7G4CsrKxi+11cXHBxcTHpmhWyanphgr///ltxd3dXHBwcFC8vL+WHH35QFEVRVq1apTg7O5c6vmfPnsrMmTPLPV9kZKQCyCabbLLJJlut2CIjIy1a7tpNDUC7du1ITEwkOzub//u//2Ps2LFs2bLF7PPNmjWLiIgIw+Ps7GyCgoJISUnB09PTEiELKysoKOD8+fM0bdoUR0e7+V+VgoICkpOTAQgKCrJabJZ6/yXPY6+fqxBCJycnh8DAQJKTk/Hy8jLst+jdP2A3335nZ2dat24NQI8ePfjrr79YvHgxISEh3Lx5k6ysLBo0aGA4PiMjAz8/v3LPV15ViX6kgbB/BQUF5OTk4OnpaVcFVUFBAR4eHgBWjc1S7//mzZuoVCrq16+Ps7Oz3X6uQojivLy8rFpe2cUogLLo2+x79OiBk5MTGzduNDx37NgxkpOTue2222wYoRA1g1arJTc3F61Wa+tQhBB2xC7S/1mzZjFs2DCCgoK4cuUKn3/+OZs3b+ann37Cy8uL8PBwIiIi8Pb2xtPTk0mTJnHbbbeV2wFQCCGEEBWziwTgwoULjBkzhrS0NLy8vOjSpQs//fQTgwYNAmDRokWo1WpGjRpFXl4eQ4YM4YMPPrBx1EIIIUTNZRcJQGxsbIXPu7q6EhMTY3+TKAghhBA1lN32ARBCCCGE9UgCIIQQQtRBkgAIIYQQdZAkAEIIIUQdJAmAEEIIUQdJAiCEEELUQZIACCGEEHWQJABCCCFEHSQJgBBCCFEHSQIghBBC1EGSAAghhBB1kCQAQgghRB0kCYAQQghRB0kCIIQQQtRBkgAIIYQQdZAkAEIIIUQdJAmAEEIIUQdJAiCEEHVQfHw8KpUKlUpl2NevXz9UKhWhoaG2C0xUG0kAhBCiltEX5GVt69atA6Bx48b07t2b3r172zZYCwoNDUWlUtGvXz+jjp80aRLBwcE4OjqiUqnw8/Mr87g9e/YwdOhQPD09cXNzo2/fvvz6668WjNw2HG0dgBBCCOtwdnamW7duxfZ5e3sDcO+993LvvffaIiy78emnn+Ls7Iy3tzf//PNPmcf8/fff3HnnnVy7dg0fHx88PT3Ztm0bQ4cOZf369QwePLiao7YcqQEQopZLTU0lNjaW1NRUW4ciqpm/vz87duwott15551A2U0AZRk7dixt2rTBw8MDZ2dnmjVrxuTJk8nJyTEco7/zbt68OStXriQwMBBPT0+mTp3K9evXmTp1Kl5eXgQFBbF8+fJi509NTWXcuHFoNBqcnZ1p2bIl8+bNo6CgwHCMvkZjzJgxREZG4u/vT8OGDXniiSe4cuUKgOHaAFu2bDG8t82bN5f73g4cOMCFCxe45557yj3m1Vdf5dq1azRv3pzTp09z5swZevfuTWFhIdOnT6/ws7N3kgAIUcvFxcVx8uRJ4uLibB2KqIG++eYbLl++TKtWrQgMDCQ5OZmlS5cSHh5e6ti0tDSee+45XFxcuHLlCosXL6ZHjx7ExcXh5eVFSkoKEydO5MiRIwBcunSJW2+9lbi4OHJzc+nQoQMpKSnMnj2bp59+utT5V69ezaJFi6hXrx5ZWVmsWrWKN998E4Bu3brh4+MDgIeHh6F5w9PTs9z3FhgYWOF7LygoMFT1Dx48GA8PDxwdHbn//vsBXQJRkxNrSQCEqOXCwsJo3bo1YWFhtg5FVLOzZ8+W6gNgqi1btnDx4kUSExM5deoUr7zyCgDr1q3jxo0bxY69efMmP//8M8ePHzcUridPnmT//v0cPnwYV1dXtFotW7ZsAeD9998nJSUFX19fTp06xf79+/m///s/QFdDcfLkyWLnd3V15ciRI5w8eZIePXoAsHHjRgC+/vprQ5NG9+7dDTUe3bt3N/k96128eJHr168D0KRJE8N+X19fw9/Jyclmn9/WpA+AELWcRqMhPDwcjUZj61BENSurD4Cpfv31Vx5//HFOnTpVrMAvKCjgn3/+KXYX3bBhQ/r06QNAUFAQKSkpdOrUiebNmwO6QjQ5OZmMjAwAdu3aBUBGRkaxAhZAURR27txJ69atDfv69+9P06ZNAWjfvj179uwxnKs6KYpS7de0BkkAhDBDRkYGCQkJvPDCCwQFBdk6HCHKpO8DYK5Vq1YZ2rn9/f0JDAzk4sWLnD59GoDCwsJixxetbnd0dCy1T18DUbIA9fDwoGPHjqWu7+bmVuxxgwYNSp3fmoWxj48P9erV4/r161y4cMGwv+jfNfn7L00AQpghISGBY8eOERsba+tQhLAaffLg4eFBUlISO3futGiv9549ewK6wnz16tWGavtffvmF5557jgcffNCk8+kThqtXr1okPkdHRwYMGADAzz//zJUrVygoKODbb78FoHPnzjW6Zk0SACHMEBISQrt27crsCGUpqampLF68uEZ3MhI1W5cuXQC4cuUKLVu2pGXLlnz55ZcWO//EiRNp2rQply9fpl27dnTt2pVWrVrRqFEjxo4da/L52rdvD8Du3bvp3Lkzt956q6ENvyz9+vWjdevWrF27FtC1+bdu3ZrWrVuzc+dOAObPn0+9evU4c+YMLVu2pHnz5uzcuRMHBwcWLFhgxru2H5IACGEGX19fJk+ebNXsPzY2lqNHj0otg7CZ8PBwIiIi8PHx4cqVK/Tr14+oqCiLnb9x48bs2LGDsLAwGjVqxKFDh7h+/Tp33HEHixYtMvl848aNY9SoUXh5eXHw4EF27txZqpmiqDNnznDq1CnDUMLCwkJOnTrFqVOnDIlDcHAwW7ZsYdCgQdy4cYNLly5x++23s379eoYOHWreG7cTKqW29GaoRE5ODl5eXmRnZ1c4LETYjlarJScnB09PT9RqNQUFBZw7d46AgABDe589KCgo4OzZswA0a9bMarElJyezaNEipk2bVqV2xhs3bpCYmEjXrl1xdXW1289VCKFTXeWV1AAIu6HVasnKykKr1do6FLug0WiYMmVKjW5jFELYL0kAhBBCiDpIEgAhTJSamsqSJUtsMv5YCCEsRRIAIUwUGxvLsWPHSEhIsHUoZpHRBUIIkARACJOFh4fTrl07QkJCbB2KWWR0gRACZCZAIUym0WiYPHmyrcMwW3h4ODk5OVadw0AIYf+kBkCIOkZGFwghQBIAIYQQok6SBEAIIYSog+wiAYiOjqZnz554eHjQpEkTHnjgAY4dO1bsmH79+pVa1/rZZ5+1UcRC2BeZRElYwtatWxk+fDgajQaVSsW6deuMel1MTAzNmzfH1dWV3r17G5b5NfUYSzDnOpW95sqVK0ydOpVmzZpRr149br/9dv766y+rxF+d7CIB2LJlCxMnTjSsApWfn8/gwYNLrej01FNPkZaWZthq+kIMQliKJADCEq5evUpwcDAxMTFGvyYhIYGIiAgiIyPZu3cvwcHBDBkypNiSucYcU55+/foRHx9vsVjMec348eP55Zdf+PTTTzlw4ACDBw9m4MCBnD9/3qi47JZihy5cuKAAypYtWwz77rrrLmXKlClmnzM7O1sBlOzsbAtEKKwhPz9fSUpKUvLz88t8bC/y8/OVkydPKidPnrRqbKa8/4qOvX79urJ9+3bl+vXrJp9XWEabNm2UW2+9Vbl27Zphn1arVXr37q289NJLNoysfIDy9ddfV3pcr169lIkTJxoeFxYWKhqNRomOjjbpmPLcddddSlxcnFExm3Odyl5z7do1xcHBQfn++++Lva579+7KK6+8YlRcpqqu8souagBKys7OBsDb27vY/lWrVuHj40OnTp2YNWsW165dK/cceXl55OTkFNuEEMIWEhIS2Lt3L9u2bTPsW7VqFWfPnuXll1827HvjjTeoX79+hVtycrIt3kKZbt68yZ49exg4cKBhn1qtZuDAgWzfvt3oY6orFnNeU1BQQGFhIa6ursVeW69ePf744w+LxV+WkmVYXl6eRc9vdwmAVqtl6tSp9OnTh06dOhn2P/bYY3z22Wds2rSJWbNm8emnn/LEE0+Ue57o6Gi8vLwMW2BgYHWEL4TdSU1NJTY2Vmb+s6Fu3brRtWtXjh49CsC1a9eYNWsW8+fPx8PDw3Dcs88+S2JiYoWbPQ3fvHjxIoWFhfj6+hbb7+vrS3p6utHHFFUyCfr999959tlnK02CTL2Osa/x8PDgtttuY968eaSmplJYWMhnn33G9u3bSUtLq+QTqprAwMBi5Vh0dLRFz293EwFNnDiRgwcPlsqsnn76acPfnTt3xt/fnwEDBnDq1ClatWpV6jyzZs0iIiLC8DgnJ0eSAFEnxcXFcfLkSeLi4pg3b56tw6mz2rZta+jcvGDBAnx8fAgLCyt2jLe3d6maT1O99NJLvPXWWxUec+TIEdq3b1+l61jLs88+y8MPP2x4/PjjjzNq1ChGjhxp2FfdSdCnn37KuHHjaNq0KQ4ODnTv3p1HH32UPXv2WPW6KSkpxZYDdnFxsej57SoBeP755/n+++/ZunUrAQEBFR7bu3dvAE6ePFlmAuDi4mLxD0uImigsLIz09PRShY2oXu3atWPr1q2cO3eOt99+mx9++AG1ungl7BtvvMEbb7xR4XkOHz5MUFBQuc+/8MILhIaGVniOli1bGh13RXx8fHBwcCi1MFZGRgZ+fn5GH1NUySSoXr16NGnShNatW1c5FnNf06pVK7Zs2cLVq1fJycnB39+fkJAQi32O5fH09CyWAFiaXTQBKIrC888/z9dff81vv/1GixYtKn1NYmIiAP7+/laOToiaTaPREB4ebldVx3WRvgbgpZdeYvDgwfTr16/UMZZoAmjcuDHt27evcHN2drbIe3J2dqZHjx5s3LjRsE+r1bJx40Zuu+02o4+prliq+hp3d3f8/f25fPkyP/30EyNGjLBY/DZh1S6GRpowYYLi5eWlbN68WUlLSzNs+h6zJ0+eVKKiopTdu3crSUlJyjfffKO0bNlSufPOO42+howCsH8le6efPXtWmTp1qnL27FkbR1acjAIQ5ti3b5+iUqkUFxcX5cSJE7YOp0xXrlxR9u3bp+zbt08BlIULFyr79u0r9h1cunSp0r9/f8Pj1atXKy4uLkp8fLxy+PBh5emnn1YaNGigpKenm3RM0RiKlgNlbQUFBWXGb8x1zIl/w4YNyo8//qicPn1a+fnnn5Xg4GCld+/eys2bN837oCtRXeWVXSQAQJmbfuhHcnKycueddyre3t6Ki4uL0rp1a2XGjBkmfTiSANi/kgXT7NmzlaFDhyqzZ8+2cWTFSQIgzHH16lVFpVIpL7zwgq1DKdemTZvK/C0eO3as4ZjIyEilWbNmxV63dOlSJSgoSHF2dlZ69eql7Nixo9S5jTlGf/7yygT9lpSUVO57qOw65sSfkJCgtGzZUnF2dlb8/PyUiRMnKllZWeXGUFXVVV6pFEVRqq++wXZycnLw8vIiOzvbqm0qwnwFBQWcO3eOgIAAHB0dSU5OZtGiRUybNq3CNs/qVlBQwNmzZwFo1qwZjo7W6UpT8vMw99gbN26QmJhI165dcXV1Nem8wnIyMzNp1KgR+/fvp0uXLrYOR9ix6iqv7KIPgBBlkVXrRG2yf/9+nJ2d6dChg61DEQKo4iiA/Px80tPTuXbtGo0bN67y8BUhhKit9u/fT8eOHXFycrJ1KEIAZtQAXLlyhWXLlnHXXXfh6elJ8+bN6dChA40bN6ZZs2Y89dRTtWKRBCGEsKSpU6eyb98+W4chhIFJCcDChQtp3rw5cXFxDBw4kHXr1pGYmMjx48fZvn07kZGRFBQUMHjwYIYOHcqJEyesFbcQQgghqsCkJoC//vqLrVu38p///KfM53v16sW4ceNYvnw5cXFx/P7777Rp08YigQohhBDCckxKAL744gujjnNxceHZZ581KyAhhBBCWJ9FxgBdvHiRnTt3UlhYSM+ePWV2PiGEEMLOVTkB+OqrrwgPD6dt27bk5+dz7NgxYmJiZN5xIYQQwo6ZPAogNze32OO5c+eya9cudu3axb59+1izZg2vvPKKxQIUQgghhOWZnAD06NGDb775xvDY0dGRCxcuGB5nZGRYbKEJIWoTrVZLVlYWWq3W1qEIIYTpTQA//fQTEydOJD4+npiYGBYvXkxISAiFhYUUFBSgVquJj4+3QqhC1Gz6BKB+/fqlloEVQojqZnIC0Lx5c3744Qe++OIL7rrrLiZPnszJkyc5efIkhYWFtG/fHldXV2vEKoQQQggLMfs25NFHH+Wvv/5i//799OvXD61Wa1hsRAghhBD2zaxRAOvXr+fIkSMEBwfz8ccfs2XLFh5//HGGDRtGVFQU9erVs3SclpOTA7IaoBBCiDrO5BqAF154gbCwMP766y+eeeYZ5s2bx1133cXevXtxdXWlW7du/Pjjj9aI1TKWL7d1BEIIIYTNmZwAxMfHs379elavXs1ff/3Fp59+CoCzszPz5s1j7dq1vPHGGxYP1GLefx+ys20dhRBCCGFTJicA7u7uJCUlAZCSklKqzb9jx478/vvvlonOGrKzYckSW0chhBBC2JTJCUB0dDRjxoxBo9Fw1113MW/ePGvEZV0LF0JWlq2jEEIIIWzG5ATg8ccfJyUlhW+++YYzZ84wYsQIa8RlPW3b6gp/qQUQQghRh5mUACQnJwPQqFEjevbsSYMGDSo8/vz582YHZjUvvaT7r9QCCCGEqMNMSgB69uzJM888w19//VXuMdnZ2Xz00Ud06tSJr776qsoBWtwDD0DHjrq+AIsX2zoaIYQQwiZMmgfg8OHDvP766wwaNAhXV1d69OiBRqPB1dWVy5cvc/jwYQ4dOkT37t1ZsGAB99xzj7XiNp+DA0RGQkgILFoEU6ZAJTUZQgghRG1jUg1Ao0aNWLhwIWlpabz//vu0adOGixcvcuLECUDXP2DPnj1s377dPgt/vYcegv/8R1cLsGiRraMRQgghqp1ZMwHWq1ePhx56iIceesjS8VQPtVpXC/Dww/DeezB1KjRsaOuohBBCiGpTd5ckGzUKOnfWTQ389tu2jkbUIGq1Gi8vL1Qqla1DEUIIs9XdBECthqgo3d+LF0N6um3jETWGWq2mQYMGsqSvEKJGq9u/YCNGQO/ecO0avP66raMRwirUajX169eXhEUIUUzd/kVQqUC/bsF//wv/m+JYiNpEEgAhRFnM+kU4e/Ys3333Hbt27bJ0PNWvf38YOBDy82HOHFtHI4QQQlQLkxOAL774grZt2zJixAhuu+02brnlFv755x9rxFZ99LUAn34Khw7ZNhYhhBCiGpicAMydO5fHHnuMo0eP8vPPPwPwkn563ZqqZ08YORIUBV591dbRCCGEEFZncgJw+vRpIiMjadu2LQMGDOCzzz5j9erV1oites2bpxsZsG4d1IamDSGEEKICJicABQUFuLm5GR63b98erVZLek0fRtexIzz5pO7vl1+2bSxCCCGElZnVCXDlypX8+eef5ObmAuDo6Mi1a9csGphNzJkDTk6wcSP8+qutoxFCCCGsxuQE4I477mD+/Pn07duXBg0a0KZNG27cuEFsbCybNm3iypUr1oizejRvDhMm6P6eORO0WpuGI4QQQliLyQnAli1byM7O5tixY3z22Wc8+OCD3HXXXSxbtowBAwbQsGFDOnToYI1Yq8drr4GnJ+zbB6tW2ToaISxOZjIUQoCZiwEBtGnThjZt2vDII48Y9iUlJbF792727dtnkeBswscHZs3Sba+8ols5sF49W0cl7ExNLkT1sQsh6jaL/nq1aNGC0aNH84Z+XH1NNWUKBAZCSopunQAhSqjJCYAQQoCdTAUcHR1Nz5498fDwoEmTJjzwwAMcO3as2DE3btxg4sSJNGrUiPr16zNq1CgyMjKsE1C9ejB/vj44qOkTHRlBq9WSlZWFVvo9CCFEnWAXCcCWLVuYOHEiO3bs4JdffiE/P5/Bgwdz9epVwzHTpk3ju+++Y82aNWzZsoXU1FRGjhxpvaCeeAK6dtUtFzxvnvWuYyckARBCiLrFLhKADRs2EBoayn/+8x+Cg4OJj48nOTmZPXv2AJCdnU1sbCwLFy6kf//+9OjRg7i4OP7880927NhR5jnz8vLIyckptplErYZ33tH9vWwZnDhRlbdYo0gyUPOkpqayePFiUlNTbR2KEMJCSpZheXl5Fj2/XSQAJWVnZwPg7e0NwJ49e8jPz2fgwIGGY9q3b09QUBDbt28v8xzR0dF4eXkZtsDAQNMDGTAAhg2DggJdp8A6wlYJQG0vxKz5/mJjYzl69CixsbEWP7cQwjYCAwOLlWPR0dEWPb/ZowAANm7cyMaNG7lw4UKpwmLFihVmnVOr1TJ16lT69OlDp06dAEhPT8fZ2blUz2VfX99yZyCcNWsWERERhsc5OTnmJQELFsBPP8FXX8Eff0DfvqafQxilaCE2d+7cKp9Pq9WSk5ODp6enXXTWs/T7Kyo8PJycnBzCw8NLPScdFoWomVJSUvD09DQ8dnFxsej5zU4A5s6dS1RUFLfccgv+/v6oVCqLBDRx4kQOHjzIH3/8UaXzuLi4WObD6tQJwsPho49g6lTdOgHyQ2oVFRVi5tDXZNSvX98uCj9Lv7+iNBoNU6ZMQaPRlHpOhv0JUTN5enoWSwAszewEYPny5cTHx/Okfv58C3j++ef5/vvv2bp1KwEBAYb9fn5+3Lx5k6ysrGI/ZBkZGfj5+Vns+uWaNw8SEmDPHli5EsLCrH/NOqiiQqw2qO3vTwhRs5h9W3Tz5k1uv/12iwShKArPP/88X3/9Nb/99hstWrQo9nyPHj1wcnJi48aNhn3Hjh0jOTmZ2267zSIxVMjXVzdDIOj6ApjaoVAIIYSwM2YnAOPHj+fzzz+3SBATJ07ks88+4/PPP8fDw4P09HTS09O5fv06AF5eXoSHhxMREcGmTZvYs2cPYWFh3Hbbbdx6660WiaFSkydDmzaQkQE1faIjIYQQdZ7ZTQA3btzgww8/5Ndff6VLly44OTkVe37hwoVGn2vZsmUA9OvXr9j+uLg4QkNDAVi0aBFqtZpRo0aRl5fHkCFD+OCDD8wN33TOzrBwIQwfDosWwVNPQatW1Xd9IYQQwoLMTgD+/vtvunbtCsDBgweLPWdqh0BFUSo9xtXVlZiYGGJiYkw6t0Xdey8MHgw//wzTp8PXX9suFiGEEKIKzE4ANm3aZMk4agaVSnf336ULrFsHv/0G/fvbOiohhBDCZLYfG1XTdOwIzz2n+3vKFMjPt208QgghhBmqNBFQVlYWsbGxHDlyBICOHTsSHh6Ol5eXRYKzW3PmwOefw8GDEBOjmx9ACCGEqEHMrgHYvXs3rVq1YtGiRWRmZpKZmcmiRYto1aoVe/futWSM9sfbG958U/f37NmQlmbbeEStJDP4CSGsyexflmnTpnH//fdz5swZ1q5dy9q1a0lKSuK+++5jal24Ix43Dnr1gitXYOZMW0dTK9X1ArCuv38hhHVVqQbgxRdfxNHx31YER0dHZs6cye7duy0SnF1Tq3XV/yoVfPYZbN1q64hqHSkAhRDCesz+ZfX09CQ5ObnU/pSUFDw8PKoUVI1xyy3wzDO6vydOlA6BQgghagyzE4CQkBDCw8NJSEggJSWFlJQUVq9ezfjx43n00UctGaN9e/11aNRI1yHw/fdtHY0QQghhFLNHAbzzzjuoVCrGjBlDQUEBAE5OTkyYMIE39R3k6gJ9h8CnnoLISAgJAVnsRQghhJ0zuwbA2dmZxYsXc/nyZRITE0lMTDSMBLD0msV2b9w46N1b1yFw2jRbRyOEEEJUqsq9q9zc3OjcuTOdO3fGzc3NEjHVPGo1LF8ODg7w5Zewfr2tIxJCCCEqZFITQEREBPPmzcPd3Z2IiIgKjzVlMaBaoWtX3YRA776rmynw0CFwd7d1VEIIIUSZTEoA9u3bR/7/errv27ev3ONMXQyo1pgzB9asgbNnISoK3nrL1hEJIYQQZTIpASi6ANDKlSsJCAgoNUZbURRSUlIsE11NU7++bm6A4cN1Swc/8QR07mzrqOosrVZLbm4uWq3W1qEIIYTdMbsPQIsWLbh48WKp/ZmZmbRo0aJKQdVo990Ho0ZBQQE8/TRI4WMzkgAIIUT5zE4AFEUpc39ubi6urq5mB1QrLF4MHh6wYwf897+2jkYIIYQoxeR5APSd/1QqFbNnzy7W87+wsJCdO3fStWtXiwVYIzVtCm+8AZMmwYsv6moFAgNtHZUQQghhYHICoO/8pygKBw4cwNnZ2fCcs7MzwcHBTJ8+3XIR1lQTJuiWDN6+Xff3d9/p1g0QQggh7IDJCYC+I2BYWBiLFy/G09PT4kHVCg4OEBurGx74ww+6ZODxx20dlRBCCAFUoQ9AXFycFP6V6dABZs/W/T1lCly4YNt4hM1otVqysrKkQ6IQwm6YnQBER0ezYsWKUvtXrFjBWzL+/V8zZ0JwMFy6pEsCRJ2k1WrJzs4ut/OsEEJUN7MTgP/+97+0b9++1P7//Oc/LF++vEpB1SpOTrBiha5JYPVq+PZbW0dkVXKnK4QQNYPZCUB6ejr+/v6l9jdu3Ji0tLQqBVXrdO8O+o6Rzz4LmZm2jceKJAEQQoiawewEIDAwkG3btpXav23bNjSyHG5pkZHQvj2kpcHkybaORgghRB1n8igAvaeeeoqpU6eSn59P//79Adi4cSMzZ87khRdesFiAtUa9ehAfD7ffDqtW6WYLfPBBW0clhBCijjI7AZgxYwaXLl3iueee4+bNmyiKQr169XjxxRd56aWXLBlj7dG7t65T4Jtv6poC7rgDfHxsHZUQQog6yOwmAJVKxVtvvcU///zDjh07+Pvvv8nMzGT27Nl1dzVAY8yZA5066YYETpxo62hENVKpVHh5eZVaQEsIIWyhyr9EycnJXLp0idOnT/PTTz/x7bff8m0t7+leJS4usHKlblTAl19CQoKtI6q1UlNTiY2NJTU11dahAKBWq2nQoIEkAEIIu2B2E8Dp06d58MEHOXDgACqVyjC+WX/3X1hYaJkIa6Pu3eGVVyAqSjdNcJ8+EBBg66hqnbi4OE6ePElcXBzz5s2zdThCCGFXzL4VmTJlCi1atODChQu4ublx8OBBtm7dyi233MLmzZstGGIt9eqrcMstcPkyjB0rywZbQVhYGK1btyYsLMzWoQgLkWGmQliO2QnA9u3biYqKwsfHB7VajYODA3379iU6OprJMsytck5OutEAbm7w22+wcKGtI6p1NBoN4eHhMiy1FpEEQAjLMTsBKCwsxMPDAwAfHx9DO2uzZs04duyYZaKr7dq2hcWLdX+//DL8b6VFIYQQwtrMTgA6derE/v37AejduzcLFixg27ZtREVF0bJlS4sFWOuFh+vmA8jPh8ceg2vXbB2REEKIOsDsBODVV181VMNFRUWRlJTEHXfcwfr161myZInFAqz1VCr46CPQaODoUbCDSZRSU1NZvHix3fSer62kOlsIYUtmJwBDhgxh5MiRALRu3ZqjR49y8eJFLly4YJgZUBipUSPd0ECA5cttsmBQ0UI/NjaWo0ePEhsbW+1x1CWSAAghbMmsBCA/P58BAwZw4sSJYvu9vb1lEiBzDRz4791/eLhuzYBqVLTQDw8Pp3379oSHh5t1LinYhBDC/pmVADg5OfH3339bOhbx+usQHAwXL0JoaLUODSxa6Gs0GqZMmWJ273lFUcjOzpYEQAgh7JjZTQBPPPGExaqIt27dyvDhw9FoNKhUKtatW1fs+dDQUFQqVbFt6NChFrm2XXFxgc8/B1dX+PlnWLq02i5d1UJfCCFEzWL2TIAFBQWsWLGCX3/9lR49euDu7l7s+YUmjGu/evUqwcHBjBs3ztCvoKShQ4cSFxdneOzi4mJe4PauY0d4913dOgEzZ0K/frpaASGEEMKCzE4ADh48SPfu3QE4fvx4sedM7QcwbNgwhg0bVuExLi4u+Pn5mRZkTTVhAvz4I3z/PYweDbt3g6enraMSQghRi5idAGzatMmScVRq8+bNNGnShIYNG9K/f3/mz59Po0aNyj0+Ly+PvLw8w+OcnJzqCNMyVCqIj4du3eDECRg/XrdokI06WGq1WnJycvD09JSFbIQQopqULLdcXFwsWvtt8q/56dOnDQv/VJehQ4fyySefsHHjRt566y22bNnCsGHDKlxwKDo6Gi8vL8MWGBhYjRFbQKNGutUCnZxgzRp4/32bhSK9+oUQovoFBgYWK8eio6Mten6TE4A2bdrwzz//GB6HhISQkZFh0aBKeuSRR7j//vvp3LkzDzzwAN9//z1//fVXhYsOzZo1i+zsbMOWkpJi1Rit4tZb4e23dX+/8ALs2mXbeIQQQlSblJSUYuXYrFmzLHp+kxOAknf/69ev5+rVqxYLyBgtW7bEx8eHkydPlnuMi4sLnp6exbYaafJkGDlSN1Xwww9DZqatIxJmkNkVhRCmKlmGWbrze41s0D137hyXLl3C39/f1qFYn0oFK1ZAq1Zw9qzdLB0szQKmsfbsimq1mgYNGkgfDSGE0Uz+tdCPwy+5rypyc3NJTEwkMTERgKSkJBITE0lOTiY3N5cZM2awY8cOzpw5w8aNGxkxYgStW7dmyJAhVbpujeHlpesH4OKiGxnwzju2jkgSABNVdXbFykgCIIQwlcmjABRFITQ01FAVcePGDZ599tlS8wCsXbvW6HPu3r2bu+++2/A4IiICgLFjx7Js2TL+/vtvVq5cSVZWFhqNhsGDBzNv3rzaOxdAWbp1gyVL4JlndEsH33or3HmnraMSRpKJloQQ9sbkBGDs2LHFHj/xxBNVDqJfv34Vjiz46aefqnyNWuGpp+D33+Gzz+CRR2DfPvD1tXVUtZYMfxRC1GYmJwBFZ+MT1UylgmXLYM8eOHJEN0nQr7+Cs7OtI7Nb+qaKJk2amFyI619bv359SQCEELWO/KrVNPXrw9q1upkBf/8dpk61dUR2S7+08YkTJ6Svgp2TPiVCVD9JAGqi9u1h1ap/awQ+/NDWEdmluLg4Tp06xZo1a2wdiqiEJABCVD9JAGqq++6D+fN1fz//PGzbZtt47FBYWBitWrVi9OjRtg5FCCHsjiQANdmsWbp+APn5MGoUnDtn64jsikajYdy4cfhKR0khhChFEoCaTKWCuDjo0gUyMuDBB+H6dVtHJYQQogYwezVAgI0bN7Jx40YuXLhQqu1uxYoVVQpMGMndHdatg549dcsGP/00fPKJzVYOFEIIUTOYXQMwd+5cBg8ezMaNG7l48SKXL18utolq1KKFbuVABwfdHAGLFtk6IiGEEHbO7BqA5cuXEx8fz5NPPmnJeIS5+veHd9/VDQucMQPattV1FKxBZDpbIYSoPmb/0t68eZPbb7/dkrGIqpo8GcLDdYsF6WcKrEEkARBCiOpj9i/t+PHj+fzzzy0Zi6gq/bwAAwbA1au6GgAZGSCEEKIMZjcB3Lhxgw8//JBff/2VLl264OTkVOz5hQsXVjk4YQYnJ/i//4M+feDwYRg+HLZuBQ+Parl8amoqS5YsISQkhGbNmlXLNYUQQpjO7ATg77//pmvXrgAcPHiw2HNVXR5YVFGDBvDDD9C7NyQm6poDvvkGHKs06MMosbGxHDt2jISEBHr16mX16wkhhDCP2SXCpk2bLBmHsLTmzeHbb6FfP1i/Xtc5cOlSqw8PDA8PJysri5CQEKteR9RNqampLF68mGnTphEUFGTrcISo0aS3VW3Wu7duWCBATAwsWWL1S2o0GiZPniyz7wmriI2N5ejRo8TGxto6FCFqvCrVCWdlZREbG8uRI0cA6NixI+Hh4Xh5eVkkOGEBo0bBggUwcyZMmwaBgTByZLVcWr/Ai6enp0169qvVaurXr1+qSUqr1ZKTk2OzuIT5wsPDycnJITw83NahCFHjmf3rt3v3blq1asWiRYvIzMwkMzOTRYsW0apVK/bu3WvJGEVVTZ8OzzwDigKPPQZbtlT5lPqq2NTU1HKPsfUKb/oEoGQhb+u49Iz5DEVxGo2GKVOmoNFobB2KEDWe2QnAtGnTuP/++zlz5gxr165l7dq1JCUlcd999zFV1qi3LyoVvP8+PPAA5OXB/ffD/v1VOqVUxVadfIZCCFuqUg3Aiy++iGORnuWOjo7MnDmT3bt3WyQ4YUGOjvDFF3DnnZCTA0OHQlKS2acLDw+nffv2ZVbFZmRksGTJErmzrURFn6EQQlib2QmAp6cnycnJpfanpKTgUU1jzoWJXF11wwE7d4b0dBg8GC5cMOtUFVXFJiQkcOzYMeLi4qoaca0m1dlCCFsyOwEICQkhPDychIQEUlJSSElJYfXq1YwfP55HH33UkjEKS2rQADZs0A0TPHkS7rkHrlyx6CVCQkJo164dYWFhFj1vTSbTHAsh7I3ZowDeeecdVCoVY8aMoaCgAAAnJycmTJjAm2++abEAhRVoNPDTT7rZAvfs0Y0KWLfOYqf39fVl8uTJaDQa0tPTLXbemkyfAAghhL0w+3bE2dmZxYsXc/nyZRITE0lMTDSMBHBxcbFkjMIa2rbVTRDk7g6//or68cchP9/WUQkhhKgmVZ4b1s3Njc6dO1siFlHdevbU9Qm4917U33yDj6LAl1+Cg4OtIxNCCGFlJiUAERERzJs3D3d3dyIiIio8VhYDqiEGDICvvkJ54AHqf/st2gkT0C5fbuuo7IJMOyuEqM1MSgD27dtH/v+qifdVsNa8LAZUw9x7L9pPP0X9+OOo4+LAzQ1eeMHqly06Ix9g9Ox81TWTX9Fx+nPnzrXadYQQwhZMSgCKLgC0cuVKAgICSv0AK4pCSkqKZaIT1UZ56CEunj9P4+nTUcfE0ECrhcWLAevdCetn5Ktfvz6A4W9jEgBjj60KmXZWCFGbmf3r2aJFCy5evFhqf2ZmJi1atKhSUMI2ro4aReH77wPQYNkyVK+/DtTsGeuqMimRjNMXQtRmZicAiqKUuT83NxdXV1ezAxK2pTzzDNp33gHAYe5ciIqq0TPW6SclqonJixBCWJPJowD0nf9UKhWzZ8/Gzc3N8FxhYSE7d+6ka9euFgtQVD/1Cy9AQQG89BJERtK0sLDG3gmHhISQkJBQI5MXIYSwJpMTAH3nP0VROHDgAM7OzobnnJ2dCQ4OZvr06ZaLUNjGiy/qFhF68UUcoqJokJ0N/6sZqEmKTkokhBDiXyYnAPqOgGFhYSxevNjQg1vUQjNngloNM2bQYPFitB4eEBWlSwyEEELUaGZPBKRf6OXw4cMkJydz8+bNYs/ff//9VYtM2Ifp0ylUFBxmzkQ9f75uX1SUbWOyQ9U1NFEIISzF7AQgKSmJBx54gAMHDqBSqQydAvVzABQWFlomQlEtKlqsRpk2jcycHLznz4f58+H6dYiOtkGU9qu6hiYKIYSlmP1LNXnyZFq0aMGFCxdwc3Pj0KFDbN26lVtuuYXNmzdbMERRHSpbrS4nPJxC/eyO776L+rnnwMZJnr7Q1Wq1No1DVJ1+rglzhmsKIcxjdgKwfft2oqKi8PHxQa1Wo1ar6du3L9HR0UyePNmSMQo7oUyaBLGxoFKh/vhjfF54waYLCEkCUHvU5LkmhKipzE4ACgsL8fDwAMDHx8eQuTdr1oxjx45ZJjphf8aNgy++QHF0pP4336B+5BG4ccPWUYkaribPNSFETWV2AtCpUyf2798PQO/evVmwYAHbtm0jKiqKli1bWixAYYdCQtB+9RWKszPqb7+F4cPhyhVbR2W0sqqbpTbBtmTWRSGqn9kJwKuvvmr4sYyKiiIpKYk77riD9evXs2TJEpPOtXXrVoYPH45Go0GlUrFu3bpizyuKwuzZs/H396devXoMHDiQEydOmBu6sADlnnvIiI9HcXeHX3+Fu++GjAxbh2WUsqqbJQEQQtQ1ZicAQ4YMYeTIkQC0bt2ao0ePcvHiRS5cuED//v1NOtfVq1cJDg4mJiamzOcXLFjAkiVLWL58OTt37sTd3Z0hQ4ZwQ6qebUatVuM6bBjKr7+Cjw/s2QN9+sDJk4ZjbN2xS61W4+XlVWp1SqluFkKIKiQAZfH29jZrKeBhw4Yxf/58HnzwwVLPKYrCe++9x6uvvsqIESPo0qULn3zyCampqaVqCkT1MYwauPVW2LYNWrSAU6dwuPNOXA4cAHRzRRjbscsayUJ5IxukulkIIaqQAPTv37/MNdIvX75scg1ARZKSkkhPT2fgwIGGfV5eXvTu3Zvt27eX+7q8vDxycnKKbcJK2raFP/+Ebt1Q/fMPmscfp96WLYSFhRl9p13beoFLk4IQoqpKlmF5eXkWPb/ZCcDmzZt5//33eeCBB7h69aph/82bN9myZYtFggNIT08HdHO6F+Xr62t4rizR0dF4eXkZtsDAQIvFJMrg5wdbtqAdOBD1tWv4P/MM3t9/z6RJk8q90y561x8WFkbLli0JCwur9FK2blowhiQAQoiqCgwMLFaORVt4ArYqNQH8+uuvpKenc+utt3LmzBkLhWQZs2bNIjs727ClpKTYOqTaz8MD7TffcOX++1EVFNBgyhS8PvgAylk6uuhdv5+fH08++SR+fn4VXkKr1RITE8ORI0cM01ELIURtlJKSUqwcmzVrlkXPX6UEwN/fny1bttC5c2d69uxplRkA9QVCRoke5hkZGRUWFi4uLnh6ehbbRDVwdubCO++Q+b9q/0b6WQPLmDDInM54Wq2WESNG0K5dO6NqC6xB7u6FENWhZBnm4uJi0fObnQDoO/u5uLjw+eefM2XKFIYOHcoHH3xgseAAWrRogZ+fHxs3bjTsy8nJYefOndx2220WvZawELWaizNncuHll1H+N2sgQ4fC5cvFDjO3M56tl/iVBEAIURuYvRiQUqJa99VXX6VDhw6MHTvW5HPl5uZyssjwsaSkJBITE/H29iYoKIipU6cyf/582rRpQ4sWLXjttdfQaDQ88MAD5oYvqkHW2LEUNmuG3wsvoPrtN7j1Vvj+e2jTxtahmazoan9CCFEbVGk1QB8fn2L7Ro0aRbt27dizZ49J59q9ezd333234XFERAQAY8eOJT4+npkzZ3L16lWefvppsrKy6Nu3Lxs2bMDV1dXc8EU1uTZgAIVbtuD44INw/Dj07g1ffaWbOKgGKbranxBC1AZmJwDNmjUrc3+nTp3o1KmTSefq169fqRqFolQqFVFRUUTJOvQ1U5cusGsXPPAA7NgBgwfDBx+AjdrwRc1V2aqVQgjjmZQAREREMG/ePNzd3Q136eVZqF86VggAX1/YtMmwmBBPP4360CF4/nlbR2YzUpiZTv+ZCSGqzqQEYN++feT/rzf3vn37yj3OnNkARc1UtG1cvyx0w4YNKSgoIDs7u/jBrq6wahV06ACzZ6NevBjf3bvh//4PqvijXjKOmkAKMyGELZmUAGzatAmA/Px81Go1y5cvp00N7NAlLKdo27g+AfD29qagoIArZa0QqFLBa69Bhw4ooaHU27YNpXdvCr/4AowouLVaLbm5ubi5uVUYhxBCiIqZ9Uvp5OTE33//belYRF3y0EMU/vEH+c2bo0pOxrl/f3x++KHSl+kTgIr6jFiKVNELIWozs3/ZnnjiiVozb7uwkU6dSPvmG7T33osqL4/W8+fjOHUq3Lxp68iA6ksAZF4BIYQtmD0KoKCggBUrVvDrr7/So0cP3N3diz0vnQCFMbSenmjXrqVw/nwc58/H8b//hb//1vULqEWr9VXUR0GaL4QQtmB2AnDw4EG6d+8OwPHjx4s9J50AhUnUagpffplTXl60mzcP1fbt0LUrfPaZbshgLSCFvBDC3pidAOg7BAphKVl9+nBz2zZcHn8c9u/XTR88axbMnQuOZv+vahT9CoPTpk0jKCio2l9vbTVxlIQQwrrkl0DYhfT0dD799FPS3N11kwU9+6xuFcE33oD+/eH8+VKvMWZZYGPb8YuuTGiOqr5ez5j+AOb0GZB+BkKIkqqcABw+fJgNGzbw7bffFtuEMEVcXBynT5/WLfHr6grLlkFCAnh4wO+/65oENmwo9ZrKCl1jEwBzVia05Ov1rJUACCFESWbXq54+fZoHH3yQAwcOoFKpDMOy9O3/hYWFlolQ1Anh4eHk5OQUL0Affhi6d9f9d98+GDYMx2nTUA0fjlarZdSoUVy5cqXKhS6YtjJhWdX95q5sKIQQtmJ2DcCUKVNo0aIFFy5cwM3NjUOHDrF161ZuueUWNm/ebMEQRV2g0WiYNGkSbm5uxe9sW7eGP/+EiRMBcFy0iC7PPEODtDTc3d2ZNGkSGo3GqOYAvfJqBSq6sy56flOq++VuXQhhr8xOALZv305UVBQ+Pj6GGeD69u1LdHQ0kydPtmSMopYqWRArikJ2dnbpwtLVFd5/H9atQ/Hxof6JE/jecw8NVq3S9RPAtDb48hKAcq//v/MfOXKEmJgYwsLCjK7ulwRACGGvzE4ACgsL8fDwAMDHx8dw59WsWTOOHTtmmeiEXbDWhDgmn3fECPJ27SKrd29UN27QeO5c/J56CjIyLNYGX57w8HDatWvHiBEj8PPzk+p+IUSNZ/YveqdOndi/fz8AvXv3ZsGCBWzbto2oqChatmxpsQCFKCpVUZjRqRNJr7yC1tkZ982bcejWjaY7d5pUKJt6Z+7n58e4ceNo3LhxVcIXQgi7YXYC8Oqrrxp+PKOiokhKSuKOO+5g/fr1LFmyxGIBCtsztxpbpVLh5eVl0ZqDuLg4Tp46xceKwvl168hr1w7VP//g8PDD+EyeDJcuGXUeU9+TqWsQmNInQQghbMHsX+YhQ4YwcuRIAFq3bs3Ro0e5ePEiFy5coH///hYLUNRMqampLF26lGvXrlUpAShZUIeFhdG6dWtCQ0O52bYt59auRfviiyhqNfW/+w6HLl3g668t9TbMpu+TEBcXZ+tQLE76NQhRO1i0Udfb21umARZANU6M4+KCdv58Cv/4g5tt2qC6cAFGjoTHHjO6NsAa9H0SwsLCbBaDtUgCIETtYNI8ABEREUYfK4sB1W1ljuu3gLi4OE6ePEl8fDzjxo3794mePUn77jsCV6xA/c478MUX8NtvsHw5PPCARWMwhn5eAD8/P9LT0w01IpVNFWzMlML2Pu2wEKJmMKkGYN++fUZtiYmJVgpX2Jvy2rqtNTFO0SaAkhQXF7Svvw7bt0OHDpCRAQ8+qJtIKC3NonGYyphZC6H8mpOid92Wql0RQtRtJtUAyAJAoqSihdHcuXOtfj2NRkN4eDj+/v5kZGSUfVCvXrB3L8yZA++8A2vWwE8/wZtvwjPPgBF9Eiy9eE5YWBi5ubmV1oiUV3NSdDVBa9Wu2JK1hpoKIcpXpW9bVlYW7777LuPHj2f8+PEsWrSI7OxsS8UmagBrj783m6urrsD/6y/o2RNycuC556BPHzhwoNKXW7qdu6wakbKuYUzNSW2cdlgSACGqn9nftt27d9OqVSsWLVpEZmYmmZmZLFy4kFatWrF3715LxijsmN0XRt266ZoElizRLSy0YwdK9+6oZs1Cdf26TUPTJwDnzp2TIYNCiGpndgIwbdo07r//fs6cOcPatWtZu3YtSUlJ3HfffUydOtWCIQpRuQrvIB0cYNIkOHwY7YMPoioowOGdd9AMGYLqu+8M0wmDbXq4G9s/AKQHvhDCcqpUA/Diiy/i6PhvNwJHR0dmzpzJ7t27LRKcEOVJS0tjyZIlhn4ARlUhBwSg/fJLMj78ECUgAKeUFBxGjoR77oHjx8nIyGDp0qUcP34crVZrKGwLCgqsWuiaurZAdnZ2qQmJKksMavLERJL0CGEdZicAnp6eJCcnl9qfkpJiWCNACGOp1Wq8vLyMnkciPj6eY8eOkZCQYPK1rg8aRN6+fWQ99xyKszNs2IBD1658O3kyxw4dMpxTvzhQQUEBOTk5KIpilcLIEs0olcVVk0cOSAIghHWYnQCEhIQQHh5OQkICKSkppKSksHr1asaPH8+jjz5qyRhFLVTyR93UTmChoaG0a9eOkJAQ8wKoX5+sGTMoTEyEe+9FlZ/P5D17uHXHDsY1aFCsWQAgIyODDz74gLfeeosTJ06UOyVwyfdlL3fedttZUwhhM2YnAO+88w4jR45kzJgxNG/enObNmxMaGspDDz3EW2+9ZckYRS1U1bs6f39/Jk+ejK+vb9UCadMGvv+ewnXraBwUxNxr1wiOisK5Xz9c9+wxHLZmzRp27drF1q1bWbNmTbmnK/m+9O37tp4S2NhRCEKIusPsBMDZ2ZnFixdz+fJlEhMTSUxMJDMzk0WLFuHi4mLJGEUdYqs7aNXw4WT/+ScXIyLQurmh3rWLpiEhNJk4kSv79vHwww/Tr18/7rrrLh5++GGjFznSt+/b45TA5iyIJAmDELWHyQmAVqvlrbfeok+fPvTs2ZOoqChat25N586dcXNzs0aMwsaqqxDWX0ffCQ+qr+1arVbTwNeX7OeeI/nXXykYNw5Frcbjp5/wufNO/vPRR7z8zDO89NJL+Pv7G91cYY1hktZYZdEYkgAIUbuY/Avy+uuv8/LLL1O/fn2aNm3K4sWLmThxojViE3aiugrh2NjYUh37bNF2XdikCQUxMZz7/nuu3nknqoICGn7yCS0GDaJhbCyqvDyTz6lfTtiYwrOsgrZoEmbtSXOkoBeibjD5F+STTz7hgw8+4KeffmLdunV89913rFq1Sn4sarHqKoTDw8NLdeyz5URDN9u25fxHH5H88cfktm6Nw5UrNF6wgMBBg1DFxUFBgdHnMjYBUKvVeHp6kpOTU+zY6lxeWBIAIeoGkxOA5ORk7rnnHsPjgQMHolKpbN7LWVhPdRXCGo3GMh37KqC/k963b5/RzRpXb7+dvR9+SOobb5Dv64tTaioOTz8NHTvqVh20YEFZ3t19yeWFbVlI28vIBiFE1ZicABQUFODq6lpsn5OTE/n5+RYLSghjmNMWrr+Tfumll4xu1rhw4QLxn37Kidtv58zPP3Nx1iwUHx84cQIeewwlOBjVN9+UGjpYlFqtpn79+sViNaUgLZmE2TIBMGXmQiGE/TI5AVAUhdDQUEaOHGnYbty4wbPPPltsnxDWULQgNactXH8n/eabb1bYrJGammqYafDLL7/kxIkTfPnllyiurmSHh1N4/DiFc+ei9fBAdfAgDg89RMCoUTj88kuZiUBZCYCp1fr6WQDLKvTNuSs/d+4cb7/9NufOnTP6NWDazIVCCPtlcgIwduxYmjRpgpeXl2F74okn0Gg0xfYJYQ2Ojo4EBAQUm4LaFPo76W7dulXYrBEXF8exY8dYs2YNo0aNIjAwkFGjRpGRkcGSJUtIvXIF5eWXOff772hfegnc3XHctw+n++/H/4EHUH//faVNAyWr9StTWFjI+fPnKSwsLPWcOR01Y2NjOXz4MB9//LFJtQl2vwCUEMIoJv+K2npCE1G36e/6C0zogGeOsLAwcnJyeOihh2jYsCGPPPIITZo04cMPP+T48ePExsby2muvofXyQjtvHupp0+DNN1GWL8fl779h9Gj827fn+rRpaO++m9TUVJYuXcq0adMICgoC/i1I/fz8SE9PLxWDvprf09MT0NW+nT17lk2bNjF16tRitQnh4eHk5OSYdFceFhZGeno6Y8eOJSsrq1QNhRCidpNvu6hxqqP9W98hEWDZsmVcunQJgNGjR9OuXbvSBW2TJrBwIXlHj3L+ySdRPDxwOXqUBs88Q8DQocRPm8bRw4dNukMv631+9913HDt2rFQibs5duUajITw8HH9/f6NfI4SoPWpEAjBnzhxUKlWxrX379rYOS9hIeSviVZVareb69essXbqU1NRUMjIymD59OgcOHOCXX34BwNfXl8mTJ5db0KYWFDDHxYWTGzeSOXkyhV5eOJ8+zbPffkvPHTsY7+wM166ZHePw4cNp166dxWcWtOehfzLqQAjrqBEJAMB//vMf0tLSDNsff/xh65BEFdjjj7parWbNmjWGO+yEhARcXV25ceMGgwYNMuoccXFxnDx5kvh167g8eTJJv/3GpRkz8G/cmKicHAJffRWaNYPISLhwweQYfXx8eP755/Hz8zN6YiFj6Fc+tMcEoCavZCiEPasxCYCjoyN+fn6GzcfHx9YhCUxfxU+vqj/q1poOt2jHvJCQEIKDg5k+fTq//PILF0oU2GUlMWFhYbRu3Zpx48bh5eWF4uFB1jPPkHf0KJdmz0Zp3hwuXoSoKBxatqTRyy+jOn68zFhKnj8jI4O4uDgyMjJMmlnQUlJTU3nvvfc4evRotV5XVjIUwjpqTAJw4sQJNBoNLVu25PHHHyc5ObnC4/Py8sjJySm2CcsrLwGoLDEo70fd0Mu+kpoBa02HW7Qt3dfXl+eff57169fz+++/s3r1alQqleHYspIYjUZDWFgY7u7ueHp6/hufmxtXwsIoPHIEvvwSevZElZeHxxdf4BIcTNMJE2DrVrIuX+bcuXMsXryYJUuWFDv/mjVrOHXqFGvWrCEtLY3Y2Ngq1aCkpqYSGxtLYmKiYcijfn9ZtTP697tixQqLN79UREYdiLqqZBmWZ8Y05BWpEQlA7969iY+PZ8OGDSxbtoykpCTuuOMOrly5Uu5roqOjiw1LDAwMrMaIRWUFdHk/6gkJCRw7dswuq3tdXFxo2LChIQkoL4lRFKXUVL4Gjo4wejTs3EnBpk1cGzQIRaXCacMGHAcMoN7tt7Ni4kSOHjwIUOz8o0ePplWrVowePZr4+HhOnjxZpVE58fHxnDlzhpdffrnYGgzl1c6Eh4cTEBBAbm6uIVkQQlhPYGBgsXIsOjraouc3bzB1NRs2bJjh7y5dutC7d2+aNWvGl19+WW614KxZs4iIiDA8zsnJkSSgBggJCSEhIcHw72rsnb5WqyUnJ6f4XbeFjBkzhvz8fCZPnkyDBg3Izc0FjLszValUZcekUkHfvlz48EMCrl1DvXgxfPopLkeP8uzRo6idnRkXGIh6zBj8NBrUajVt2rRh3LhxNGnShGHDhrFhw4Zi03Kbavz48Vy5coWHH36Y1atXG9ZgKG9IoUajwd3dnaNHj5KQkMAtt9xiGKYowweFsLyUlBTDMGDQ3YRYUo381jZo0IC2bdty8uTJco9xcXHB09Oz2CbsX8le9qYkAMb0Yjen6cDX15dx48bRtGnTSo9NTU1lxYoVhjtk/fUcHR3Lv27btmhjYjj7xx9cfPll/Fq2ZM7NmwTFxRFw992ohw+H9etJT0019AH48ccf8fT0ZP369Ua9h7Led9FJkfSJV2pqarHEpujnqtVque+++2jbti0hISGVfub2PLJAiJqgZBkmCQCQm5vLqVOnZPyyMJk5CYBarcbV1dWo18TFxRna6Y25brEC1suL7HHjuHngABnx8WiHDUNRqVBv2IB6+HDW3n03qb/9xv/FxREaGkrr1q2NHg5Y2fsur+mlZALg5ubGc889Z9SCTZIACGHfakQCMH36dLZs2cKZM2f4888/efDBB3FwcODRRx+1dWh1grXXn69NwsLCDO30lSlv6V/Uaq7fdRfab79FdeIESkQEWi8vns3J4dYzZ5gZH0+LKVOY3rYtmiZNDOeqyr9RSEhI2RMc2YgkD0JYX434RT937hyPPvoo7dq14+GHH6ZRo0bs2LGDxo0b2zq0OsHeEoDqnENA34av7/hnTOfGcePGGXWHbNTn2qoVhW+9xZk//oDoaKZ26UJTRcHxxx9pN3MmLm3awMsvoz592uR/o6JDCSub4EhfIAN4eHgUGw1hDZIACGF99vGLXonVq1eTmppKXl4e586dY/Xq1bRq1crWYQkLU6vVeHl5VVq4WHJimMoKGn1M+oLVlGTIkomKUq8eOQ88wKm4OM5s2EDe5MncbNAAVXo6REdDmzbQvz+sXAkVjI4pSp8AAOV+7vr3cO7cOU6cOMGKFSu4fv16lZJBUwp3SQSEsJ4akQCIusHYwtWSE8NUpYApGm9Z51mxYgUHDx7k448/Nup85c2BkJqaSnR0NAsWLODq1asUtm5N3rx5/JmQQN7nn8OwYbpRBZs2QWgo+PrCY4/Bjz+CEYsmVfS5F12yWD8PQXx8fIXns2TnQGtN+yyEkARA1EDVOTGMVqvlypUruLm5VTjZUVmF2tixY2nevDmhoaGVXic1NZWIiAj2799famx/bGwsmzdvZuvWrYbOhWq1GldPT5QHH4T16+HMGZg3D6VtW7h+Hb74Au65B5o2halTYfduKFGIpqen8+mnn5a5EqFe0ZkR9fMQjBs3rsJETe7ahagZJAEQdYKphZK+6h90o040Gg2OjpVPm6FWq7l27Rrvv/8+oCtAjUlU4uLiDOsOlOzZHx4eTr9+/bjjjju477778PDwKF34BgXBq69SePAgqevWoZ04ERo31q03sHgx9OwJHTvC66/DiROGa54+fbrCyYRKzoyoHw5ZUQJgiaYPe1wrQojaRhIAUStUVmCUTAAqa24o+rwp7f5qtZrvv/+e48ePV1pVXjTue+65h+DgYBYuXFgqYdBoNMyaNYuZM2eWnmK4JJWKm8HBaN97D86fh++/h5AQcHWFo0fh1VehbVvo1o3x16/T0dfXqKYUfU3ItWvXKk2iijYbVHbOspIyrVZLTEwMR44cqdJMh0KIikkCIGoFUzsGVjQuv6qL7OiH1BlT9a+Pe/369UyePNkweqCyZEOtVlO/fv1SxxVLhJyc4N57YfVqyMiAFStg8GBwcIDERALffpt3V64k8P77Ub/1Ft6XLlVYrX/16tVibfHlJV1Fmw3KO1dubi4FBQXlJgAjRoywyrLHQoh/SQIgagVLdQwsmgCUHJVgbDOCfkidMRNV6eMODQ0tlngYM6NhWQlAyUTIEHP9+hAWBj/9BOnp8NFHaAcNQnFwQLV/P+pXX8XzlltQ9+iB+s038U5LQ11kVEBaWhrx8fHk5uYarhkbG8uRI0eIiYkpFm9AQACRkZGo1eoyEwT9Z3z+/Plya20qG5YohKg6SQCETVh6bgFzOwZW1HSgb89funQpqampRiUAxg5lLBm3r68vubm5ZfZ2LzoUsaxCv2hcJROhMmP28YHx42HDBnKOHUP73/8WqxlQv/YanrffjrpdOxxmzMDn4EFW/q+/wM8//2y4fnh4OO3atWP48OFkZmaWal6Ji4srVStTNJ74+Hijam3MWfpZOiIKUTlJAIRNlNeDvronHSp5x1yykLVU00JV6M+ZmZnJihUrSvXaL/r5mZIIqdVqvFq1Qv3007qagYwM+Phj3egBFxc4dQrVe+9Rf/hwJnzwAb3OnePppk3hyhWysrLw8/Nj8uTJNGnShOzs7FKFbVm1MlqtlhMnThAbG8uwYcPKrLU5d+4c77zzDmlpacXevyQAQliWJADC5qojAaisvbro6oNFEwBjmxaqo9f6F198waFDh/j444/Nul6lhWKjRhAeDj/8ABcvwldfwdix0KgRTbOzWXDqFN1nz8ZJo8FlxAhUS5filJRUanihXlnJSGpqKtOnT+fw4cP8+OOPZSYr+qaFb775hvT0dBkNIISVSAIg6oTy7uQru2M29o66rPOnpaURGxtrscJr9OjRNGvWjNDQ0HLfT8kEqmiiYNJdcf36MHIkxMdDejp5v/zCnhEjmO3lRdrNm9TbuhWn6dMJGjSI5oMG4RMZier77+F/MwuWJzY2FldXV/Ly8sqdTyAsLIw2bdowYcKEMpsRhBCWIQmAqHHMqSWw9p28vj189OjRhgI2Pj6ekydPmjyUrbwY9OPwNRpNsfdT9PiSn01FTRhGJwSOjih9+/JBo0bs7NWLZc8+S+ZLL1HYrx/nHR2JTknh6qpVODz4oK4WYeBAeOcdOHiwVO1AeHg4wcHBREZGlttJUv/+Sr5PIYRlSQIgapyyEoDKhu9V5U7eGBqNhkmTJlGvXj1DDMYs2Vv0veg7u1XWOU6/LO+kSZPQaDQVxlx0SF7Jwr6sBKC8pCA1NZVr167RNCCAsBkzyHnmGfJ//JG3n3qKnV27srRzZ5TmzeHmTdi4EWbMwLFbNwJ79UL95JO6vgVJSWg0Gp5//nl8fHwqnAfg+vXrpfo0SLu+EJYlCYCoFSwxfh8su86Av79/pTMB6hOACxcusGrVKq5fv8748eMrjKHk/PgVxawvQJs0acK+fft47733KqzdKK+QjYuLIy0tjfr16wMYahwefuIJ/Lp1459evUjZvBmOHdPNPDhsGIqrKw4XL6JevRqeegpatsShTRu006ax+pVXSN+/v8wYzp8/T2xsLOfPnzcqtrLITIJCVE4SAGFzlvixLm+IXEXKKlAqqikwNk59oe7o6Gh0TEWn5S0vBpVKVeb5jKndOH/+PC+++CL79++vtHajrM8lPDycDh068PDDDxva5ePi4vD19aV+/fqcO3eO2BUrdLMMTp4M69dTePEiaV98gfaVV6BPH3B0RHXmDAnr15O6ZQvL+/ZlRY8enH3ySbRffKGbuRBYuXIlycnJrFy5stLPrbzYP/74Y6v2HZDaCFEbSAIgbM4Sy/uamwBUttJc0Sp6Y+M0JwGorOYhNTWVN998k6VLl1a4eE954uPjcXFx4caNG4ZrlJfQKIpSalhfQEAA06dPx9/fn7CwsGIz/Y0ePZp27dqVjt3Fhbxbb0U7Zw788Qdcvgw//kjY+PG0b9QIBTiUmclnq1ejfuwxCAiAli155tAhuisK4XfeWe4Ig/LoC+bQ0FCr9h2QBEDUBpWvbiKElYWHh5OTk2OTjl4qlarCufX1hTlYN87K7uJjY2PZsmULoGtaMGaa4aJCQ0NJT08nNDTUcI2iCc3cuXNNjtXPz4/z588bP2tf/fowdChN+vVjeHg4jRwc+HDOHB52d0c5eRLV/v2QlERgUhLvAQwZgtKoEU26d0c1cCDK7bfrFjgyMcasrKyK108Qoo6SBEDYXHUs71veyAFTRhRU5zLEJYWHh3Pu3DlUKhVhYWHFai3K6xSZk5ND/fr1adCgAa6urkyePBk3N7di56xKQqNWq2nYsCEFBQXklhj+V9Edsr62pmnLljy7ZAkAhc2a4XjtGuzYQfK337Js7Vqez8yk6aVLuP3yC/zyCwDNnJ0hOBhuvRV694bevdG2aEHOlSt4enoaajUmTZpUbJIpU2uHrEX/7yIJibAH8n+gqBOqY4bBjIwMlixZYpWOZxqNhpdffplZs2YBsGTJEjIyMoDyE4CsrCwAQyfDFStW8M8//xQ7Z1kT9RQ9d0XUajXe3t54e3uXOT1xeQVduc01np4weDAfNWzIjnbtWB4RQd7mzZx9/nkK77sPpVEjVDdvovrrL1i6FJ54Atq0QeXri/ODD6LMnUvsrFkc/ftvVqxYYZFOoZYmTQfCnkgCIGoFczoSWrKnuFqt5ptvvuHYsWNWn7QmLi6OY8eO8c033xid0MTGxhr1Gv1xCQkJpZ4zdU5+Y5Zcvn79umGtBb3w8HA6duxI+NNPo/TuTdrjj5O/Zg0pe/bwQmgoZxYtIicsDKV3bxQXF1SXLuG2eTMOUVE88/nn3Prbbzy1ciV+U6fi8M47uP7+u25mQwsy9/8dGZ0g7IkkAKJWMKcjYVmvMbemQK1WM3HiRDp06FBsWmFzz1XR68LCwujQoQMTJ040+tz6XvyVvSY8PJyAgAByc3NLFVKmvB9jCjq1Ws2aNWtKJU0BAQHMmDGDgICAYrUFsStWcDg9ndhLl8icPZvCP/6gMDOTc2vX8s/s2Wgfewz/Nm2IBJzOnSNm40YORUby3zFjuODvD0FBMGIEzJ0L334L586Z3MlQz9yOq5bo8CqEpUgCIGoFc8bvl/WaqjQVlKxSt1YCUFbVfcmq5ZIFcFmvKdrBsehiQu7u7oax+ObEB8YXdJX9u6WnpxMbG0t6errh2Hvuueff9+bsTF6XLuSMGYN25UoKDx9mx48/MqpdO7Y2b06Ejw+Jbm58BJCSoiv458zRJQKBgWgbN0YZPBhmzoRPP4V9++DGjTJjqWjVRWPJzIbCnkgCIGoFczro2bJTn6WVTACMKYDLSgAAQkJCDMP6SvYlUKvVpVZyLKut39iCzpjRD0ePHuXjjz82zH64fv16w3tLTU3ljTfeIDo62jAV8tebNuHs68tVX1+iv/2W5k8/TdjBg7B1K7z3HowZA507ozg4oL50CdUvv8Dbb+v2d+8O7u7Qvj2MHg1RUbB2LZw4gTY/36xVF4sKCAggMjKSgIAAk14nhDXIKABhc9W9BLC5rBmnMeeuaMii/o5/2rRpBAUFVamHf9FhfefOnWPRokVMmzaNgICAUoW9PgHQV9frVVRAFn2v+hEE5XWK07+P0NBQsrOzAV0TSG5uLuHh4cWGRzZt2pS5c+cyYcIErly5wrPPPkv79u3x9/dHExAAjo5wxx2Gc6ccP86SyEimdO5MYGoqHDig2y5f1s1oeOwY/N//GY53qFcP/9atUXftitKxI/UaNYLbb9dNfuTkZNRnWzTpsiQZXSDMIQmAsDlr/ShamjXjrOzc+h73np6eZRaWJcf0W6p2o+h5X3vtNYsMqSv6XvUFV2ZmJk2aNClVw6B/H1qtliVLlhASEkL37t0N763o8Eh9suPn58eYMWNo27Ztxc0Uq1ZxKCuLj/PymPv++7qdigJpabpE4ODBf/976BBp16/z0YEDPHXgABrAV38iR0do3VpXa6Df2rXT/bea/r8umqgFBQVVyzVFzScJgBA1QMlCs6Sq3PEXrT0omTBUdt6SNQ/mXHvp0qXUq1ePOXPmEBQUVOZ71Y98SEhIoHv37oamCT8/P15++WUAk5OdMt+bSgUajW4bMuTf/YWF/HfKFPZs24bWy4s5Gg2Fhw/jcOIEqmvX4OhR3VaSr++/CUGrVrpEoVUr3fa/dRX077Mqd/DmTuok6jZJAESdZekq/epqyijrOpXd8VcUW9HCIzIyEi8vL3Jyckqdt6zCvryCx9jPIi4uju3bt+Po6Fhm4aU/j76wDgkJQa1W4+npaZjoqKT09HQ++eQT3NzcaN++fbnX1rfHe3p6VhgjAA4OhM+cSc6iRYRPm4YqKEj346nV6tYwOHpU12SgTwSOHtXtz8jQbf9rpijG19eQFCgtWpDfoAHanj1Rt2unW1ZZpao8rv956qmnuHnzJk899ZTRrxFCpVQ0EXotkpOTg5eXF9nZ2cZ94UWNUtfbQAsKCjh37hwBAQE4OjqWelyR5OTkYtXHRT9LwPD33Llz2bVrF7169TIU1CVfa6ozZ84we/ZsHB0dDTUAZSnZGVFf5T1p0iTDrIjNmjXD0dGRyMhIdu7cSbt27Xj77be5du2aRf6/qOwzLfX/4JUrcPz4v8nBqVO67eRJuHSp2GtTgY+ApwANgJeXLjlo1uzfrXnzf/9u2NCkBEHULNVVXkkCIEQtULLwMSUhMjZZKKuwNyXRKO/aycnJeHp6ljmjYHkiIyPZtWsXPXv2ZOrUqcC/oxSSk5OJiopCURQiIyMt1iZe2Xs1KRnKzv43GTh1ijlffkni2bN0LShgzpUrlQfj4VF+ctCsGTRpAlZKhOt6sl0dqqu8kiYAIWqBkp0IrdFh0RrDJot2bjSlMNE3CYwfPx5vb+9Scbq7uxsmGLJUm3hlzRomtcN7eemGHHbvDsC4xx9n0aJFjJs2Tbfg0enTuu3sWd125sy/f1+4oKtdOHhQt5XFyUnXj6FpU90WEFD6b40GXFxM/hykw2HtITUAQtRxxt7RlXWcuTUAVb2LrOj1Wq2Ww4cP8/HHHxMREVFthVRVmkNM+hyvXYPk5LKTg7NndX0PjP1Zb9y4dILg51d8a9KkWKKgr30p2hQkLEuaACxMEgBhj8wtCO2lGtbcOKradGCtuGx1TVM+j3PnzhETE8PEiRPLnlAoP183lPH8ed127lzZf+flGR9gw4aGhOCchwcxKSk8N2QIge3b6zoz6pOFxo3BwcGk9y5KkwTAwiQBEPbI3ILQ2gWotVW182BtY0ryYJE7cEXRdUQsmRScPw/p6bqRC/r/5ucbf16VCry9wcdHlwwU/W95f7u7S4fGEqQPgBB1gLnj6Ks6/t7WZNx6cab02ajKnA8GKtW/BXFwcPnHKYpuZsT0dMOmTU/n5tmzuFy+jOrChX+f++cf3bDIS5d027FjxsXi4lI6MfD21tU6NGxY/t/16kniUEVSAyCEDZl7N1fT22GlBsB89tL8U0phoW7Z5YsXdclARX//849uM6UZoiQXl3+TgYoSBS8v8PQs/V9nZ8u9dwuTGgAh6gBz7+YschdoQ7VpIabqZrdTZzs46PoD+PpWfizoaheuXePc33/zwX//y3MDBhAAusTg8mXdlpkJly+TmpHBR8nJPKXVosnJ0SUbeXn/1j6Yw9W1eEJQVpJQ8r/6rX593ebuDm5uNbYmQmoAhLAhe+1EZ212exdbh1XaudBKjKnNmjNnDnv37qV79+7MiYyE3FxDclA0USjzcU6Obt4F/X+vXrXsG1CpdImAPikouVX0XMlj/pdQ5BQW4uXnJzUAQtRm5t7N1ZQVFMtjt3exddhHH31EYmIiH330UbU2KxlTmzV+/HhiYmIYP368rsD18Ph3MiRTFRTo5lEomRgY+9/cXN2mTyQU5d99NYzUAAghhLBZDYClaoOqPX6tVjcnQ9GEQP93WVtlz+fm6s6Xl0cO4AUyDLComJgY3n77bdLT0wkODmbp0qX06tXLqNdKAiCEELVXsWaCOXNsHY75CgvJSU/HKyDA6uVVjak/TEhIICIigsjISPbu3UtwcDBDhgzhwoULtg5NCCGEjT399NP06NGDp59+2tahVI2Dg655oxrUmBqA3r1707NnT95//31AV20UGBjIpEmTeOmll0odn5eXR16RISbZ2dkEBQWRkpIiNQBCCCHsVk5ODoGBgSQnJ+Pl5WXY7+LigosZ6zeUS6kB8vLyFAcHB+Xrr78utn/MmDHK/fffX+ZrIiMjFUA22WSTTTbZasUWGRlp0bK1RowCuHjxIoWFhfiWGF/q6+vL0aNHy3zNrFmziIiIMDzWarX06NGDvXv3oqohYzZ79uzJX3/9VWuvb+nzW+J8VTmHOa819jX6OwKpwTKerb8/prJ1vPJ9t/5rjX2Nvsb67NmzxUbLWPTun1o8DLCsqhIXF5di1Sn2zsHBwaY/9ta+vqXPb4nzVeUc5rzW1Nd4enpKAmAkW39/TGXreOX7bv3XmvqaBg0aSCdAHx8fHBwcyMjIKLY/IyMDPz8/o88zceJES4dmVbaO19rXt/T5LXG+qpzDnNfa+t+4Nqtpn62t45Xvu/Vfa+t/45JqVCfAXr16sXTpUkBXpR8UFMTzzz9fZidAIWoTGcYqRN0hawGUEBERwdixY7nlllvo1asX7733HlevXiUsLMzWoQlhdS4uLkRGRlq8DVAIYX+q6/teY2oAAN5//33DREBdu3ZlyZIl9O7d29ZhCSGEEDVOjUoAhBBCCGEZNaIToBBCCCEsSxIAIYQQog6SBEAIIYSogyQBEEIIIeogSQCEqGFSUlLo168fHTt2pEuXLqxZs6bC/UKImqu873VWVha33HILXbt2pVOnTnz00Ucmn1tGAQhRw6SlpZGRkUHXrl1JT0+nR48eHD9+nJycnDL3u7u72zpkIYSZyvu+u7q6kpeXh5ubG1evXqVTp07s3r2bRo0aGX3uGjMRkBBCx9/fH39/fwD8/Pzw8fEhMzOTwMDAMvdLAiBEzVXR993NzQ2AvLw8FEXB1Pt5aQIQoppt3bqV4cOHo9FoUKlUrFu3rtQxMTExNG/eHFdXV3r37s2uXbvKPNeePXsoLCwkMDDQqP1CiOplze97VlYWwcHBBAQEMGPGDHx8fEyKTRIAIarZ1atXCQ4OJiYmpsznExISiIiIIDIykr179xIcHMyQIUO4cOFCseMyMzMZM2YMH374oVH7hRDVz5rf9wYNGrB//36SkpL4/PPPSy2YVylFCGEzgPL1118X29erVy9l4sSJhseFhYWKRqNRoqOjDftu3Lih3HHHHconn3xS7LXl7RdC2J6lv+9FTZgwQVmzZo1J8UgNgBB25ObNm+zZs4eBAwca9qnVagYOHMj27dsBUBSF0NBQ+vfvz5NPPmk4rrz9Qgj7VJXve0ZGBleuXAEgOzubrVu30q5dO5OuL50AhbAjFy9epLCwEF9f32L7fX19OXr0KADbtm0jISGBLl26GNoTP/30U7Kzs8vc37lz5+p8C0III1Xl+379+nWefvppQ+e/SZMmmfxdlwRAiBqmb9++aLXaMp8rb78Qomaq6PuemJhYpXNLE4AQdsTHxwcHB4dSnXkyMjLw8/OzUVRCCGuw9fddEgAh7IizszM9evRg48aNhn1arZaNGzdy22232TAyIYSl2fr7Lk0AQlSz3NxcTp48aXiclJREYmIi3t7eBAUFERERwdixY7nlllvo1asX7733HlevXiUsLMyGUQshzGHX33eTxgwIIaps06ZNClBqGzt2rOGYpUuXKkFBQYqzs7PSq1cvZceOHbYLWAhhNnv+vstaAEIIIUQdJH0AhBBCiDpIEgAhhBCiDpIEQAghhKiDJAEQQggh6iBJAIQQQog6SBIAIYQQog6SBEAIIYSogyQBEEIIIeogSQCEEEKIOkgSACGEEKIOkgRACCGEqIMkARBCCDuQkpJCv3796NixI126dGHNmjW2DknUcrIYkBBC2IG0tDQyMjLo2rUr6enp9OjRg+PHj+Pu7m7r0EQtJTUAQpShX79+TJ061dZh2L3a8jldunSJJk2acObMGZvF4O/vT9euXQHw8/PDx8eHzMxMw/OPPPII7777ro2iE7WRJACi1klPT2fSpEm0bNkSFxcXAgMDGT58OBs3brR1aBYRGhqKSqVCpVLh5OSEr68vgwYNYsWKFWi12mqNZe3atcybN8/w2JIJwV133YVKpeKLL74otn/p0qVoNBqLXEPv9ddfZ8SIETRv3tyi5zXXnj17KCwsJDAw0LDv1Vdf5fXXXyc7O9uGkYnaRBIAUaucOXOGHj168Ntvv/H2229z4MABNmzYwN13383EiRNtHZ7FDB06lLS0NM6cOcOPP/7I3XffzZQpU7jvvvsoKCiotji8vb3x8PCw+HkVRWHfvn34+/vz1VdfFXtuz549dO/e3WLXunbtGrGxsYSHh1vsnOXp2rUrnTp1KrWlpqYajsnMzGTMmDF8+OGHxV7bqVMnWrVqxWeffWb1OEUdoQhRiwwbNkxp2rSpkpubW+q5y5cvK4qiKDdu3FAmTZqkNG7cWHFxcVH69Omj7Nq1q9ixd911lzJlyhTD42bNmimLFi0qdkxwcLASGRlpOP75559XpkyZojRo0EBp0qSJ8uGHHyq5ublKaGioUr9+faVVq1bK+vXrS11n0qRJyowZM5SGDRsqvr6+hnOWZ+zYscqIESNK7d+4caMCKB999JFhX2FhofLGG28ozZs3V1xdXZUuXbooa9asMen6a9asUTp16qS4uroq3t7eyoABAwyfb9HPaezYsQpQbEtKSlJWrlypeHt7Kzdu3Ch23hEjRihPPPFEme/x2LFjCqDExMQobm5uytWrVw3Pde7cWZk9e3aFn5Ep1qxZozRu3LjU/sLCQuWtt95SWrVqpTg7OyuBgYHK/PnzFUUx/9+7Mjdu3FDuuOMO5ZNPPinz+blz5yp9+/Y1/U0KUQapARC1RmZmJhs2bGDixIlldpxq0KABADNnzuSrr75i5cqV7N27l9atWzNkyJBi7a3mWLlyJT4+PuzatYtJkyYxYcIERo8eze23387evXsZPHgwTz75JNeuXSv1Ond3d3bu3MmCBQuIioril19+Mfn6/fv3Jzg4mLVr1xr2RUdH88knn7B8+XIOHTrEtGnTeOKJJ9iyZYtR109LS+PRRx9l3LhxHDlyhM2bNzNy5EiUMvoOL168mNtuu42nnnqKtLQ00tLSCAwMZPTo0RQWFvLtt98ajr1w4QI//PAD48aNK/O97NmzB1dXV8aPH4+npyc//vgjADdu3ODIkSMWrQH4/fff6dGjR6n9s2bN4s033+S1117j8OHDfP755/j6+hqeN/ffuzyKohAaGkr//v158sknyzymV69e7Nq1i7y8PPPerBBF2ToDEcJSdu7cqQDK2rVryz0mNzdXcXJyUlatWmXYd/PmTUWj0SgLFiww7DOnBqDonVlBQYHi7u6uPPnkk4Z9aWlpCqBs37692HVK3tH17NlTefHFF8t9D+XVACiKooSEhCgdOnRQFEV3N+nm5qb8+eefxY4JDw9XHn30UaOuv2fPHgVQzpw5U+b1Sn5OJR/rTZgwQRk2bJjh8bvvvqu0bNlS0Wq1ZZ53+vTpSq9evQyvfeSRRxRFUZQdO3YogJKcnFzm68wxYsQIZdy4ccX25eTkKC4uLsVqU4oy99+7Ir///ruiUqmU4OBgw/b3338XO2b//v0V/nsIYQpHWyYfQliSYsSI1lOnTpGfn0+fPn0M+5ycnOjVqxdHjhyp0vW7dOli+NvBwYFGjRrRuXNnwz793eOFCxfKfR3oeoOXPMZYiqKgUqkAOHnyJNeuXWPQoEHFjrl58ybdunUz6vrBwcEMGDCAzp07M2TIEAYPHsxDDz1Ew4YNTYrrqaeeomfPnpw/f56mTZsSHx9v6MxYlr179xru8keOHMnIkSPJy8tj7969NG7cuFjnuIULF5KZmcn8+fMBKCgooFGjRly6dAlHx8p/4q5fv46rq2uxfUeOHCEvL48BAwaU+zpz/73L07dv30o7cdarVw/A6FoFISoiCYCoNdq0aYNKpeLo0aMWP7darS6VYOTn5xd77OTkVOyxvpd+0cdAqR/5sl5nbm/+I0eO0KJFCwByc3MB+OGHH2jatGmx41xcXIy6voODA7/88gt//vknP//8M0uXLuWVV15h586dhusYo1u3bgQHB/PJJ58wePBgDh06xA8//FDu8Xv37uXRRx8FdCMLnJyc+Omnn8rsAHjw4EHuueeeYp9B8+bNjSr8AXx8fLh8+XKxffqCtiLm/ntXhb6ZqnHjxhY7p6i7pA+AqDW8vb0ZMmQIMTExXL16tdTzWVlZtGrVCmdnZ7Zt22bYn5+fz19//UXHjh3LPXfjxo1JS0szPM7JySEpKcmyb6CKfvvtNw4cOMCoUaMA6NixIy4uLiQnJ9O6detiW9E76MqoVCr69OnD3Llz2bdvH87Oznz99ddlHuvs7ExhYWGZz40fP574+Hji4uIYOHBguTGcPn2arKwsQ0Hv6OjI/fffz1dffcXevXtLtdcfPHiQ4OBgw+PExMRid+KV6datG4cPHy62r02bNtSrV8/uho4ePHiQgIAAfHx8bB2KqAWkBkDUKjExMfTp04devXoRFRVFly5dKCgo4JdffmHZsmUcOXKECRMmMGPGDLy9vQkKCmLBggVcu3atwmFg/fv3Jz4+nuHDh9OgQQNmz56Ng4NDNb6z4vLy8khPT6ewsJCMjAw2bNhAdHQ09913H2PGjAHAw8OD6dOnM23aNLRaLX379iU7O5tt27bh6enJ2LFjK73Ozp072bhxI4MHD6ZJkybs3LmTf/75hw4dOpR5fPPmzdm5cydnzpyhfv36eHt7o1br7jMee+wxpk+fzkcffcQnn3xS7jX37NmDs7MznTp1MuwbNWqUoUPdK6+8YtivKApHjx4t1sxx+fJlZs2aVel70xsyZAizZs3i8uXLhqYNV1dXXnzxRWbOnImzszN9+vThn3/+4dChQ9UyXLA8v//+O4MHD7bZ9UXtIgmAqFVatmzJ3r17ef3113nhhRdIS0ujcePG9OjRg2XLlgHw5ptvotVqefLJJ7ly5Qq33HILP/30U4Xt2rNmzSIpKYn77rsPLy8v5s2bZ9MagA0bNuDv74+joyMNGzYkODiYJUuWMHbsWEOBCzBv3jwaN25MdHQ0p0+fpkGDBnTv3p2XX37ZqOt4enqydetW3nvvPXJycmjWrBnvvvsuw4YNK/P46dOnM3bsWDp27Mj169dJSkoyTK7j5eXFqFGj+OGHH3jggQfKvebevXvp1KkTzs7Ohn2DBg2isLCQmzdvFmsCSEpKom3btuzevduwb9iwYXTu3Jk1a9bwyy+/cPnyZSIjI4slFEV17tyZ7t278+WXX/LMM88Y9r/22ms4Ojoye/ZsUlNT8ff359lnnzXmY7OKGzdusG7dOjZs2GCzGETtImsBCCGqzYABA/jPf/7DkiVLLHK+b7/9ljVr1vDpp58a9jVr1oytW7fSrFkzQNck8PPPPzNz5sxyz/PDDz8wY8YMDh48WCyBsifLli3j66+/5ueff7Z1KKKWkBoAIYTVXb58mc2bN7N582Y++OADi5334MGDxfpuXL16lezsbEPhr9Vq+eCDD5g9e3aF57n33ns5ceIE58+fN6l/RHVycnJi6dKltg5D1CKSAAghrK5bt25cvnyZt956i3bt2lnsvAcPHuThhx82PD569KghIdBqtUybNo0JEyYQEBBQ6bnsfVGj8ePH2zoEUctIE4AQolZasGAB69ato0uXLgwaNMgwOkIIoSMJgBBCCFEH2WdvFyGEEEJYlSQAQgghRB0kCYAQQghRB0kCIIQQQtRBkgAIIYQQdZAkAEIIIUQdJAmAEEIIUQdJAiCEEELUQZIACCGEEHWQJABCCCFEHSQJgBBCCFEHSQIghBBC1EH/D0Bg9OzUxGtRAAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 512x384 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Polarization fraction as a function of the column density\n",
    "PvN_plot_Step2 = HtoH.PvN(CatP, CatdP, CatN, Nscale=[2e+21,10e+22], Pscale=[0.0,30.0], showfit='true', weighted='false', errorbars='true', Source=FullName)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "c49996f1",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAgAAAAF2CAYAAAAY6yC7AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/P9b71AAAACXBIWXMAAA9hAAAPYQGoP6dpAABZRUlEQVR4nO3de1xU1d4/8M8gAqIMiqBgAuL9giIioqloal5OZabn/Mxjiql1KhSVvICZppagmKbm0afydjpeyh61Tpl15CiWdxEvqJEoBh7vFwZRRGDW74+emRiYGeay5/55v17zUvbsWfu75rLXd6+99toyIYQAERERuRQ3WwdARERE1scEgIiIyAUxASAiInJBTACIiIhcEBMAIiIiF8QEgIiIyAUxASAiInJBTACIiIhcEBMAIiIiF+QyCUBpaSnee+89lJaW2joUyThjnQDnrJcz1glgvRyJM9YJcM56WatOMnuYCnjNmjVYs2YNrly5AgDo0KED5s6diyFDhgAAHj9+jLfffhvbtm1DaWkpBg0ahL///e9o3LixwdsoKiqCr68vFAoF5HK5Japhdc5YJ8A56+WMdQJYL0fijHUCnLNe1qqTXfQANG3aFKmpqcjMzMSJEyfQr18/vPjiizh37hwAYNq0afjXv/6F7du3IyMjA9euXcPw4cNtHDUREZHjcrd1AADwwgsvaPz9wQcfYM2aNThy5AiaNm2KdevWYcuWLejXrx8AYMOGDWjXrh2OHDmC7t272yJkIiIih2YXCUBlFRUV2L59Ox4+fIgePXogMzMTZWVlGDBggHqdtm3bIiQkBIcPH9aZAJSWlmqcPyksLAQAKBQKi8ZvTUVFRRr/OgtnrJcz1glgvRyJM9YJcM56qdopVbul4unpCU9PT+k2JOzEmTNnRN26dUWtWrWEr6+v+O6774QQQmzevFl4eHhUWz86OlrMnDlTZ3nz5s0TAPjggw8++ODDKR7z5s2TtN21mx6ANm3a4NSpU1AoFPjqq68QFxeHjIwMk8tLTk5GYmKi+m+FQoGQkBAUFBQ4zUARaygvL8d///tfNG7cGDdv3sRTTz0Fd3fpvjaq8qUu15jtWTuGmthbPERkXUVFRQgODkZ+fj58fX3VyyU9+ocdnQLw8PBAy5YtAQBRUVE4fvw4VqxYgZEjR+LJkycoLCxE/fr11evfvHkTgYGBOsvT1VUil8uZABihvLwcRUVFkMvlKC4uhlKpRL169eDmJs340crlWysBqLo9a8dQE3uLh4hsw9fX1/mvAtBGqVSitLQUUVFRqF27NtLT09XP5eTkID8/Hz169LBhhK5HCAGFQgGlUmnrUIiIyEx2cXiRnJyMIUOGICQkBA8ePMCWLVuwf/9+/PDDD/D19cWECROQmJgIPz8/yOVyTJ48GT169OAVAERERCayiwTg1q1bGDt2LK5fvw5fX1906tQJP/zwA5599lkAwPLly+Hm5oYRI0ZoTAREREREprGLBGDdunV6n/fy8sLq1auxevVqK0VERETk3Ox2DAARERFZDhMAIiIiF8QEgIiIyAUxASAiInJBTACIiIhcEBMAIiIiF8QEgIiIyAUxASAiInJBTACIiIhcEBMAIiIiF8QEgIiIyAUxASAiInJBTACIiIhcEBMAIiIiF8QEgIiIyAUxASAiInJBTACIiIhcEBMAIiIXtHHjRshkMshkMvWyvn37QiaTYdy4cbYLjKyGCQARkZNRNeTaHrt27QIABAQEICYmBjExMbYNVkLjxo2DTCZD3759DVp/8uTJiIiIgLu7O2QyGQIDA7Wul5mZicGDB0Mul8Pb2xu9evXC3r17JYzcNtxtHQAREVmGh4cHIiMjNZb5+fkBAJ577jk899xztgjLbnz++efw8PCAn58fbt++rXWdM2fOIDY2Fo8ePYK/vz/kcjkOHjyIwYMHY/fu3Rg4cKCVo5YOewCIiJxUUFAQjhw5ovGIjY0FoP0UgDZxcXFo1aoVfHx84OHhgdDQUCQkJKCoqEi9jurIu1mzZti0aROCg4Mhl8sxdepUlJSUYOrUqfD19UVISAjWrl2rUf61a9cwfvx4NGnSBB4eHmjevDkWLlyI8vJy9TqqHo2xY8di3rx5CAoKQoMGDfDKK6/gwYMHAKDeNgBkZGSo67Z//36ddTt79ixu3bqFP/3pTzrXmTNnDh49eoRmzZrh8uXLuHLlCmJiYlBRUYHp06frfe/sHRMA0isrKwujRo1CVlaWrUMhIhv4+uuvcf/+fbRo0QLBwcHIz8/HqlWrMGHChGrrXr9+HW+99RY8PT3x4MEDrFixAlFRUdiwYQN8fX1RUFCA+Ph4XLhwAQBw9+5ddO/eHRs2bEBxcTHatWuHgoICzJ07F6+//nq18rdt24bly5ejTp06KCwsxObNm5GamgoAiIyMhL+/PwDAx8dHfXpDLpfrrFtwcLDeupeXl6u7+gcOHAgfHx+4u7tj6NChAH5PIK5du2bAu2ifmACQXjNmzEBeXh6SkpJsHQoRGem3336rNgbAWBkZGbhz5w5OnTqFS5cu4Z133gEA7Nq1C48fP9ZY98mTJ/jxxx/x66+/qhvX3NxcnD59GufPn4eXlxeUSiUyMjIAAB9//DEKCgrQuHFjXLp0CadPn8ZXX30F4PceitzcXI3yvby8cOHCBeTm5iIqKgoAkJ6eDgDYuXOn+pRGly5d1D0eXbp0MbrOKnfu3EFJSQkAoFGjRurljRs3Vv8/Pz/f5PJtjWMASK+0tDQkJCSos2wichzaxgAYa+/evRg9ejQuXbqk0eCXl5fj9u3bGkfRDRo0QM+ePQEAISEhKCgoQHh4OJo1awbg90Y0Pz8fN2/eBAAcO3YMAHDz5k2NBhYAhBA4evQoWrZsqV7Wr18/PPXUUwCAtm3bIjMzU12WNQkhrL5NS2ACQHpFRkZi69atCAwMxH//+19bh0NERlCNATDV5s2b1ee5g4KCEBwcjDt37uDy5csAgIqKCo31K3e3u7u7V1um6oGo2oD6+Pigffv21bbv7e2t8Xf9+vWrlW/Jxtjf3x916tRBSUkJbt26pV5e+f8hISEW276l8RQAuRSOaSAynCp58PHxQV5eHo4ePSrpqPfo6GgAvzfm27ZtU3fb//vf/8Zbb72Fl156yajyVAnDw4cPJYnP3d0d/fv3BwD8+OOPePDgAcrLy/HNN98AADp27IgmTZpIsi1bYAJALkU1pmHGjBm2DoXI7nXq1AkA8ODBAzRv3hzNmzfHl19+KVn58fHxeOqpp3D//n20adMGnTt3RosWLdCwYUPExcUZXV7btm0BACdOnEDHjh3RvXt39Tl8bfr27YuWLVtix44dAH4/59+yZUu0bNkSR48eBQC8//77qFOnDq5cuYLmzZujWbNmOHr0KGrVqoUlS5aYUGv7wQSAbMraR+RpaWkICwtDWlqaVbZH5MgmTJiAxMRE+Pv748GDB+jbty8WLFggWfkBAQE4cuQIXn31VTRs2BDnzp1DSUkJevfujeXLlxtd3vjx4zFixAj4+voiOzsbR48erXaaorIrV67g0qVL6ksJKyoqcOnSJVy6dEmdOERERCAjIwPPPvssHj9+jLt37+Lpp5/G7t27MXjwYNMqbidkwllGM9SgqKgIvr6+UCgUei8LIU3l5eW4evWqxhiA0NBQ9fk3c/Xt2xe//PIL2rZtq/d6Xamo6tO0aVN1HbQtsyV7i4eIrMta7RV7AMimeERORGQbPLwgm1JdZdC0aVNbh2ISpVKJoqIiyOVyuLkxnyYix8E9FpEZlEolCgsLoVQqbR0KEZFRmACQy+OlgUTkipgAkEGysrIwevRoZGdn2zoUyfHSQCJyRUwAyCBJSUnIy8tDSkqKrUORHAciEpErYgJABklNTUVYWBiSk5NtHYrkVAMRzZ0znYjIkfAqADJIZGQkNm/ebOswiIhIIuwBICIickFMAIiIiFyQXSQAKSkpiI6Oho+PDxo1aoRhw4YhJydHY52+fftCJpNpPN544w0bRUxEROTY7CIByMjIQHx8vPo2kGVlZRg4cGC1Wzq+9tpruH79uvrh6HdicnacJIeIyH7ZxSDAPXv2aPy9ceNGNGrUCJmZmYiNjVUv9/b2RmBgoLXDIxOpEoB69epxmlwiIjtjl3tlhUIBAPDz89NYvnnzZvj7+yM8PBzJycl49OiRzjJKS0tRVFSk8SAiInIUVduw0tJSScu3uwRAqVRi6tSp6NmzJ8LDw9XL//rXv+Kf//wn9u3bh+TkZHz++ed45ZVXdJaTkpICX19f9SM4ONga4Tu17OxsjB49mlPmEhFZQXBwsEY7JvVEbHZxCqCy+Ph4ZGdn4+eff9ZY/vrrr6v/37FjRwQFBaF///64dOkSWrRoUa2c5ORkJCYmqv8uKipiEmCmlJQU9ZS5+/fvt3U4REROraCgAHK5XP23p6enpOXbVQ/ApEmT8O2332Lfvn013h42JiYGAJCbm6v1eU9PT8jlco0HmSc5OZlT5hIRWUnVNkzqBMAuegCEEJg8eTJ27tyJ/fv3IywsrMbXnDp1CgAQFBRk4ehIJTw8HJs3b0ZoaKitQyEiIjPZRQ9AfHw8/vnPf2LLli3w8fHBjRs3cOPGDZSUlAAALl26hIULFyIzMxNXrlzBN998g7FjxyI2NhadOnWycfRE0uLtiYnIGuwiAVizZg0UCgX69u2LoKAg9eOLL74AAHh4eGDv3r0YOHAg2rZti7fffhsjRozAv/71LxtHTiQ93p6YiKzBbk4B6BMcHIyMjAwrRUNkW2lpaUhISOBYCyKyKLtIAIjoD6rbE9c0EJaIyBx2cQqAiIiIrIsJABERkQtiAkBEROSCzBoDUFZWhhs3buDRo0cICAioNnc/ERER2SejewAePHiANWvWoE+fPpDL5WjWrBnatWuHgIAAhIaG4rXXXsPx48ctESsRERFJxKgEYNmyZWjWrBk2bNiAAQMGYNeuXTh16hR+/fVXHD58GPPmzUN5eTkGDhyIwYMH4+LFi5aKm4iIiMxg1CmA48eP48CBA+jQoYPW57t164bx48dj7dq12LBhA3766Se0atVKkkCJiIhIOkYlAFu3bjVoPU9PT7zxxhsmBURERESWJ8lEQHfu3MHRo0dRUVGB6Oho3qCHiIjIzpmdAPzv//4vJkyYgNatW6OsrAw5OTlYvXo1Xn31VSniI9JJqVSiqKgIcrkcbm68opWIyBhG7zWLi4s1/p4/fz6OHTuGY8eOISsrC9u3b8c777wjWYBEuiiVShQWFkKpVNo6FCIih2N0AhAVFYWvv/5a/be7uztu3bql/vvmzZvw8PCQJjoiIiKyCKNPAfzwww+Ij4/Hxo0bsXr1aqxYsQIjR45ERUUFysvL4ebmho0bN1ogVCIiIpKK0QlAs2bN8N1332Hr1q3o06cPEhISkJubi9zcXFRUVKBt27bw8vKyRKxEREQkEZNHTo0aNQrHjx/H6dOn0bdvXyiVSnTu3JmNPxERkQMw6SqA3bt348KFC4iIiMBnn32GjIwMjB49GkOGDMGCBQtQp04dqeMkIiIiCRndA/D222/j1VdfxfHjx/G3v/0NCxcuRJ8+fXDy5El4eXkhMjIS33//vSViJSIiIokYnQBs3LgRu3fvxrZt23D8+HF8/vnnAAAPDw8sXLgQO3bswKJFiyQPlIiIiKRjdAJQt25d5OXlAQAKCgqqnfNv3749fvrpJ2miIyIiIoswOgFISUnB2LFj0aRJE/Tp0wcLFy60RFxERERkQUYPAhw9ejQGDx6My5cvo1WrVqhfv74FwiIiIiJLMqoHID8/HwDQsGFDREdH19j4//e//zU5MCIiIrIcoxKA6Oho/O1vf8Px48d1rqNQKPDpp58iPDwc//u//2t2gERERCQ9o04BnD9/Hh988AGeffZZeHl5ISoqCk2aNIGXlxfu37+P8+fP49y5c+jSpQuWLFmCP/3pT5aKm4iIiMxgVA9Aw4YNsWzZMly/fh0ff/wxWrVqhTt37uDixYsAfh8fkJmZicOHD7PxJ4fh5uaG+vXr85bCRORSTJoJsE6dOvjzn/+MP//5z1LHQ2R1qgSAiMiV8JCHiIjIBTEBIL3YPU5E5Jy4Vye9mAAQETkn7tWJiIhckEmDAH/77TecOXMGjRs3Rrdu3aSOiYiIiCzM6ARg69atGDduHMrKyiCTydS3/w0ICLBEfERERGQBRp8CmD9/Pv7617/il19+wY8//ggASEpKkjwwIiIishyjewAuX76MPXv2oFmzZmjdujX++c9/IioqCuvWrbNEfERERGQBRvcAlJeXw9vbW/1327ZtoVQqcePGDUkDIyIiIssx6SqATZs24dChQyguLgYAuLu749GjR5IGRkRERJZjdALQu3dvvP/+++jVqxfq16+PVq1a4fHjx1i3bh327duHBw8eWCJOsjE3Nzf4+vpCJpPZOhQiIpKA0WMAMjIyAAAXL15EZmYmTp48iZMnT2LNmjVISUmBm5sbWrVqhQsXLkgeLNmOakIgVa+Po+LERkREvzNpHgAAaNWqFVq1aoWXX35ZvSwvLw8nTpxAVlaWJMERSY03/iEi+p2kh0FhYWH4y1/+gkWLFhn1upSUFERHR8PHxweNGjXCsGHDkJOTo7HO48ePER8fj4YNG6JevXoYMWIEbt68KWX4RARAqVSisLAQSqXS1qEQkQXZRT9oRkYG4uPjceTIEfz73/9GWVkZBg4ciIcPH6rXmTZtGv71r39h+/btyMjIwLVr1zB8+HAbRk3knJgAELkGk08BSGnPnj0af2/cuBGNGjVCZmYmYmNjoVAosG7dOmzZsgX9+vUDAGzYsAHt2rXDkSNH0L1792pllpaWorS0VP13UVGRZStB1WRlZSEhIQErV65EdHS0rcMhPZRKJYqKiiCXy20dChH9n6rtlqenJzw9PSUr3y56AKpSKBQAAD8/PwBAZmYmysrKMGDAAPU6bdu2RUhICA4fPqy1jJSUFPj6+qofwcHBlg+cNMyYMQN5eXmYMWOGRcrPysrCqFGjbDrmxB5ikAKP+onsT3BwsEY7lpKSImn5dpcAKJVKTJ06FT179kR4eDgA4MaNG/Dw8Kg2eKtx48Y6JyBKTk6GQqFQPwoKCiwdutMzdgR9WloawsLCkJaWZpF4LJlgGNogWjrJISLXVVBQoNGOJScnS1q+WacA0tPTkZ6ejlu3blXbUa5fv96kMuPj45GdnY2ff/7ZnNAk7yoh40fQR0ZGYuvWrWjatKlF4klLS0NCQoJFEgxVAlCvXj29CY8lYyAi1yaXyy16Ws7kBGD+/PlYsGABunbtiqCgIEkmiJk0aRK+/fZbHDhwQKPRCAwMxJMnT1BYWKjRAN28eROBgYFmb5cck6UTDEeJgYjIFCYnAGvXrsXGjRsxZswYs4MQQmDy5MnYuXMn9u/fj7CwMI3no6KiULt2baSnp2PEiBEAgJycHOTn56NHjx5mb5+IiMjVmJwAPHnyBE8//bQkQcTHx2PLli34+uuv4ePjoz6v7+vrizp16sDX1xcTJkxAYmIi/Pz8IJfLMXnyZPTo0UPrFQBERESkn8mDACdOnIgtW7ZIEsSaNWugUCjQt29fBAUFqR9ffPGFep3ly5fj+eefx4gRIxAbG4vAwEDs2LFDku0TERG5GpN7AB4/foxPPvkEe/fuRadOnVC7dm2N55ctW2ZwWUKIGtfx8vLC6tWrsXr1aqNjJSIiIk0mJwBnzpxB586dAQDZ2dkaz/GOcURERPbN5ARg3759UsZBREREVmR3EwERERGR5Zk1EVBhYSHWrVuHCxcuAADat2+PCRMmwNfXV5LgyPkZO7sgERFJw+S97okTJ9CiRQssX74c9+7dw71797B8+XK0aNECJ0+elDJGcmJMAKrje0JE1mByD8C0adMwdOhQfPrpp3B3/72Y8vJyTJw4EVOnTsWBAwckC5LIlRg75TIRkSlMTgBOnDih0fgDgLu7O2bOnImuXbtKEhwRERFZhsl9jHK5HPn5+dWWFxQUwMfHx6ygiIiIyLJMTgBGjhyJCRMm4IsvvkBBQQEKCgqwbds2TJw4EaNGjZIyRiIiIpKYyacAli5dCplMhrFjx6K8vBwAULt2bbz55ptITU2VLEAiIiKSnskJgIeHB1asWIGUlBRcunQJANCiRQt4e3tLFhwRERFZhlnzAACAt7c3OnbsKEUsREREZCVGJQCJiYlYuHAh6tati8TERL3rGnMzICIiIrIuoxKArKwslJWVqf+vC28GREREZN+MSgAq3wBo06ZNaNq0abXZyoQQKCgokCY6IitRKpUoKiqCXC7nDHxE5BJM3tOFhYXhzp071Zbfu3cPYWFhZgVFZG1KpRKFhYVQKpW2DoWIyCpMTgCEEFqXFxcXw8vLy+SAiIiIyPKMvgpANfhPJpNh7ty5Gpf9VVRU4OjRo+jcubNkARIREZH0jE4AVIP/hBA4e/YsPDw81M95eHggIiIC06dPly5CIiIikpzRCYBqIOCrr76KFStWQC6XSx4UERERWZbJEwFt2LBByjiIrIqj/onI1Zm850tJScH69eurLV+/fj0WL15sVlBElsZR/0Tk6kxOAP7nf/4Hbdu2rba8Q4cOWLt2rVlBEZHlMQkicm0mJwA3btxAUFBQteUBAQG4fv26WUERkeUxASBybSYnAMHBwTh48GC15QcPHkSTJk3MCoqIiIgsy+RBgK+99hqmTp2KsrIy9OvXDwCQnp6OmTNn4u2335YsQCIiIpKeyQnAjBkzcPfuXbz11lt48uQJhBCoU6cOZs2ahaSkJCljJCIiIomZnADIZDIsXrwY7777Li5cuIA6deqgVatW8PT0lDI+Ip3c3NxQv359XsZHRGQCkxMAlfz8fNy9exdPnjzB5cuX1cuHDh1qbtFEeqkSAEvIyspCQkICVq5ciejoaItsg4jIlkxOAC5fvoyXXnoJZ8+ehUwmU98cSCaTAfj9vgBEjmrGjBnIy8vDjBkzsH//fluHQ0QkOZP7TqdMmYKwsDDcunUL3t7eyM7OxoEDB9C1a1fuME3Ey7LsR1paGsLCwpCWlmbrUIiILMLkHoDDhw/jP//5D/z9/eHm5oZatWqhV69eSElJQUJCgvqmQWQ4VQJQr149nte2scjISGzduhVNmza1dShERBZhcitTUVEBHx8fAIC/vz+uXbsGAAgNDUVOTo400REREZFFmNwDEB4ejtOnTyMsLAwxMTFYsmQJPDw88Mknn6B58+ZSxkhEREQSM7kHYM6cOepz1QsWLEBeXh569+6N3bt3Y+XKlZIFSESWl5WVhVGjRvHUHZELMTkBGDRoEIYPHw4AaNmyJX755RfcuXMHt27dUs8MSGQLHExpvMpXPTAZIHINJiUAZWVl6N+/Py5evKix3M/PT30ZIJGtMAEwjFKpRHFxMZRKpcZVD5WTASJyXiYlALVr18aZM2ekjoWIrKhyAqC66iEyMpKXQBK5CJNPAbzyyitYt26dlLEQkR2onAwQkfMy+SqA8vJyrF+/Hnv37kVUVBTq1q2r8fyyZcsMLuvAgQNIS0tDZmYmrl+/jp07d2LYsGHq58eNG4dNmzZpvGbQoEHYs2ePqeETERG5NJMTgOzsbHTp0gUA8Ouvv2o8Z+w4gIcPHyIiIgLjx49XDyysavDgwdiwYYP6b950iIiIyHQmJwD79u2TLIghQ4ZgyJAhetfx9PREYGCgZNskcnRKpRJFRUWQy+WcOZKIjGb0XuPy5cvqG/9Y0/79+9GoUSO0adMGb775Ju7evat3/dLSUhQVFWk8iJwJr3Ygcm5V27DS0lJJyzc6AWjVqhVu376t/nvkyJG4efOmpEFVNXjwYPzjH/9Aeno6Fi9ejIyMDAwZMkTvHQdTUlLg6+urfgQHB1s0RiIiIikFBwdrtGMpKSmSlm90AlD16H/37t14+PChZAFp8/LLL2Po0KHo2LEjhg0bhm+//RbHjx/Xe9fB5ORkKBQK9aOgoMCiMZJj4WQ3RGTvCgoKNNqx5ORkSct3yBOHzZs3h7+/P3Jzc3Wu4+npCblcrvEg63Jzc0P9+vXt8vy0o052w25/ItdRtQ2TevC70XtmmUxWbZS/tWf/u3r1Ku7evYugoCCrbpeMY88JgKNOdsMEgIikYvRVAEIIjBs3Tp2JPH78GG+88Ua1eQB27NhhcJnFxcUaR/N5eXk4deoU/Pz84Ofnh/nz52PEiBEIDAzEpUuXMHPmTLRs2RKDBg0yNnwiAH9MdtO0aVNbh0JEZBNGJwBxcXEaf7/yyitmB3HixAk888wz6r8TExPV21qzZg3OnDmDTZs2obCwEE2aNMHAgQOxcOFCzgVARERkIqMTgMqT8Uilb9++ei8t/OGHHyTfJpEhVF3ujnKtPecGICJDcQ9BpIejnXN3tHiJyHaYADgQ7twth5cFEpGrYQLgQJgAWI6jXhZIRGQqJgBEcNzLAomITGXyzYCInAkvCyQiV2NWApCeno709HTcunWrWrf0+vXrzQqMiIiILMfkBGD+/PlYsGABunbtiqCgIKvPBkhERESmMzkBWLt2LTZu3IgxY8ZIGQ8RERFZgcmDAJ88eYKnn35ayliIHI493++AiEgfk/daEydOxJYtW6SMhcjhMAEgIkdl8imAx48f45NPPsHevXvRqVMn1K5dW+P5ZcuWmR0cERERWYbJCcCZM2fQuXNnAEB2drbGcxwQSGT/srKyMGnSJHz88ceIjo62dThEZGUmJwD79u2TMg4isrKkpCTk5+cjKSkJ6enptg6HiKyMJy6JXFRqaipCQkKQmppq61CIyAbMmgiosLAQ69atw4ULFwAA7du3x4QJE+Dr6ytJcK4mKysLCQkJWLlyJbtkyeIiIyOxevVq9ak8InItJvcAnDhxAi1atMDy5ctx79493Lt3D8uXL0eLFi1w8uRJKWN0GbwhjW1ou8kSR/cTkbMzee82bdo0DB06FFeuXMGOHTuwY8cO5OXl4fnnn8fUqVMlDNF18IY00jD21r5KpRIKhQJCCPUye00AlEoliouLeUdIIjKbWT0As2bNgrv7H2cR3N3dMXPmTJw4cUKS4FyN6oY0kZGRtg7FoTlzT0pmZibi4uKQmZlp61CIyMGZnADI5XLk5+dXW15QUAAfHx+zgiIyhzP3pFQeuW+urKwsxMfHG9xTQkTOxeQEYOTIkZgwYQK++OILFBQUoKCgANu2bcPEiRMxatQoKWMkMooz96RIOXJfymSCiByPyVcBLF26FDKZDGPHjkV5eTkAoHbt2njzzTd5WRHZvcrn+B3pfLqUI/dTU1MxadIk/l6JXJTJCYCHhwdWrFiBlJQUXLp0CQDQokULeHt7SxYckaWoEgAADpUASImXARK5NrPmAQAAb29vdOzYUYpYiIiIyEqMSgASExOxcOFC1K1bF4mJiXrX5c2AiIiI7JdRCUBWVhbKysrU/9eFNwMiIiKyb0YlAJVvALRp0yY0bdq02kQpQggUFBRIEx2RC1AqlSgqKoJcLjd74iFOJ01EhjJ5bxMWFoY7d+5UW37v3j2EhYWZFRSRo9I2rXBN65aXlxv8mpo48yRIRCQtkxOAytOmVlZcXAwvLy+TAyJyZKYkAKY2/NqmPHbmSZCISFpGXwWgGvwnk8kwd+5cjcv+KioqcPToUV5WRE7DnrvUKx/t79+/H8AfkyA1bdrUtsERkd0zOgFQHW0IIXD27Fl4eHion/Pw8EBERASmT58uXYRENqStkbUXaWlpSEhIkORo31EnRiIi0xmdAKgGAr766qtYsWIF5HK55EER2QspG1mpSXm0z4mRiFyPyRMBbdiwAQBw/vx55Ofn48mTJxrPDx061LzIiOwAu9SJyFmZnADk5eVh2LBhOHv2LGQymXpQoGoOgIqKCmkiJCIiIsmZfBVAQkICwsLCcOvWLXh7e+PcuXM4cOAAunbtanfnSolclblXGhCR8zI5ATh8+DAWLFgAf39/uLm5wc3NDb169UJKSgoSEhKkjJHIaVi7QTZle5UHBBKR8zL5F15RUQEfHx8AgL+/P65duwYACA0NRU5OjjTRkQZt132TY3GEI3ImAESuweRfeHh4OE6fPg0AiImJwZIlS3Dw4EEsWLAAzZs3lyxA+gNneSMiIqmYnADMmTNHfRSzYMEC5OXloXfv3ti9ezdWrlwpWYD0B87yRkREUjH5KoBBgwap/9+yZUv88ssvuHfvHho0aMC7AVoIL0kzj66b7mRlZWHy5MlITk5GaGioDSMkIrIeSU/y+fn5mdT4HzhwAC+88AKaNGkCmUyGXbt2aTwvhMDcuXMRFBSEOnXqYMCAAbh48aJEUZOr0HX+XXVqJSUlxUaRERFZn8kJQL9+/TB//vxqy+/fv49+/foZVdbDhw8RERGB1atXa31+yZIlWLlyJdauXYujR4+ibt26GDRoEB4/fmxS7ESVqU6tLFq0yO4Hvp0+fRoJCQnq8TdERKYyeW+3f/9+fPzxxxg2bBgePnyoXv7kyRNkZGQYVdaQIUPw/vvv46WXXqr2nBACH330EebMmYMXX3wRnTp1wj/+8Q9cu3atWk8BkSkiIyOxbds2xMbG2iQBMObqjlmzZiE/Px+zZs2yQmRE5MzM2tvt3bsXN27cQPfu3XHlyhWJQtKUl5eHGzduYMCAAeplvr6+iImJweHDh3W+rrS0FEVFRRoPIlNY+tK9pKQkg6/u4EBQItdRtQ0rLS2VtHyzEoCgoCBkZGSgY8eOiI6OtsgMgDdu3AAANG7cWGN548aN1c9pk5KSAl9fX/UjODhY8tjINQghoFAoDEoAjDmaV607btw4NGvWDHPnzq1xG6qBoJGRkQbHT0SOKTg4WKMdk3qckskJgGqwn6enJ7Zs2YIpU6Zg8ODB+Pvf/y5ZcOZITk6GQqFQPwoKCmwdErkAY+ZqUK27ceNGbN68GaGhoeoEwBEmDCIiyyooKNBox5KTkyUt3+QEQHXzH5U5c+Zg8+bN+PDDD80OqrLAwEAAwM2bNzWW37x5U/2cNp6enpDL5RoPIkszpIte1bgvXrwYYWFhSE1N1bmOJRIAJhdEjqFqG+bp6Slp+SYnAHl5efD399dYNmLECBw5cgTr1683OzCVsLAwBAYGIj09Xb2sqKgIR48eRY8ePSTbDpEUDOmiV51SiIiIsEl3PhMAIgLMmAhI14Qp4eHhCA8PN6qs4uJi5Obmqv/Oy8vDqVOn4Ofnh5CQEEydOhXvv/8+WrVqhbCwMLz77rto0qQJhg0bZmr4RERELs2oBCAxMRELFy5E3bp1kZiYqHfdZcuWGVzuiRMn8Mwzz2hsBwDi4uKwceNGzJw5Ew8fPsTrr7+OwsJC9OrVC3v27IGXl5cx4RNRJUqlEiUlJewJIHJRRiUAWVlZKCsrU/9fF2NnA+zbt2+1MQVVy1uwYAEWLFhgVLlEpBsTACLXZlQCsG/fPgBAWVkZ3NzcsHbtWrRq1coigbki3oZVGnwfiYhqZtIYgNq1a+PMmTNSx+LyVA0XmYfvIxFRzUw+RHrllVewbt06KWMhcmq2GH2vVCpRXFysdZvu7u5o2LAh3N1NHgtMRA7M5F9+eXk51q9fj7179yIqKgp169bVeN6YQYBEzkx1G2Jvb28oFAq9413MVfX0R00JQOPGjZkAELkok3/52dnZ6NKlCwDg119/1XjOlFsCEzkr1ZG/Na5a4ekPIjKUyQmAakAgkStQNeJ+fn46Bxdygh0iciQcJk1kAENuCKRUKi3exU9EJBWzT/6dP38e+fn5ePLkicbyoUOHmls0OSnVOXG5XM5L9YiIbMTkBODy5ct46aWXcPbsWchkMvVRj+r8f0VFhTQRktNRdZXXq1fPqglAVlYWEhISsHLlSkRHRxv0Gjc3N/j6+uL+/fsWjs4ysrKyMGnSJHz88ce8dwYRaTB57ztlyhSEhYXh1q1b8Pb2xrlz53DgwAF07doV+/fvlzBEImkYc6telcqD6qqOpjflnH92djbGjBmD7OxsjW34+vpqDJ7NysrCqFGj9M64aYikpCTk5+cjKSnJrHKIyPmYnAAcPnwYCxYsgL+/P9zc3ODm5oZevXohJSUFCQkJUsZIJAlDbtWrixBCkgQgJSUFV65cQUpKinqZtpkLTUlWtElNTUVISIjWWw4TkWszOQGoqKiAj48PAMDf3x/Xrl0D8PtdAnNycqSJjkhChtyq15KUSiWmTp2KgIAAPHz4UO/RvTnJSmWRkZFYvXq1zepMRPbL5AQgPDwcp0+fBgDExMRgyZIlOHjwIBYsWIDmzZtLFiCRpZlz7wBVL0B5eTkKCwshhICvr6/Wsk6ePIklS5YAAG7dulWtW75yj4I1khWpTjMQkWMyOQGYM2eOuutzwYIFyMvLQ+/evbF7926sXLlSsgCJLM2cBEB1eWB5eTkUCgVkMpnOsmbPno2CggIAQFhYGJYsWaKxriGXGlZlztwDSUlJkpxmICLHZPJVAIMGDVL/v2XLlvjll19w7949NGjQgDMBEmmxaNEiJCQk4J133kHnzp0RGhqqnobXlAbczc0NcrkcRUVFJiUwqampSExMNPs0AxE5JkmvwfLz82PjT6RD586dsWrVKrRv397g1+g7wje250JVVmZmJkaNGgUANh0TYQucrZHoD0b1ACQmJhq8Lm8GRM5CdS39G2+8gaioKKtuu+qcCVKMV5g5cyby8vKQlJSETZs2WSBq+2WrOSiI7JFRCYChg4XYC0C2YolZBmfMmIErV65g1apV6iNnc1W9t0BWVhYmT56M5ORkhIaG6nydFDf7UXX989JAItdmVALAGwCRvbPEEV5aWpq6B8AUqhkI09LS1MmxasCf6mhedd1/SkoKnnvuOUni1kV1hUFgYCBu3Lhh0W0Rkf0yaw9ZWFiIDz/8EBMnTsTEiROxfPlyKBQKqWIjsguRkZFISkrCqlWrNHrBsrKyMHr0aPWsftnZ2Rg9enS1njJV4z579myd21Bd95+cnGyZShARVWFyAnDixAm0aNECy5cvx71793Dv3j0sW7YMLVq0wMmTJ6WMkcjmUlNTUVBQoHHtfuWjduD3Wf60XVanatwXLVqks/zIyEhs3rwZ4eHh1Z4zZ+BaVlYW4uPjea0/EVVjcgIwbdo0DB06FFeuXMGOHTuwY8cO5OXl4fnnn8fUqVMlDJFUzBkARuZJSkpCcHCwxnnzqkftycnJWmfvU3W5d+7c2eDtKZVK9dTDpswPUDlu3guAiLQxqwdg1qxZ6uuYAcDd3R0zZ87EiRMnJAmONDEBsJ0OHTpg1apVGpfMVT1qDw8Px+bNm6tdVle5MTeUKa/RhvcCICJdTG5J5HI58vPzqy0vKChQ3yOAyFmcO3cOkydPNqkrXarG3BSuei8AXu9PVDOTE4CRI0diwoQJ+OKLL1BQUICCggJs27YNEydOlOxSKSJ7oW0MgClkMhmEEBJF5Xis1TAzASCqmckJwNKlSzF8+HCMHTsWzZo1Q7NmzTBu3Dj8+c9/xuLFi6WMkeyQuQPTHOUmNKp6zpo1q9oYAED3yH9tTp8+jUmTJuHChQuWCtfusWEmsh8mJwAeHh5YsWIF7t+/j1OnTuHUqVO4d+8eli9fDk9PTyljJDtkzo5cqnvdW4NSqYRCoUC7du2qjQEA/hj5b0jPwJw5c1BQUGBQguzm5gZvb28UFRXpfI+lbkzZOBO5FqMTAKVSicWLF6Nnz56Ijo7GggUL0LJlS3Ts2BHe3t6WiJGcjFT3urcHqpH/48aN05gTQJv3338fwcHBmDVrVo3lurm5oW7duiguLtZ5ykCpVOLgwYP461//KklvChMAItdidALwwQcfYPbs2ahXrx6eeuoprFixAvHx8ZaIjZyUNe51by2qkf8bN27UmBNAm4iICJNuBqRvAKGqB2LatGnV1qnaoFc99cKrSohcm9G//H/84x/4+9//jh9++AG7du3Cv/71L2zevJlHDeTSUlNTDZ7JTyaTwdfXFzKZrMajbiEEHj58qLMXIDk5Gc2aNcOMGTO0JgBFRUXq16pOvahOVzhjAsBeDCLDGf3Lz8/Px5/+9Cf13wMGDIBMJsO1a9ckDYwIkHaHbk6DJ5PJ9N5fQN9MftriuHz5MkaPHo3MzEx1/YwZUKgSHh6Ozz//3KDtqk69OPOcADV9XxxpACqRpRm9JywvL4eXl5fGstq1a6OsrEyyoIhUpEgAVGUAMDkBcHNzk/QGQ0lJSbh8+TKmT5+urpuuqYSl4oinXqQ+onekAahElmb03kwIgXHjxmH48OHqx+PHj/HGG29oLCOyF5bqFnZzc1N35RtLdcpg+vTp6i56XVMJ//rrrxgzZozTH7Vq+5yk/uycaQAqkbmMuh0wAMTFxVVb9sorr0gSDJG9UiqVKCkp0WiIVKcUiouLjS5PdcqgoqJCvUw1oDA0NFRj3Y8//hjXrl3DjBkzsH//fpPrYO+kvJVzVlYWEhMTsXLlSkRHR6uXq3pBmjZtam64RA7P6ARgw4YNloiDSBJZWVlISEiotuM3p7zJkydjxowZCAwM1HskKpPJIJfLqzVeWVlZmDRpEpYtW6Zx7wxDnD9/HiUlJQgICHDpo1ZVcqDt/dUmKSlJ3dXvzEkTkTmcZ/gvEaQ/x6sqTzV5j74uaV2DDJOSklBQUIA5c+boTBJ0Wbx4Me7evYu6devqPHd/7tw5nYMHaxq86Ciqvu81nRpQnWJx5aSJqCaOvVcgqkLqc7yq8lST95hyTnrp0qVo0aIFUlJSDL4SISsrC8OGDcPdu3fRsGFDrTMNZmVlYfTo0Xj33Xd1Jj3aBi9WjsFRL5urKW5HHPBIZG1MAMipGLrjN7ThU52rN2bynqrbCQsLw+bNm9G5c2eDXzdjxgzk5OSgoKAAderUQYcOHarFreqdAGBU0qMrAbDlXQstwVGTGyJrYQJALsmad6VTKBQQQuDUqVOIj483aDR/Wloa2rRpg9atW2PSpEka5aniVvVOLFy4EJs3bzb7aNceEwBjr9uvvL4QAgqFwq7qQ2RPmAAQWcns2bORn5+v88ZBlRvgyMhI7Nq1C9u3b0fr1q21rq/qnVD1DlSlOkWg7/4EUrJEUlV19kKp1ydyZQ6RALz33nuQyWQaj7Zt29o6LJfGGdV+Z8zsgosWLUJISIjOmfiEEJIegasaQ333J6gsKyvL4B4KbQxJAIz93hg7e6ErzHZIJBWHSAAAoEOHDrh+/br68fPPP9s6JJdmrRnV7D3RMCYB6NKlCzZt2oSoqCiDyz9//jwSEhJw7tw5o2NTNYaLFi0y+NI5fT0UUtD3van6WVcePwFA4zld34vKY0BMmVrZVBxvQI7IYRIAd3d3BAYGqh/+/v561y8tLUVRUZHGg6RjrRnVpEg0zE0iZDIZfHx8TJrxr7KqI/Kzs7MxZswYvV30ixcvxtWrV006oo2MjMS2bdsQGxtrUAKQmpqqt4dCCvq+N1U/68rjJypf169tXW0sPbVyZUwAyBKqtmGlpaWSlu8wCcDFixfRpEkTNG/eHKNHj0Z+fr7e9VNSUuDr66t+BAcHWylS12Cty6ykSDTMTSLc3NyQn5+PqVOnYvv27ZL1SKSkpODixYvqaX7lcnm1JGPWrFlo2rSpWUfl+u6FoC85skSjpu97o++zrnpdvyHfC11TKxM5iuDgYI12zNDTeQYTDmD37t3iyy+/FKdPnxZ79uwRPXr0ECEhIaKoqEjnax4/fiwUCoX6UVBQIAAIhUJhxcidV1lZmcjLyxNlZWUWfa2x29G2/rFjx0T37t3FsWPHjC5XtV5sbKxo3Lix8PX1FY0bNxZ9+vQxKM6ysjKRm5srcnNzRUlJiXqdiooKkZ6eLnx8fET9+vVFbGysKCkpEYcPHxYlJSXq150/f14cOHBA5OTkiLKyMo3tlJWViUuXLombN2+Ku3fvioqKCq0xVd5uVX369FHXJzY2Vvj7+4vY2FjJ3ntj1qn63LFjx0RMTIzYtWtXtTroe79VdVa975XXMed7q4+lyiXXpFAoBABRUFCg0Y49fvxY0u04RA/AkCFD8Je//AWdOnXCoEGDsHv3bhQWFuLLL7/U+RpPT0/I5XKNB7kmc3orVOf4lyxZgrCwMCxfvlzrUaUxYwGUSiWKiorQq1cvbNq0qcZBa/pmD3Rzc4Ofnx/8/PxMmu2v8pG0JU8BmNKbYOwgxsrMuVGTKex9rAo5pqptmKenp6TlO0QCUFX9+vXRunVr5Obm2joUMoExjaU5VI1OeXl5tdH1hu6wVbFGRUVh69atGDVqlNZkwtgEQNUYdujQAatWrUJkZKTOUfgymUxdtqFxGzPRkao+kZGRWL16tdZEydzTAaa83thBjJVZ6zumwtsMkyNyyASguLgYly5dQlBQkK1DIRPYOgFQKpWYNm0aLl++bFc7bENG4Rva0FQeQCcFWwxyM3YQY1XavmeW+u7xNsPkiBwiAZg+fToyMjJw5coVHDp0CC+99BJq1aqFUaNG2To0ckCqaXTDwsKwePFiizdsMpkMvr6+NTY6VbvgtV19wIbGcMYmAOYkObz3ADkih0gArl69ilGjRqFNmzb4f//v/6Fhw4Y4cuQIAgICbB0aOajw8HBs3rwZERERFk8ADD3qrNoFr+11lRsaa/SkGNMoSnEevPL29NVP13PmvCe8lI9cjUMkANu2bcO1a9dQWlqKq1evYtu2bWjRooWtw3J43OGZztD3TlujqGvOfW137tNHX2NnyjTA2rYv/m8+/czMzBobdynOg4tK8/dbOwGQCn9X5CgcIgEgy3D2HZWq8T116pTkZRv63mmb3ObevXsoKiqCUqnEuXPnMHnyZGRlZRmdABiyXWNG0OvbftWJeLSR4vSEvuTIEl33lZnTg+EMt1gm18MEgJxK5R2xqhGcPXu2Rco3RNVGUQiBBw8eqAfnpaamoqCgQPLpd1XbTU5O1rueofWpOhGPNlKcB8/Ozsarr76KzMxMg+OUqsE1pwfD2j0PptaZyQlVxgSAnErlHXHly8h0Mfaoz9gdfdVz9lWvTU9KSkJwcHC1a+/NbVAiIyOxZcsW9OzZE6dPn9ZZR0NPIxjbuJvS0Li5uWHp0qVGJ0SVP0NzGjhHGmDJBICkwASALM7QnY7UOydVo9W5c2eN5ZUbNkO6tqWKUVtjW3keAG3ruru7m3xaQFXGrFmzTDqyNeY0QtX3perfhiRaqgTA2MmIKh+5cyQ/keGYAJDFmZMAWKJrtXLDZkjXduU4AGitizENj+q8f9WxCbouWzN3XICpR7aGnkYAap53QPWe13Rkr28yopridIQj98rMSVY48yBJgQkA2TUpEoAzZ84gISEBp0+fBqDZsFU96tO1U9YXR3l5Oa5cuYI7d+7o3Zmrylad9589ezZkMpm60bTUeWRDj2yrbj8yMhLJyclISUmpsaGp6aoD1XtuiWmGHfXI3ZwEgDMPkhSYAJBJ7OFyK0MlJSXhypUrmDlzJgDN8+NV49e3U9b1XHl5Oa5fv25QD4dCocDMmTMRHByMRYsWQQihd776qu+zoe+7KZ+Ptl4O1S11k5KS9JZX0+kCY5IQY3s8pPouOtJ32lF7Pci+2P83neySJXeWUo8FWLRoUbWBdjXFr+2c9r1793D//n0IIczqgm3fvj1WrVpVbWyCNtZMAFQq1111S92lS5fqLU+qqw7MTQCkupSvJvq+o9bonje114OnDqgyJgBkd6RKAFQ79C5dumicV658rbm2nb6qsb93755GAqBQKJCdnY3Ro0cjPj7e5C7Y8+fPq8cAyGQyya79twTVjIk1NTSRkZHYvHkzwsPDAehuaCx9lC1l17gpvUFSxyA1e46NrM8+9zpkFc5+NKCrsakpARBCoKioSOuAttTUVOTl5QGAugs2KysLkydPxrlz52qMSXWdu2p+Aikn/5GKKd+Lqpc4mtrQmPudlLJr3NRE1Brd86YmUjx1QBqEi1AoFAKAUCgUtg7FbvTp00c0btxY9OnTx6LbKSsrE3l5eaKsrEzveseOHRPdu3cXhw4dMmh9Qx06dEh06dJFHDp0SAghRElJiTh8+LAoKSnRGmtubq7IyckRubm56hhUy3fs2CFiYmI0YoyNjRX+/v6iW7du1WKuXPeysjIRExMjGjRoIOrVqyfWrVsnYmJixK5duySrqzZlZWXi0qVL4vbt26KioqLG9VXfi9jYWJGbm6vxPuhTUVEh7t+/LyoqKtSf5bFjx9QxGPKZmvudNHQ75pZl6nOWitVQttgmGc9a7ZX9HHaQ1VnraMDQozpDLxXTRdcRW9Xb7Bpy1K1rStoOHTpU6xJPSUlBgwYN8PDhwxrrmJycjJYtW2LFihX47LPPkJeXh7S0NKOO5Ew5MhWV5tiviakj9isflZp6jtqVjlAt0QPHiX7IGEwAXJgpO2lTdjCGdgebe6mYrtiq3mbXkARACKGRAKguc9PWzd+pUyd4enri1q1bmDZtmt73Jjw8HJs2bUKrVq3w9ttvIywsDMuXL7d4AmAMW15WZ0+X9Fn6FJklzseb+91gAuFamACQUbQNkKuJoUd1qoFkDx48wMsvvyzZjteYyWVUDf23336rvkkP8PvO+uLFi4iLi6t2rfupU6dQUlKCgIAA9Wx0KlVHqI8ePRrnz58H8PvVAIYMsDOXtimIa1pfNQthgwYN0KBBA7sao2Atlh4wZ4neDmOSFm2NPRMA1+J6v2oXJeUP25juZMC4ozohBN59912bDCAD/tjpL1y4UGNO+rS0NJSVlcHd3b1al/3s2bNx//59eHl5oX379hrlabs50eLFi02OzxRubm7Iy8vD6NGjDXpvKicAfn5+8PPzc5hr7KXcjr4GWorvmiV6O2pKWnjXQqqMCYCLcKQfe1JSkklHRoYesenbeaelpaFRo0YICAhAw4YN1acNIiMj8fnnn6N169bVuuwXLVqEJk2a4LXXXtM5Fa6q7LCwMMyaNcuoeknB1pd/OWICoK+BtvX7qUtNvQqONNkRWR6/BWSQygmE6lp4S50b1TbQzhCGdqnq23lHRETAy8sLDx48gLe3NyIiItT1rnxNfOUdaZcuXfDhhx+idevWererOsVRtZcAsHyCZu3Bdc7e0Oh7P21Zd3saQ0H2zzl/nSS5yg2UanpYWx396GosVTu/yo02UH3Qn76dt1KpxFtvvYWAgAA8evQIJ0+erPHeAKdPn8asWbNw8eJFg+KXyWSoW7eu+ty8Id2x5iYI1m4YnD0B0Pd+2rLuzv6+k7T4LXER5p6zVL1+69atePjwIRo1amTy0aS+G+4YMlitphv2AH/MZa+6nK9JkyZwd3cH8PtR/qeffoqIiAit5bdu3Rp169bF3bt3MXv27Brrk5SUhKtXr+LDDz/EmDFjanyP3dzc8NtvvyEuLg55eXlGTz3r7BM42QtHbEzNjZnfLdfiON9sMou55yxVr582bRpu3bqFunXrmnw0qa8Bz8vLw5gxY3TeVU4fXdP6FhUVQS6Xq5dnZmYiLi4OmZmZOsuaPn26+oY9qrJ1JSepqalo2LAhrl69itzc3BrfY5lMhqVLl5r8edjr+Wcp2FOja+lYbF1XbY29M3+3qDrb/8rIKrR1exuzA1K9fvny5QbfI94UNd1VTgpVJwaqTDWtLwCsWrUKADBq1CicPn1a53sVGRkJb29veHl5oaysTOM9rprsyGQyNGjQAMuWLUPz5s1N6kVx5slybN0oWpOt66qtsTfmu+VIA4tJB4vOM2hHXH0qYHOnAFW9vqSkxKjpYbVtu+o0sZUdO3ZMxMTEiPT0dJ3T1lYur/L0s5XLUJWvrd5VpwauTDWtb3R0tDh16pTo3bu3elrcu3fvirt371aLq6SkRKxfv15ER0eLXbt2idLSUnVMVacC1vb/mt4TQ+pUE23vE5nOGabU1fadM6ZeNa3L75zprNVeMQFwEebusFQ/5tLSUqPmla/8WtX6+uZ7NyTOmhrDyuVr28npuxfAoUOHRGRkpPjqq69Ebm6uOHToUI33JygtLRWnTp0SFy5cELm5uaKkpES9rq5Yq8Zd0xz4lZ93hsbH0RnTuNlrQ6jteyRlAsDvqel4LwCyK5W7K43tuqy6vrFd2CdOnEDv3r1x4sQJADWfp6xcvrZ19U0FrJo1MCYmBjKZzKDR8/rKq7x9fe9bTe+JIe+Zq3TJ2kM9jfkN2EO8llDTgEEOKHQAFk0v7Ah7AKTJxqUox9g7qVU9Oq6pB6DyMm09APqOyFS9Aw8ePFCf8sjLy1P3BGjroq+oqBC3b98Wv/76q0avwbFjx3R27VeN25ijKV3rusoRl6PVs6bTO7Zibg+AMb1WhrDXnhJb4CkAibl6AmDoj6um9aTY+erbhrby9d1atqbzmMbGqxof8NNPP2kkALGxsXp3ZpVvuVt5XUMbaykSgJoaGmfZwdprg6qLtW67bSxzEwBjxq0YwtE+V0tiAiAxV08ADGXr83rayq/acFVep/LOVbVeWVmZ1kF4hoiNjRUBAQEiNjZWo6wjR44Y3LhW3pHpanT11amm8nWVWVND42hHzrrYa4Oqi702bNq+R8YkiVLvKxztc7UkJgASYwJgGFsfRRpSvq4eAHOPaITQXX8pB0dJ9ZqqpD4is1eGXCliT5yl56WqmuplbL2d5fspBSYAEmMCoJ0xI/TtReWYzTnfr40U59dtdZleTdt1hM/WEGVlZWLXrl0iJiaGjYUTcdZEyRS8CoCsIjMzE0OGDFHPiucIk8zoGoGta8S/sZOtaBu17QyzwjnCZ2sINzc3pKWlccY6J2PriZFcEd9pF1e10XTku4lJ0cC5ublBLpejqKjI5ATAVjuymrYbFRWF77//HlFRUVaNS2pubm5Yvny5yTMpEtHvmAC4OEc/Kqzc6EnRwEnReFsjAag6N4Ih23WmIyxHT2YsPTeAs849QNJy/D0BmaXqEb+jNRLmTFCkiyNMYDJ9+nRcvHgR06dPt3UoNuFo39OqmACQPXDMXw9JpuqO1NF3rFJwhDuiOXrPjauzdJLpCEks2Z5MCCFsHYQ1FBUVwdfXFwqFAnK53NbhkB07fvw4EhISsHLlSkRHR9s6HK203eaYHEffvn3xyy+/oG3btti/f7/DlU+WZa32insOskvaznFbiyMMhLRkT83JkyfxzDPP4OTJk5KXTb+zdA8Oe4jIEOwBILtkyyMYRzm6PnHiBKZNm4bly5eja9eukpXbt29f5OTkoE2bNjx6tJDy8nJcvXoVTZs2hbu7u+TlO8p3mLRjDwC5NFsewTjKOAhLDQRcunQpWrZsiaVLl0paLv3BGeaVIMfHbwfZJUe/zMsaLJUkde3aFT/99JNRvQq2PGXjiNhAkz1wqG/f6tWr0axZM3h5eSEmJgbHjh2zdUhkIdxB1syekqTp06cjNzfXZS9LJHJEDrN3/eKLL5CYmIh58+bh5MmTiIiIwKBBg3Dr1i1bh0ZkE/aUJC1btgxt27bFsmXLbB0KERnIYQYBxsTEIDo6Gh9//DGA3we5BAcHY/LkyUhKSqrx9RwESEREjsBa7ZX0w08t4MmTJ8jMzERycrJ6mZubGwYMGIDDhw9rfU1paSlKS0vVfysUCgC/v7FERET2StVOqdotFU9PT3h6ekq2HYdIAO7cuYOKigo0btxYY3njxo3xyy+/aH1NSkoK5s+fX215cHCwRWIkIiKSUkhIiMbf8+bNw3vvvSdZ+Q6RAJgiOTkZiYmJ6r8LCwsRGhqK/Px8+Pr62jAy6RQVFSE4OBgFBQVOdVrDGevljHUCWC9H4ox1ApyzXgqFAiEhIfjtt99Qv3599XIpj/4BB0kA/P39UatWLdy8eVNj+c2bNxEYGKj1Nbq6Snx9fZ3mS6Iil8udrk6Ac9bLGesEsF6OxBnrBDhnverXr8+JgDw8PBAVFYX09HT1MqVSifT0dPTo0cOGkRERETkmh+gBAIDExETExcWha9eu6NatGz766CM8fPgQr776qq1DIyIicjgOkwCMHDkSt2/fxty5c3Hjxg107twZe/bsqTYwUBdPT0/MmzdP8nMotuSMdQKcs17OWCeA9XIkzlgnwDnrZa06Ocw8AERERCQdhxgDQERERNJiAkBEROSCmAAQERG5ICYARERELsghE4ADBw7ghRdeQJMmTSCTybBr1y6N58eNGweZTKbxGDx4cI3l2vp2wzXVq2qdVA9994N/7733qq3ftm1bC9fkDykpKYiOjoaPjw8aNWqEYcOGIScnR2Odx48fIz4+Hg0bNkS9evUwYsSIapM+VSWEwNy5cxEUFIQ6depgwIABuHjxoiWrolZTne7du4fJkyejTZs2qFOnDkJCQpCQkFBtXu+qTP3eSsWQz6pv377VYnzjjTf0lmvPn9WVK1d0/q62b9+us1xbf1Zr1qxBp06d1JPf9OjRA99//736eUf7TQH66+Sovymg5s/Kpr8p4YB2794t3nnnHbFjxw4BQOzcuVPj+bi4ODF48GBx/fp19ePevXt6y9y2bZvw8PAQ69evF+fOnROvvfaaqF+/vrh586YFa6KppnpVrs/169fF+vXrhUwmE5cuXdJZ5rx580SHDh00Xnf79m0L1+QPgwYNEhs2bBDZ2dni1KlT4k9/+pMICQkRxcXF6nXeeOMNERwcLNLT08WJEydE9+7dxdNPP6233NTUVOHr6yt27dolTp8+LYYOHSrCwsJESUmJpatUY53Onj0rhg8fLr755huRm5sr0tPTRatWrcSIESP0lmvK91ZKhnxWffr0Ea+99ppGjAqFQm+59vxZlZeXV/tdzZ8/X9SrV088ePBAZ7m2/qy++eYb8d1334lff/1V5OTkiNmzZ4vatWuL7OxsIYTj/aaE0F8nR/1NCVHzZ2XL35RDJgCV6UoAXnzxRaPK6datm4iPj1f/XVFRIZo0aSJSUlIkiNJ42upV1Ysvvij69eund5158+aJiIgI6QIz061btwQAkZGRIYQQorCwUNSuXVts375dvc6FCxcEAHH48GGtZSiVShEYGCjS0tLUywoLC4Wnp6fYunWrZSugRdU6afPll18KDw8PUVZWpnMdU763lqStXn369BFTpkwxuAxH/Kw6d+4sxo8fr7cce/ushBCiQYMG4rPPPnOK35SKqk7aOOJvSqVyvWz5m3LIUwCG2L9/Pxo1aoQ2bdrgzTffxN27d3Wuq7rd8IABA9TLarrdsK3dvHkT3333HSZMmFDjuhcvXkSTJk3QvHlzjB49Gvn5+VaIUDtVl52fnx8AIDMzE2VlZRrvfdu2bRESEqLzvc/Ly8ONGzc0XuPr64uYmBibfF5V66RrHblcDnd3/XNvGfO9tTRd9dq8eTP8/f0RHh6O5ORkPHr0SGcZjvZZZWZm4tSpUwb9ruzls6qoqMC2bdvw8OFD9OjRwyl+U1XrpI0j/qZ01ctWvymHmQnQGIMHD8bw4cMRFhaGS5cuYfbs2RgyZAgOHz6MWrVqVVvflNsN29qmTZvg4+OD4cOH610vJiYGGzduRJs2bXD9+nXMnz8fvXv3RnZ2Nnx8fKwU7e+USiWmTp2Knj17Ijw8HABw48YNeHh4aNzxCvj9vb9x44bWclTLtX1eul5jKdrqVNWdO3ewcOFCvP7663rLMvZ7a0m66vXXv/4VoaGhaNKkCc6cOYNZs2YhJycHO3bs0FqOo31W69atQ7t27fD000/rLcsePquzZ8+iR48eePz4MerVq4edO3eiffv2OHXqlMP+pnTVqSpH+03pq5ctf1NOmQC8/PLL6v937NgRnTp1QosWLbB//37079/fhpFJZ/369Rg9ejS8vLz0rjdkyBD1/zt16oSYmBiEhobiyy+/NOgoR0rx8fHIzs7Gzz//bNXtWlJNdSoqKsJzzz2H9u3b13gfb3v63uqqV+UdbseOHREUFIT+/fvj0qVLaNGihVVjNFZNn1VJSQm2bNmCd999t8ay7OGzatOmDU6dOgWFQoGvvvoKcXFxyMjIsMq2LUVXnSonAY74m9JXL1v+ppz2FEBlzZs3h7+/P3Jzc7U+b8rthm3pp59+Qk5ODiZOnGj0a+vXr4/WrVvrfC8sZdKkSfj222+xb98+NG3aVL08MDAQT548QWFhocb6+t571XJbf1666qTy4MEDDB48GD4+Pti5cydq165tVPk1fW8tpaZ6VRYTEwMAOmN0lM8KAL766is8evQIY8eONbp8W3xWHh4eaNmyJaKiopCSkoKIiAisWLHCoX9Tuuqk4qi/qZrqVZk1f1MukQBcvXoVd+/eRVBQkNbnHe12w+vWrUNUVBQiIiKMfm1xcTEuXbqk872QmhACkyZNws6dO/Gf//wHYWFhGs9HRUWhdu3aGu99Tk4O8vPzdb73YWFhCAwM1HhNUVERjh49apXPq6Y6qeIZOHAgPDw88M0339TYU6NNTd9bqRlSr6pOnToFADpjdITPSmXdunUYOnQoAgICjN6OtT8rbZRKJUpLSx3yN6WLqk6qeBztN6VL5XpVZdXflFFDBu3EgwcPRFZWlsjKyhIAxLJly0RWVpb47bffxIMHD8T06dPF4cOHRV5enti7d6/o0qWLaNWqlXj8+LG6jH79+olVq1ap/962bZvw9PQUGzduFOfPnxevv/66qF+/vrhx44Zd1EtFoVAIb29vsWbNGq1lVK3X22+/Lfbv3y/y8vLEwYMHxYABA4S/v7+4deuWxesjhBBvvvmm8PX1Ffv379e4zOXRo0fqdd544w0REhIi/vOf/4gTJ06IHj16iB49emiU06ZNG7Fjxw7136mpqaJ+/fri66+/FmfOnBEvvvii1S5ZqqlOCoVCxMTEiI4dO4rc3FyNdcrLy7XWydDvrS3rlZubKxYsWCBOnDgh8vLyxNdffy2aN28uYmNjNcpxpM9K5eLFi0Imk4nvv/9eazn29lklJSWJjIwMkZeXJ86cOSOSkpKETCYTP/74oxDC8X5TNdXJUX9TNdXL1r8ph0wA9u3bJwBUe8TFxYlHjx6JgQMHioCAAFG7dm0RGhoqXnvttWoNeWhoqJg3b57GslWrVomQkBDh4eEhunXrJo4cOWLFWumvl8r//M//iDp16ojCwkKtZVSt18iRI0VQUJDw8PAQTz31lBg5cqTIzc21cE3+oK0+AMSGDRvU65SUlIi33npLNGjQQHh7e4uXXnpJXL9+vVo5lV+jVCrFu+++Kxo3biw8PT1F//79RU5Ojl3USdfnCEDk5eVprZOh31tb1is/P1/ExsYKPz8/4enpKVq2bClmzJhR7ZplR/qsVJKTk0VwcLCoqKjQWY49fVbjx48XoaGhwsPDQwQEBIj+/furG38hHO83JYT+Ojnqb6qmetn6N8XbARMREbkglxgDQERERJqYABAREbkgJgBEREQuiAkAERGRC2ICQERE5IKYABAREbkgJgBEREQuiAkAERGRC2ICQERE5IKYABAREbkgJgBERBK5e/cuGjVqhCtXrhj1updffhkffvihZYIi0oEJALmccePGQSaTQSaToXbt2mjcuDGeffZZrF+/HkqlUrLt9O3bF1OnTtW7jioOXY/33ntPsnjsjSHvjyX98MMPNb7/P/74o1FlfvDBB3jxxRfRrFkz9bJx48Zh2LBhGut99dVX8PLyUjf6c+bMwQcffACFQmFutYgMxgSAXNLgwYNx/fp1XLlyBd9//z2eeeYZTJkyBc8//zzKy8utFsf169fVj48++ghyuVxj2fTp060Wi1SePHniENuLjY3VeK8bNmyId999V2NZ//79DS7v0aNHWLduHSZMmKB3vc8++wyjR4/GmjVr8PbbbwMAwsPD0aJFC/zzn/80qS5EJjHjLodEDikuLk68+OKL1Zanp6cLAOLTTz8VQvx+a+Xly5drrBMREaFxu+Xt27eL8PBw4eXlJfz8/ET//v1FcXGxiIuL03vbUm02bNggfH19tT5XUVEhFi1aJJo1aya8vLxEp06dxPbt29XP9+nTR0yaNElMmTJF1K9fXzRq1Eh88sknori4WIwbN07Uq1dPtGjRQuzevVuj3D59+oj4+HgRHx8v5HK5aNiwoZgzZ45QKpUGbbdyGVOmTBENGzYUffv2FUII8f3334uePXsKX19f4efnJ5577jn1raj1vT81ve+6tmdIrLpcvXpVANC4pa6xtm/fLgICAqotr/x9W7x4sfDy8tK4t7vK/PnzRa9evUzePpGx2ANA9H/69euHiIgI7Nixw6D1r1+/jlGjRmH8+PG4cOEC9u/fj+HDh0MIgRUrVqBHjx547bXX1EeTwcHBJseWkpKCf/zjH1i7di3OnTuHadOm4ZVXXkFGRoZ6nU2bNsHf3x/Hjh3D5MmT8eabb+Ivf/kLnn76aZw8eRIDBw7EmDFj8OjRI42yN23aBHd3dxw7dgwrVqzAsmXL8Nlnnxm8XVUZHh4eOHjwINauXQsAePjwIRITE3HixAmkp6fDzc0NL730EpRKpdnvj7btGRqrNllZWQCALl26GBxDVT/99BOioqJ0Pj9r1iwsXLgQ3377LV566aVqz3fr1g3Hjh1DaWmpyTEQGcXWGQiRtenqARBCiJEjR4p27doJIWo+Es3MzBQAxJUrV7SW1adPHzFlyhSD49LVA/D48WPh7e0tDh06pLF8woQJYtSoUeptVT56LC8vF3Xr1hVjxoxRL7t+/boAIA4fPqwRY7t27dRH/EIIMWvWLNGuXTuDtqsqIzIyssb63b59WwAQZ8+eVb9O2/tjSA9A1e0ZGqsu8+fPF8HBwTWu179/f/Hrr79qfe7FF18U48ePr7Y8Li5OeHh4CAAiPT1dZ9mnT5/W+30ikpq7TbMPIjsjhIBMJjNo3YiICPTv3x8dO3bEoEGDMHDgQPz5z39GgwYNJI0pNzcXjx49wrPPPqux/MmTJ4iMjFT/3alTJ/X/a9WqhYYNG6Jjx47qZY0bNwYA3Lp1S6Oc7t27a9S5R48e+PDDDw3eLgCtR74XL17E3LlzcfToUdy5c0c9wDI/Px/h4eEG1V2XqtszJlZtTp48adDR/8WLF9GiRQutz5WUlMDLy0vrc506dcKdO3cwb948dOvWDfXq1au2Tp06dQCgWg8NkaUwASCq5MKFCwgLCwMAuLm5QQih8XxZWZn6/7Vq1cK///1vHDp0CD/++CNWrVqFd955B0ePHlWXIYXi4mIAwHfffYennnpK4zlPT0/1/2vXrq3xnOoqh8p/AzD4SgdDtwsAdevWrfb6F154AaGhofj000/RpEkTKJVKhIeH1zhor6b3Xdv2jIlVm5MnT2LixInVlp87dw7jx49HSUkJRo4cicDAQLi5aT9z6u/vj/v372t97qmnnsJXX32FZ555BoMHD8b3338PHx8fjXXu3bsHAAgICKgxXiIpMAEg+j//+c9/cPbsWUybNg3A7zvi69evq58vKipCXl6exmtkMhl69uyJnj17Yu7cuQgNDcXOnTuRmJgIDw8PVFRUmB1X+/bt4enpifz8fPTp08fs8qo6evSoxt9HjhxBq1atzNru3bt3kZOTg08//RS9e/cGAPz8888a6+h6fwx536syJ9Y7d+6goKCgWg9AaWkpRo4cia1bt6Jjx44YNmyYRi9LVZGRkXpH8YeGhiIjI0OdBOzZs0cjCcjOzkbTpk3h7+9vVPxEpmICQC6ptLQUN27cQEVFBW7evIk9e/YgJSUFzz//PMaOHQvg90GBGzduxAsvvID69etj7ty5qFWrlrqMo0ePIj09HQMHDkSjRo1w9OhR3L59G+3atQMANGvWDEePHsWVK1dQr149+Pn56Tx61MfHxwfTp0/HtGnToFQq0atXLygUChw8eBByuRxxcXFmvRf5+flITEzE3/72N5w8eRKrVq3Chx9+aNZ2GzRogIYNG+KTTz5BUFAQ8vPzkZSUpLGOrvenpvdd6vfo5MmTAKoPANy1axf69OmjPo3Srl07NGnSRGc5gwYNQnJyMu7fv6/zNFBwcDD279+PZ555BoMGDcKePXsgl8sB/D6IcODAgXrrSSQlJgDkkvbs2YOgoCC4u7ujQYMGiIiIwMqVKxEXF6dupJOTk5GXl4fnn38evr6+WLhwocaRqFwux4EDB/DRRx+hqKgIoaGh+PDDDzFkyBAAwPTp0xEXF4f27dujpKQEeXl5GhPEGGPhwoUICAhASkoKLl++jPr166NLly6YPXu22e/F2LFjUVJSgm7duqFWrVqYMmUKXn/9dbO26+bmhm3btiEhIQHh4eFo06YNVq5cib59+6rX0fX+1PS+62JqrFlZWWjcuHG1xv3s2bPo3Lmz+u/MzEwMHjxYZzkdO3ZEly5d8OWXX+Jvf/ubzvWaNm2qkQT88MMP8PDwwK5du7Bnz54a60kkFZmoerKNiFxG37590blzZ3z00Ue2DsXuLFu2DL/99htWrFiBH3/8EYMHD8bdu3f1DvL87rvvMGPGDGRnZxvV27NmzRrs3LnT6JkHiczBHgAiIi1eeeUVDBkyBJGRkQgPD0doaGiNV3g899xzuHjxIv773/8aNa9B7dq1sWrVKnNDJjIKewCIXBh7AIhcFxMAIiIiF8SpgImIiFwQEwAiIiIXxASAiIjIBTEBICIickFMAIiIiFwQEwAiIiIXxASAiIjIBTEBICIickFMAIiIiFwQEwAiIiIX9P8B6z4dtA6R9hsAAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 512x384 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Polarization fraction as a function of the dust column density\n",
    "PvT_plot_Step2 = HtoH.PvT(CatP, CatT, Tscale=[15,35], Pscale=[0.0,30.0], CatdP=CatdP, errorbars='true', Source=FullName)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "491d0620",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Saving figures\n",
    "# PNG\n",
    "PvI_plot_Step2.savefig('Plots/'+SourceName+'_PvI_plot_Step2.png',dpi=300)\n",
    "PvN_plot_Step2.savefig('Plots/'+SourceName+'_PvN_plot_Step2.png',dpi=300)\n",
    "PvT_plot_Step2.savefig('Plots/'+SourceName+'_PvT_plot_Step2.png',dpi=300)\n"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "PoLiteWIP",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.10.8"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}
